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1.0 INTRODUCTION 


"Modeling and Analysis of Power Processing Systems," (Contract 
NAS3-21051), provides basic analytical tools. In the form of mathematical 
models, and the computer-aided techniques required to solve those models, 
to guide engineers In the various aspects of power processing equipment 
and system design. Early work In this area Included a feasibility 
study done by TRW Systems^^^ and General Electrlc^^^, and the Initial 
modeling and analysis performed by TRW Systems^^^ and California 
Institute of Technology^^^. The present effort, performed by TRW Systems 
and Virginia Polytechnic Institute and State University, Is an extension 
of that early work. 

The four tasks completed on this contract Include the following: 

e Performance Analysis of Buck, Boost, and Buck-Boost DC-DC 
Converters using dual-loop feedback control system^®^’^®^. 

e Design Optimization for Boost and Buck-Boost DC-DC Converters. 
(Design Optimization for Buck DC-DC Converters is contained in 
Ref""ence [3].) 

a Investigation of Current- Injected, Multiloop-Controlled Switching 
Regulators. 

a Application of Design Optimization to U.S. Am\y VSTOL Emergency 
Power System. 

A technical summary is presented in the following chapters on each of 
these tasks. 
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2.0 DC-OC CONVERTERS. 


2.1 Introduction . 

Th« following sections art both tutorial* and application oriented. 

Because of the variety of operating converter power and control schemes* 

the tutorial Is necessary to understand the various analytical procedures 

and their use. Once understood* the material may be applied to a designers 

specific needs through the various analysis and simulation subprograms 

which are provided. Starting with a description of the basic characteristics 

of OC-DC converters* a general evaluation of discreet time domain analysis* 

and the performance analysis of the buck* boost, & buck-boost converters Is given. 

Due to the finite flux capacity of the Inductive elements* a DC-DC 
converter must be oscillatory In nature. The oscillation Is achieved by 
cyclically operating the power switch of the converter In conduction and 
non-conduction state. Consequently, the converter control system must be 
able to accept an analog signal obtained from the sensing circuit and the 
reference, and to convert It Into discrete time Intervals In controlling 
the conduction and non -conduction of the power switch. 

The electrical performance of a DC-DC converter depends primarily on 
the quallt'' of Its control system. The performance characteristics of In- 
terest to a converter d'*s1gner Include stability as well as the converter- 
output response to step and sinusoidal disturbances, both from the line and 
the load. 


2.1 Introduction. (Cont.) 


Functionally, a OC-OC regulated converter can be divided Into two 
parts: A power circuit, and a control circuit. By definition, the power 
circuit handles the energy transfer from the source to the load. Three 
most commonly used power circuits are the buck, the boost, and the buck- 
boost. 

The control circuit manages the rate of the source-load energy 
transfer as a function of the load demands. During nominal steady-state 
and transient operations, the control objectives are associated with (A) 
the tracking of a certain controlled quantity In accordance with a given 
reference, and (B) the compliance to converter specifications such as the 
system response to step or sinusoidal line and load disturbances, and to 
external command signals. During transient operations, the control ob- 
jective Is to limit electrical -stress for all the elements associated with 
the converter, providing effective protection against catastrophic/degra- 
dation types of failures. A control circuit thus serves the multiple 
functions of regulation, command, and protection. 

A generalized standardized control module (SCM) has been developed 
to Implement the above control functions (Reference 5,6,7;. For the 
purpose of chls report, the SCM Control circuit has been selected to 
accommodate the power circuits mentioned above. 

A SCM-controlled DC-DC converter Is shown In Fig. 2.1. The power 

circuit, occupying the upper half of the block diagram, processes the 

transfer of energy from a raw Input to a regulated output V^. Three 

basic power stages are shown here: Buck, boost, and buck-boost. The 

control circuit regulates the rate of energy transfer. It receives an 
analog signal (V^) from the power-stage output, and delivers a discrete-time 
Interval signal (d) to achieve the required on-off control of the switch In 
the power stage. The discrete-time voltage or current pulses generated In 
the power stage are averaged by an LC filter having a much longer time 
constant than the discrete-time pulse Intervals. The averaged output therefore 
contains negligible switching-frequency components., and can be regarded as an 
analog signal containing only lower- frequency Information. 
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2.1 Introduction (Cont.) 


The SCN. occupying the lower half of the block diagram, performs the 
control function within the converter. It contains an Analog Signal Pro- 
cessor (ASP) and a Digital Signal Processor (DSP). Implementations of 
both ASP and DSP are standardlied: They combine to provide the required 
analog signal to discrete-t ine Interval conversion. The key feature of 
the SCM Is the utilization of an Inherent AC switching-frequency signal 
within the power stage. This utilization Is In addition to the con- 
ventional DC sensing of output V^. The sensed AC signal and the DC error 
are processed by the ASP. As a result, an adaptive stability Is obtained 
which Is Independent of the filter parameter changes. The SCM control 
function Is completed by the DSP, which processes the control -signal output 
from the ASP In conjunction with a prescribed control law, and operates 
the "ON-OFF" of the power switch via a duty-cycle signal, d. 

As stated previously, the control -circuit functions also Include 
command and protection. The command function generally requires the 
converter to respond to an external signal capable of overriding control 
signal (d) In determining the on-off of the power switch. The protection 
function Includes power-component peak-stress limiting and the converter 
shutdown in the event of a sensed abnormality such as overvoltage, under- 
voltage, o» ^vercurrent beyond a predetermined, tolerable level and duration. 
These functions are performed within SCM by the DSP^®*®*^^. 

The three basic functional blocks of an SCM-controlled converter are 
shown in the block diagram of Figure 2.2. 



DIGITAL SIGNAL 
PROCESSOR 


ANALOG SIGNAL 
PROCESSOR 


Fig. 2.2. switching regulator block diagram 


DC LOOP 



2.2 PTfofwno Analysis and Simulation Ttchnlouts . 


2.2.1 Nonlinear Ootratlon of Switching Rtaulator . 

A dc-dc switching-regulated converter Is Inherently a nonlinear de- 
vice. The first major nonlinearity exists In the power stage, and Is due 
to the operation of the power switch. Different circuit topologies 
correspond to the respective on and off time Intervals In the switching 
cycle. The second major nonlinearity exists In the Digital Signal Processor 
(DSP). Harmonic frequencies, which are multiples of the Input disturbance 
frequency, are contained In the DSP output. Because of such system non- 
linearities. difficulties are encountered In reaching performance assess- 
ments of various system performances such as stability, attenuation of 
sinusoldal/step line disturbances, and response to sinusoldal/step load 
disturbances. 

The Power Stage nonlinearity will be elaborated here, and various 
analytical approaches capable of treating this type of nonlinearity will be 
discussed. In light of the stated objective of this program, which Is the 
performance analysis and simulation of dc-dc switching regulators, a 
specific approach will then be selected as the basic analytical tool for 
the entire program. 


2.2.2 Power Stage Nonlinearity . 

Each of the power stages shown In Figure 2.1 can be divided as a 

function of the output filter Inductor MNF status. In Figure 2.3(A). 

often referred to as "continuous -conduction” or Node 1 operation, the 

NNF ascends during on time T^^ when the power switch In ON and the diode 

Is OFF, and descends during Tp^ when power switch Is OFF and the diode 

Is ON. Notice Chat the NNF never vanishes In the output Inductor. In 

Figure 2.3(B). often referred to as "discontinuous-conduction" or Mode 2 

operation, the MMF ascends from zero MMF at the beginning of T _, and 

on 

descends back to zero during Tp^ . An additional off time Tpg exists when 
both the power switch and the diode are OFF. during which the Inductor MMF 
remains zero, and load current Is supplied entirely by the output-filter 
capacitor. 
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2.2.2 Power Stage Nonlinearity - (Con't) 


Topologies of the buck» boost, and buck-boost power stages correspond 
to Tp^, and Tp 2 » are Illustrated in Figures 2.4 to 2.6, respectively. 
Even though a given power stage Is linear for each time Interval, the combi- 
nation of all different linear circuits for the purpose of analyzing a 
complete cycle of switching-regulator operation becomes a pi ecew1se-1 Inear 
nonlinear analysis problem. 


The difficulty Is Integrating these different topologies and collect- 
ively evaluating their responses to various line/load disturbances having 
a much longer time period than Individual T^^, Tp^, Tp 2 . 

The basic modeling approaches for conducting performance analysis 
and simulation include the following: 


e Discrete Time-Domain Analysis 
e Average Time-Domain Analysis 
a Exact Frequency-Domain Analysis 

• Discrete Time-Domain Simulation 


Small Signal Analysis 

.... Large Signal Analysis 


Techniques ^or all four approaches have been established In the 
previous MAPPS program phases^^ . Their respective utility for a given 
application depends on the analysis objective, the desired accuracy, the 
control-circuit type, the nature of the disturbance, and perhaps most 
Important, the user's analytical backgoround. A summary description of 
all four approaches is presented In Table 2-1. 


The exact frequency domain analysis will not be pursued In the pro- 
posed program due to its difficulty In Incorporating the Input filter which 
often causes major complications In the design for the required performances 
In regulator stability as well as In audlosusceptiblllty and output Impedance. 
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The performance requirements associated with the analysis and 
simulation effort are: 

e Regulator control-loop stability (local) 

e Audio susceptibility (attenuation of source small signal disturbance) 

e Transient response (response to source/load large-signal step change) 

For the first two performance categories, the nature of the disturbances 
Is such that the regulator can be regarded as a time-invariant system without 
a significant loss In analytical accuracy, and linearized about Its equilib- 
rium state to obtain a linear analytical model for small-signal performance 
evaluations. For the last category, the generally varying duty cycle sub- 
sequent to a step llne/load change represents a time-varying nonlinear system. 
Since any practical system Is likely to be higher than second-order, perfor- 
mance evaluation Is through tactics closely Identified with simulation 
techniques. 




2.2,3 The D1screte»T1me Domain Analysis. 


Realizing that a switching regulator Is Inherently a highly non- 
linear circuit containing analog-to-dlscrete-time conversion. It Is only 
natural that It can be more accurately analyzed through discrete time- 
domain modeling and analysis. Therefore, discrete time-domain analysis 
has been selected as the approach to be utilized for the performance, 
analysis and simulation program. 

In this approach, state-space techniques are employed to characterize 
regulators exactly through the formulation of nonlinear discrete time- 
domain equations In vector form. Newton's Iteration Is then used to solve 
for the equilibrium state of the regulator. The system Is then linearized 
about Its equilibrium state to arrive at a linear discrete time model. 

The closed-loop regulator Is thus modeled as a single entity rather than 
the three separated functional blocks shown In Figure 2.1. The stability 
Is studied by examining the eigenvalues of the linear system. The analysis 
can be extended, through Z-transform, to determine frequency-related 
performance characteristics such as audlosusceptiblllty. The modeling and 
analysis approach makes extensive use of the digital computer, making 
automation In regulator analysis possible. 

2.2.4 22 ete Time-Domain Analysis Objective . 

Discrete time-domain analysis Is the most accurate and straight for- 
ward of the different mathematical modeling techniques available. The user 
need only be proficient In the use of state-space analysis. Discrete time- 
domain analysis Is applicable to all types of power and control circuit 
configurations, operating in either continuous or discontinuous mode. Thus, 
It clearly stands as the best approach available. 

The performance analysis and simulation objective Is twofold: (1) the 

creation of generally applicable, practical, analysis subprograms and (2) a 
tutorial role of providing power processing designers with an effective 
analytical tool. 
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2.2.4 Discrete T1me»Doroa1n Analysis Objective - (con't) 

While time-dofflain analysis readily fulfills the second objective, 
the subprograms can only be applied by the prospective user if the user's 
regulator power and control circuitry is identical to the circuit confi- 
guration upon which the subprograms are based. Consequently, in pursuing 
time-domain analysis, the objective is to achieve the following: 

(1) To present clearly the tutorial information needed, by a user 
conducting time-domain analysis, for adopting the analysis to 
his specific application. 

(2) To create subprograms for regulators with standardized mul- 
tiple-loop feedback control . 
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2.2.5 Genera] Description of Performance Analysis and Simulation 
Techniques. 

Time-domain analysis may be applied to energy-storage converters to 
yield transient and steady-state solutions* stability, and audlosusceptiblllty. 
This approach has been applied here to buck* boost* and buck/boost converters 
with constant frequency and constant OFF time control. By way of Introduction* 
a general step-by-step procedure for performing time-domain analysis of energy 
storage converters will now be pres<^nted. 

Select the n X 1 system state variable vector x > [x] * 

Normally* the state variables are the voltage across the capacitors 
and currents through the Inductors. However, for the convenience 
of each Individual problem* state variables can be chosen differently. 

Write the system equations according to the modes of operation of the 
converter which are defined as follows: 

Mode 1 Operation ; The current through the Inductor Is always greater 

than zero. The period of each switching cycle can 
be clearly divided Into two time Intervals. Tqj^ and 
Tp^ . During the power transistor Is "ON" and 
the diode Is "OFF", and during Tp-j , the power tran- 
sistor Is "OFF" and the diode Is "ON". 

I ’ 

I * 


i 


Step 1: 


S. 1 2; 
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Mode 2 Optra t1 on ; The current through the Inductor reduces to zero 

and resides at zero for a time Interval Tp 2 > 

During this time Interval, both the transistor and 
the diode are "OFF". The time Intervals Tqj^ and 
Tp^ defined In the Model operation also exist In 
the Mode 2 operation. 


The system 

representation for 

the Mode 1 operation Is 


• 

% ■ 

FI X ♦ 61 u 

during T^,^ 

(2-1) 

• 

X ■ 

F2 X ♦ 62 u 

during Tp^ 

(2-2) 


where 


Ir 

X 


L 

m 


dx, 

3T 


dx 

3T 


and 


_3T 


The column vector ^ Is a (mxl) Input vector, containing the 
Input voltage Ep the reference E|^, the saturation voltage 
drop across the power transistor, and the forward voltage 
drop across the diode. The nxn matrices FI and F2 and the 
nxm matrices 61 and 62 are constant matrices containing 
various circuit parameters. 
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In th« Nodt 2 op«ratlon, aquation (2-3), Must b# addad to tho 
syttan roprosontatlon: 

i • F3 X ♦ 63 u during (2-3) 

Tho dimensions of F3 and 63 art not ntctssarlly the same as those of 
FI and 61. respectively. 

The converters, which are basically nonlinear switching circuits, 
are accurately described by the plecewUe-llnear representations 

(2-1) to (2-3). 

Step 3 The general solution of the linear differential equation 


X - F^ X 4 6^ u 1 • 1, 2, 3 (2-4) 


Is 


x(t,^4T) . e^(T)x(t^) 4 o,(T) u (2-5) 

where e^(T) • e^^^ 1 • 1, 2, 3 

0.(T) • e^<^ [ / ^ e’^^^ds] 6. 1 - 1. 2, 3 

0 

The tenns e^(T) end D^(T) can be computed either analytically or numerically. 
If they are computed numerically, the following Taylor series 
expansion may be used: 



- 1 4 F^T 4 


CF,T)2 

"TT'* 


(F^T)^ 


1 ■ 1, 2, 3 


Step 4 Write the discrete-time-domain equation for the converter. The 
discrete- time-domain equation can be expressed as 

x(t,j 4 ^) • 0 x(t,j) 4 V u (2-6) 


where t|^ end tj^^^ correspond to Instants of time at the 
beginning of the It^ cycle and the k4l ^ cycle, respectively. 


Mode 1 Operation 

It may be shovm, for Mode 1 operation, that the terms ^ and V 
of equation (2-6) are given by - 




V • ♦^(TpJ) DKTqJ) ♦ 02(Tpjf) 


(2-7) 


( 2 - 8 ) 


It Ic 

where T^jJ and TpJ represent the and 

Tp^ Intervals during the fc th cycle. The time 
k k 

Intervals Tq^ and Tp^ can be determined through 
the following two conditions: 

Condition 1 A threshold condition, which Is determined 

by the particular type of digital control 

signal processor en^l'^yed In the converter, 
and may be expressed as - 

t^(x(tj^). TqI;. Tplf) • 0 (2-9) 


Condit ion 2 A condition which specifies whether the 

converter Is operating at a constant frequency, 
or a constant ON time, or a constant OFF time, 
or a constant voltage- second, and Is expressed as - 


«2<^o5- ® 


( 2 - 10 ) 


Mode 2 Operation 


For mode 2 operation, the terms ^ and V of eq. (2-6) are given by - 


♦ ■ •a^^Fi^ ♦i^^oN^ 

V ajdpk) a^dpH) D,(T^) 




(2-11) 


(2-W) 
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In ordtr to dotemine tho Hint Inttrvalt 

^o5’ ^F1 ^2* • third condition, In 

addition to (2-9) and (2-10), should bt 

Included to detect the time Instant when 

the Inductor current reduced to zero. 

This condition may be written - 

tj(*,(t|^)» Tqi^, Tp^ , Tp2) * 0 (2-13) 


Of course, In the Mode 2 operation, the time 

Interval 1^2 should also be a parameter 1n (2-9) 
and (2-10). Thur , (or Mode 2 operation. 


C,(»(t|,). 4. Tf,. rfj) . 0 .ni • 0 


Equations (2-6 to 2-13) are the complete rep 
resentatlon of the converters. 


Step 5 Solve for the approximate steady state x.* 

The approximate solution Is employed later as an Initial guess 
toward solving the exact steady state thi*ough Newton's Iteration 
method. In the steady state, equation (2-6) may be written as 

X* • ♦ X* ♦ V u (2-14) 

The a matrix and V matrices can be computed for the given Tq^^, 

Tpi , and Tp2* f^or given Input-output requirements of the converter, 
the approximate time Intervals, Tq|^, Tp^ and Tp2 can be determined. 
Reference 1 gives a detailed list of duty cycle formulae for differ- 
ent power stages and different control schemes. 
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If the (l-f) it non-singular, equation (2-14) may 

be solved for x* 

X* • (I • V u (2-15) 

However, in many cases, the matrix (I -a) Is singular. In ordei’ 
to solve for x*, equation (2-6) together with (2-9) is required. 


Step 6 Solve for the exact steady state x*. 


Newton's iteration method is employed to find the steady- 
state solution with the initial guess Equations (2-6) through 

(2-13) are used in the iteration process until a specified state- 
•"atching condition is satisfied. The state-matching condition can 
be defined as 





< c 


(2-16) 


for an arbitrarily small positive number, or it can be defined 
simn''y as 

l*i(Vi) - < c (2-17) 


Step 7 Linearize the discrei time-df/main equation, Eq. (2-6), about 
the equilibrium state x_* for studying stability, audiosuscepti- 
bility, and transient response to a small step change of the 
input or the load. 


Equation (2-6) may be written in the form 

i^Vl^ - f U^. Ify Tp2) (2-18) 

where it the Input voltage. 
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Form the term 6 x (t,^ ^ 1 ) by taking partial derivatives of 
Eq. (2-18) with respect to x and and writing 

* i * li *i <‘k> ♦ Is;- *“i 

• ♦ 6 jK (tj^) ♦ r (2-19) 

where << Is a (nxn) matrix and r Is a (nxl) column matrix. 


The differentiation of (2-18) can be performed analytically. 

If the problem Is simple, or numerically by difference quotients. 


Step 8 Analyze the stability of the converter. 

The linearized system (2-19), Is stable If and only If all the eigen- 
values of the matrix >p are absolutely less than unity, I.e., 

jx^ ) < 1 1 « 1 n 

The c envalues are evaluated by the computer. Changes of eigen- 
values as a function of system parameters can be plotted In the 
complex plane, the x-plane. The location of the eigenvalues in 
the x-plane Indicates not only the stability but also the transient 
behavior of the system, I.e., damping and rapidity of response. 

Step 9 Analyze the Audio-Susceptibility of the converter. 

Audiosuscept1b1l1ty may be defined as the frequency response of 
the output voltage to a small amplitude sinusoidal perturbation 
of the Input voltage u^ . The Z-transformatlon may be used with 
Equation (2-19) to derive the audlosusceptibll Ity as a frequency 
domain transfer function. 
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( 2 - 20 ) 




The output voltage may be exprestod as 

where C Is a constant (Ixn) row matrix. 

The Z transformation of (2-19) Is 

a X (z) - (iz - s-)"^ rau^(z) (2-zi) 

The frequency-domain transfer function can be derived after 

replacing Z with e’^'^^P (2-21X and combining (2-20) and (2-21). 

WJ-) * i^TOSTTTT • 4 - ♦> ' 

where Ej and E^ are the dc average of the Input voltage and the 
output voltage, respectively. 

Step 10 Study the transient response of the Converter. 

Thp linearized system remain: valid for a small step change 

of Input voltage or load, since the system still continues 

to operate about Its equilibrium state. The behavior of the 

transient with respect to damping, oscillatory nature, decay 

time, and overshoot of the equilibrium position Is governed 
by the location of the eigenvalues on the X-plane. 



2.2.6 Performance Analysis and Simulation Objectives 

Since the principles and procedures for performing the time-domain 
analysis are the same regardless of the circuit configuration used, both 
objectives stated In the previous section can be achieved by conducting 
the time-domain analysis on regulators using SCM multiple-loop control. 
Certain Performance Analysis Subprograms (PAS's) using this control were 
generated In the Initial MAPPS phase II effort under NAS3-19690^^^. A 
summary on what has been done Is presented In Table 2-II. 

In Table 2-II, the three basic power stages, buck, boost, and 
buc’.-boost, are separated Into Mode I and Mode II operations with 
constant on time, constant volt-second, constant off time, and con- 
stant frequency control. To cover all possible categories, there 
are a total of twenty-four different power/control configurations. 

For the particular multiple-loop control configurations marked by 
"x", time-domain analysis has been applied and completed through 
previous contracts. Those marked by "NA" are configurations Incom- 
patible with the control Implementation, and therefore are not 
recommended for hardware design or time-domain analysis. 

Multip1<^-loop control senses the rectangular voltage across the 
Inductor, and integrates it to form a triangular ramp output. The 
triangular ramp possesses a negative slope during on time Tj^ and a 
positive slope during off time Tp. Since the regulator control 
determines the point at which the ramp intersects the fixed threshold 
level, the following conclusions become apparent: 

( a ) In constant on time or constant volt-second control, regu- 
la tion is achieved by controlling off time T p. The threshold 
level therefore prescribes the peak of the triangular ramp 
as the intersection of the ascending ramp with the threshold 
level marks the end of off-time Interval Tp. 



guratfon Not Reo 






2.2.6 Performance Analysis and Simulation Objectives (Cont.) 

( b ) In constant off time or constant frequency control , regu - 
lation Is achieved by controlling the on time Tqj^. The 
threshold level therefore prescribes the valley of the 
triangular ramp, as the Intersection of the descending 
ramp with the threshold level marks the end of the on-time 
Interval 

This rule governing the ramp-threshold Interface Is further Illustrated 
In Figure 2.7 for continuous 1nductor-MMF operations. Here, the Inductor 
voltage and the integrator output voltage arc- shown for constant-on-tIme 
control. Notice the difference in the relative position of the threshold 
level with respect to the triangular ramp. 

When the inductor MMF becomes discontinuous. In Figure 2.7(B), 
the Inductor voltage vanishes for a certain Interval. During this 
Interval, the integrator output exhibits a flat top. In the constant- 
on-time control, this flat top coincides with the threshold level. Any 
slight noise in either the flat top or the threshold level Is likely 
to trigger the on-off control, making the constant-on-tIme control high- 
ly suscept’ 0 to noise when engaging In dIscontInous-MMF operation. 
Conversely, in the constant-off-time control, the flat top caused by 
the zero inductor voltage level is above the threshold level, thereby 
eliminating this noise susceptibility problem. Since the MMF Is dis- 
continuous when the regulator output load Is light, the constant-on-time 
and the constant- vol t-second control, shown In Figure 2-II, are con- 
sidered to have limited utility, and therefore are not to be included 
in the analytical effort. 

Hence, this effort concentrates on the time-domain analysis of 
multiple-loop control based on constant off time and constant frequency. 
The analysis includes all three power stages, with both continuous and 
discontinuous MMF operations. 
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Figure 2.7 Thre$ho1d-R«ip Interface for (A) Continuous- 
If!F and (8) Discontinuous WiF 







2.2.6 P erformance Analysis and Simulation Objectives (Cont.) 


Perform Discrete Time-Domain Analysis 

In performing the discrete t1me*doma1n analysis, three subprograms, 
one for each of the three basic power stages, are create... Each subprogram 
includes both continuous and dIscontInuous-MMF operations as well as the 
constant off time and constant frequency control schemes. Each subprogram 
(buck, boost, and buck-boost) has the capability of: 

(1) Either constant-off-time or cons tant-frequenc)^ control 
In one composite subprogram, or 

(2) Either continuous- or dIscontInuous-MMF operation In the 
same subprogram. 

In dIscontInuous-MMF operation, the existence of a third time Inter- 
val during which the inductor MMF vanishes complicates the composite 
subprogram to a certain extent. Two separate computer subroutines are needed 
to compute the exact equilibrium state of the system for both continuous and 
discontinuous modes. The Information Is then fed Into a common linearization 
subroutine to numerically derive the linearized system for small-signal anal- 
ysis . 

An information flow chart Is presented In Figure 2.8 to show how the 
converter, , wO both continuous and discontinuous current operations together 
with two types of duty-cycle controllers, can be Implemented In a single 
computer program. The flow chart presented here Is sel f-expl an' ''ry, therefore 
no description will be given. 
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2 . 3 Formulation of State Equotions 

Mathematical models have been generated for the Buck, Boost, 
and Buck-Boost DC-DC convertlons using the constant frequency and 
constant off-time pu1se-w1dth>modu1at1on techniques, and are 
discussed In detail In the following sections. 

2.3.1 Buck DC-DC Converter 

2. 3. 1.1 Constant Frequency Buck Regulator 

The buck dc to dc power converter topology Is shown In Figure 
2.9. The analysis approach Is based on the ensuing mathematical model 
and the capability to consider both continuous and discontinuous In- 
ductor current operation. 
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t 



1 


Fig. 2-10 Uavtfonn of for Discontinuous Inductor- 
Current Operation. 

EQUIVALENT DISCRETE TIME SYSTEM 


The Maveform of e^ vs. time shown In Fig. 2.10 is used to establish 
some notation regarding the time instant t|^ when each cycle starts and 
each switching action occurs: 

In steady state operation, 


At t^. 
At t|[. 
At t^, 


the clock pulse turns the power twitch 

♦I 

k4l 

t| • t| the threshold condition turns the power 

' ' switch 'OFF' 

k^1 k^2 

• t« the zero Inductor current condition turns 

off the power diode. 


The time Intervals t{ -tp and hre defined as tJ^. 

Tp^ and Tp 2 » respectively. I^se time intervals may vary from cycle to 
cycle. However, the time Interval between t|^ and Is a constant equal 
to the period of oscillation Tp, i.e., 




(2-23) 
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•nd. 


Sj • O3 


0 

‘‘d 

*?r. 


The vectors X and U art state variable and forcing function vectors, 
respectively. 


X t 


U t 




w m 



•0 

*2 

• 

1 

*3 


•c 

L J 

» dlk 

“l 


h 

“2 


h 


The constraints of the system are governed by the threshold 
conditl jn at t » t^i , t2 . . . . 

*3 “il > ■ 't- 

It k^l 

the zero Inductor current condition at t • t2 • • • • 

Xg ,(tj ) • 0, 


(2-27) 


(2-28) 


and the constant frequency condition . . . 


^ for .11 k. 


(2-29) 


Each of thf lintar equations (2-24) to (2-26) a<ftaU« a cloaad fom lolution of 
the form 


KtJ) • ^ li 

(2-30) 

X(tJ) ■ X(tJeTj^) • #2(T^,) X (tj) ♦ 02(T^,) 0 

(2-31 ) 

X(V^) • X(tJeT^2^ • ^ ^ 03(Tp2) H 

(2-32) 

where ei(T) ■ e^^^ 1 • 1* 2, 3 

(2-33) 

01 (T) • ] 61 1 • 1. 2, 3 

(2-34) 


The structuras of the matrices ai and D1 for 1 ■ 1» 2. 3 have the following 
forms: 


el 


01 


01 


♦M) 

el 12 

**2) 

el 22 

-**31 

el 32 

«11 

0 

^*31 

0 

■dlj, 

d1 


0 0 

0 0 

0 d1 


for 1 « 1 . 2 . 3 
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Using equation (2-30) In equation (2-31), and the result In equation (2-32) 
the equivalent discrete time system for the constant frequency buck converter 

Is given by: 

X(t^^,) • e3(T^2^ •2(T^,) al(Tj^) X (t^) ♦ #3(T^2) ♦2(Tj‘,) OI(tJ^) U 

♦ #3(T^2^ ®2(T^,) U ♦ 03(TjS2) U (2-35) 

Also required In this description are the threshold conditions derived 
from (2-27) 

*1^H^ * *^32^^W^ *2^^k^ * *^33^^0N^ *3^^k^ * ^^31 *^1 

♦ dlj2(T5,^) U2 • £y. (2-36) 

the zero Inductor current condition, derived from (2-28) 

♦22i(Tpi) X^(t,^j) 4 ♦ 222 ( 1 ^ 1 ) Xgttjjj) • 0. (2-37) 

and the constant frequency condition 

Tp2 ■ Tp-4-Tn • <2-38) 
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Equilibrium Solutions 

In steady>state. the following conditions prevail: 


• 1* ■ constant 
for all k 


-k^l 

'on 

rk4l 


■ T 


ON 




constant 


■ Tp^ ■ • constant 


'FI 

pk^l 

T2 




const^int ■ Tp-Tgj^-T ^1 


The Approximate Steady State 


(2-39) 


(2-40) 

(2-41) 

(2-42) 


The approximate T^. T^^ and It can be computed using the following 
formulae: 


where 


^5n 


•/S 


0^0 


_ . Wo 

T?2 • '^p‘^8n“^F1 

Lq ■ the energy storage Inductance 
F^ ■ output power 
E| ■ Input voltage 
Eg * output voltage 


(2-43) 

(2-44) 

(2-42) 
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.J 


Substituting (2-39 to 2-42) Into (2-35), one obtains 


1* • O(Tfj) «a(T^,) ♦l(TJ^) 1* * {♦JCTfj) «(Tf,) 01 (TJ^) 

♦ .J(T|j) 02(T^,) ♦ D3(Tfj) } 11 

• * i* * i *^8 n' J1 

where (TJh- T?,. T^) . »3(Tfj) rtUf,) and V (TJir T?,. 

. *3(Tf2) *2(T^,) DUT{„) * *3(T^2> OZiTf,) * D3(T^2) ■ 

The magnitudes of Tg^j, , and •''*® obtained from equations (2-42) to (2-44) 


Equation (2-45) may now be solved for the state vector X*. Using this solution 
we may now solve for the states X*. X*. and X*. Solution begins by arbitrarily 

setting X| = 0. Then: 




♦ (♦3n*^ll^’ll^''^n*^12‘*^21 * *^12*^21^^l'^*^12*^22^^21^ 


♦ (e3^^d2^^ ♦ e3^2^^21^ 


4 d3„E, 


Equation (2-46) may be solved for Xif 




<a3^^[42^^d1^^4*2^2<^l2^4d2^^] ♦ ®3^2^*^21^^1^*^22^^21 


4 d22^] * d3,, ) E 


n ' -I 


(2-46 


(2-47; 
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(2-48) 


Using this result* the state may now be derived from equation (2-36): 

*3 * - ♦'31**1 - *'32*1 ■ <"31^1 ■ "'32 h 

In this approximation, the threi:ho1d condition where the Inductor- 
current X* equals zero may not be satisfied. This approximation Is merely 
employed as a starting point In order to search for the exact steady state. 

The Exact Steady State 

Define the system state with the power switch off as 


I(t|.) 4 X (2-«) 

4 i (V^0N*^F1* 

where, using Fig. 2-32, it is clear that 

Y 3 (t|^) » Ej for all k (2-51) 

Z2(t,^) « 0 (2-52) 

In steady-state operation 

Y* = ♦1(T5,^)X* + D1(T*^)U = fl(T*^, X*) (2-53) 

Z* = «2(Tf^)Y* ♦ D2(T*^)U - f2(T*^ . Y*) (2-54) 

and 

X* « ♦3(T*2)Z* + D3(T*^) U 


It Is tmportant to note that and T|^ are functions of X* and u via the 
threshold conditions (2-51) and (2-52). 

If Tgi^ and T^^ are the exact steady state values, then the steady-state 
calculated from (2.46 to 2-48) has to satisfy the following two matching 

conditions: 
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(1) the zero*1nductor*currtnt condition 

Wi** »2 • « 

(2) the state matching condition 

* ♦*32*^F2* *2 * H 

♦ dSjjCTfj) Uj - XJ . 0 

Newton's method may now be used ',o find the Tg^i and which satisfy 

the matching conditions. 

The step-by-step procedure is described as follows: 

Step 1 Use the approximate given In (2-43 & 2-44) to 

the approximate state X* from (2-46) to (2-48). 

Step 2 Find a new Tpi by Newton's method such that for the 
given X* and ^ON together with the new Tp^ , the zero- 
current condition (Tqj^, X*. ) = 0 will be 

satisfied. 

T* ■ T* - ^«Mtch^^ON^* i* 

[aB /aT 1 • 

Step 3 Check if S„,.^. = 0 is satisfied. 

^ match 


(2-56) 


(2-57) 


derive 
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step 4 


If ^ satisfied, modify Tgj^ according to 

the equation 



T* ^match^^ON* 


Step 5 Use the new and Tfr calculated in Step 2 and Step 4, 

to derive a new approximate state X*. Then to go to Step 2 
and repeat tho process until the state matching condition, 

^n,atch ’ 


A flow diagram for determining the steady state is presented in 
Fig. 2.11(A) and (B). A subroutine is developed to search for 

a proper Tp^ to satisfy the zero-inductor condition shown in (2-56) 


This subroutine is contained in another subroutine S . . which 

match 

ultimately computes the state matching condition given in (2-57) 
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r 

approximate 


v_ 

X* 



COMPUTE 
• S«.tch<TjN> 


“ • *«tch<T8N*‘To«) 
OS . (SE-$t)/4Tj^ 


^ON ’ 


COMPUTE 

"^match 


|S1| < c 

'tT 




Z2 

USE APPROXIMATE 
STEADY STATE 


Continut with the 
remainder of the program 

FIGURE 2.11A nOH DIMMAN FOR OETERNININS THE STEROY STATE 
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Analysis of Linearized Discrete Time System 


The analysis of stability* audio susceptibility* ond transient response 
due to step change In the Input voltage and the load* Is presented. The 
analysis Is based on a discrete system linearized about Its tgulllbrium 
state. 

Derivation of the Linearized System 

The linearized system can be derived by perturbing the system at the 
cycle. After the perturbation the nonlinear discrete time system vquatlon 
(2t 35) can be rewritten as: 


»h*p* { (X , U). .3 (tfj) «2 (rf,) »1 « 

♦ .3 (T^j) 01 (tJ^) U 

♦ .3 (tfj) 02 (t‘,) U 

♦ 03 (t{^) U (2-55) 

This system can be linearized about Its equilibrium state. 

If the following two terms are defined: 


«JL “ i* • 

and 

6‘Ji (t,^) • U, (t^) - 
It follows that: 

«x (t,^ ^ ♦ «x (t^) ♦ Vj, (t^) (2-56) 




- 45 - 


where ♦ 


(2-57) 


and r • «) y* (2-58) 

The matrix ♦ Is 3x3 and the matrix r Is 3x1. 

The partial derivatives may be approximated by difference quotients, 
and evaluated numerically. 


For sufficiently small a x^,(1 ■ 1 ^)* 

f,(«/nx,) - . . . f,(X3«*A«3) - f,(X 3 «) 


« 'x*. U* 


ax, 


ax, 


f,(x/*ax,) ■ f3(x^*) . . . f^ (x3*»ax3) - ^3(>^3*j 


ax, 




(2-59) 


since X does not eppe.r explicitly In f. the ch.nge of T . Tp,. end 
“ u Af 1 ■ 1 . .» 3, must be determined first 

doe to e '•'•"S' ^ ,, .wording to the 

In evaluating (2-59). The new «nu ip^ 

threshold conditions (2-36) and (2-37). 

Similarly, 

r 

f,(U,+aU,) - f,(U,) 

‘W, 


il 


(2-60) 


1 'x*. y* 


fj(U^+aU^) - ^3 ((^t) 


“aU 


X*, y* 
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It Is Important to select the appropriate Increments axj and aU^. 

Some experimentation with the Increment size Is advisable* since the accuracy 
of the partial derivatives depends on It. If the linearized system shows 
high sensitivity, that Is, changes Its behavior rather rapidly as It moves 
away from equilibrium, then the results obtained for the linearized system 
are valid only for very small perturbations about the equilibrium state. 


The Stability of the Linearized System 
The linearized system 

♦ «x (t^) ♦ r6U^(t^), 

Is stable If and only If all the eigenvalues of are absolutely less than 
unity, i.e., 

|X^I <1 1 - 1, 2, 3, 4 (2-61 

The eigenvalues are evaluated by the computer. Changes of eigenvalues 
as a function of system parameters can be plotted In the complex plane. The 
location of the eigenvalues In the complex plane Indicates not only the 
stability but '’'>0 the transient behavior of the system, I.e., damping and 
rapidity of response. 
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2. 3. 1.2 Buck Converter with Constant Off-T<me Duty Cycle Control 

The basic structure* and analysis approach, for the buck PAS program with 
constant off- time duty cycle control are the same as for the constant frequency 
buck PAS program. The differences between the two converter schemes exist In 
the procedures for the computations of the converter switching times. 

Constant off-time, Tp, as defined In the context of previous converter 
analyses, Is the total time that the switching transistor remains off during 
a switching cycle. Therefore, 

Tp 5 Tp^ + Tp 2 r constant (2-62) 

Arbitrarily, Tp has been assigned the value of the sum of Tp^ and Tp 2 » 
determined for the constant frequency buck converter operating In MODE 1. 


The expression of conservation of power for the buck converter Is 
given as: 


(2-63) 


Substituting Tp » Tqj^ + Tp in the above equation results In the following 


quadratic equation for T^^: 


hn ■ 2 L. P. - 2 P^ T, 


■0 *^0 'BN ■ " •-0 *^0 'F 
Solving this quadratic equation for T^j^ gives the following expression: 

, . '•0 >■„ C-o '•0 'o Tf) 

E (E -E ) 

Where the sign has been chosen before the square root term since 


(2-64) 


(2-65) 


•-o ''o ^ Ljj Pq Tp) 


( 2 - 66 ) 
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As In the constant frequency buck control, the clock pulse Initiates the 
period. Following the computation of T^|^, the other switching times 
can be readily computed: 


FI 


(2-67) 

F2 

■ V-Tpi 

(2-68) 

P 

" * ^F 

(2-69) 


The sequence of testing for the duty cycle scheme and Inductor MMF mode of 
operation Is Illustrated In the computational flow chart presented In 
Figure 2.12. As In the constant frequency buck PAS program, the error cri- 
terion for MODE 2 operation Is that Tp 2 Is greater than EPS ■ l.E-6. 

The buck PAS program (Appendix A, Volume II) Is written such 
that one computer program package may be used to analyze both duty cycle 
control schemes operating In either continuous or discontinuous Inductor 
MMF mode. 


















2.3.2 Boost DC- DC Converter 


2. 3. 2.1 Boost Converter with Constant Frequency Duty Cycle Control 

Given the boost circuit of Figure 2.13, and irie state and Input vectors 
defined above, the state equations can be readily determined. To facilitate 
the development of the state equations, two duirmy variables, the voltage e^ 
and the current 1 q, as defined In Figure 2.13, are Introduced. The resulting 
generalized state equations for all modes of operation within each switching 
cycle can be expressed In terms of the state and dummy variable and Inputs. 


"c " “ 

wv 

»<L 

" * ''c ^ (Hl + R 5 ) • 


(2-70) 

^ “ ■ 


. ^ . e^ + ^ . Ej 

0 0 


(2-71) 

®R ‘ - 

C 2 R 5 * ‘ 

1 

*R * C 2 R^ * ’ ''c 

1 ( "s“l 

I- <2-72) 

®c “ 

1 

. ^l'^5 

(^rj. 

*c * ^ 

( 1757 )- 'R 

+ 

1 

^iRg 

1 / \ 
^1*^3 1 *^ 1 '*’ *^ 2 /J 

'd ♦ (cfiq)- ‘1 



^ ‘ * * ^R (2-73) 


Similarly, the output voltage, e^, can be written: 

®o " • ■‘d ^ • ''c (2-74) 


The general form of the state equations Is: 
X . ^ X + ^ u 


where i =1,2 and 3 refer to the times T„, Tp , and Tp , respectively. 

on r 1 In 
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During the power switch S.| Is on and the power diode $2 

Is off. Therefore, the dunn^y variables are set to the following 
values: 


«( ■ Eq 
<0 ■» 


Consequently, FI and G1 can be expressed es follows 
FI • 


where 


G1 


fl,. 



0 

0 0 


, 

0 

"22 

0 0 



fij, 

0 

"33 ° 


(2-75) 

_E'41 

^^42 

"43 1 



91 2, 
0 

0 

0 

0 

0 0 

9^23 ° 

0 0 


(2-76) 

£'41 

9142 

9>43 ^ 



1 


fi 

1 s 



flj2 • 

Eo 

'31 

^*^5 ^s/ 

1 

^2*^R 

E'41 ■ 





nR 


^^3 “ ^ 9^21 “ “ 9^23 “ 

9^1 • - 9I42 - 9^3 “ 
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The output voltage becomes: 


" TTTX 


. V. 


During the power switch Is off and the power diode $2 
Therefore, e^ and Ip are: 


«1 1?^ • ''c 


. 1 + E, 


^0 “ ^ 


The resulting F2 and 62 are: 



F2 - 


f2i2 

0 

0 



f22i 

f222 

0 

0 



^^31 

^^32 

f2 

33 0 



_«41 

f242 

f2 

43 ° _ 


62 > 

"o 

0 

0 

0 1 



922, 

0 

0 

CSJ 

CM 

cn 



0 

0 

0 

0 



9*41 

9*42 

0 

9244 

where 

1 


f*,j . 

1 

*L 


rSjT- 


TR[nqy 



*’s*L ^ 

i 1 

f*3. 

- ^ (^ 

- - V^O - 

'C ' 

"07 VI 


Is on. 





42 ^1^4 [^1^3 
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The output voltage e^ is now 


e 


0 



. i 


( 


During Tp^» the inductor current has gone to zero and both the 
power switch S-| and the power diode $2 are off. The dummy variables 
take on the following values. 


e 

i 


i 

D 



0 


The F3 and 63 matrices discribing this mode of circuit operation are 


f3i^ 

0 

0 

0 

0 

0 

0 

0 

^^31 

0 

^^33 

0 

^^41 

0 

^^43 

0 
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G3 


0 

0 

0 

0 



0 

0 

0 

0 


0 

0 

0 

0 


( 2 - 
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where 


"n ■ - C,(n’* *s) i ”31 • ^ 

*h^ ■ - ipq ‘ «41 • - ■ (iq^ - 

^^43 " rjV^ ’ 9^42 " 

The output voltage, is the same as during 



By transforming the state equations into the state transistion format (Ref. ) 
the preceeding equation sets may be represented in standard form: 

Ik.l • Jl***"!!! 

During each time duration within a switching cycle, the following 
state transition and control distribution matrices are generated 



«1 - 

♦111 

0 

0 

0 

D1 * 

0 
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0 

♦122 

0 
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dl2i 
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dl23 
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♦^43 

J 


^41 

<"42 

«’43 

0 


For Tp , 


■ 

CM 

0 

«11 
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0 

<1224 
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0 
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♦^41 

♦242 
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0 
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0 

D 3 - 
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0 
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0 
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0 

«4, 

0 

”43 
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0 
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where 


♦i ■ e 
D1 - e 
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1 - 1. 2. 3 


By combining the above state transition matrices, the nonlinear discrete 
state transition equation for a complete switching cycle Is obtained. 


A = *X (t,^) + V u 

where tj^ and t|^^^ correspond to the time Instances at the beginning of 
the kth cycle and the (k-t-Uth cycle, respectively. The « and V matrics have 
the following structures: 


♦11 

♦12 ° 

0 

V « 

Vll 

0 

Vl3 

*14 

^21 

♦22 ° 

0 


V21 

0 

^23 

*24 

♦31 

♦32 ^33 

0 


V31 

0 

V33 

*34 

*41 

^42 *43 

1 


V41 

V42 

^43 

*44 


From the structure of 0, It Is apparent that the next values of X^, X2 
and X3 are Independent of past values of X^, a property which Is useful In 
determining the steady state solution. 
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Analysis 

In order to approximate the steady-state solution necessary for the 
linearized system analysis* the approximate switching times TFl, and TF2 
must be computed. The following formulae may be used: 


Ton - 


¥ 


'FI 




'F2 


Tp * Tqn 


'FI 


( 2 - 86 ) 

(2-87) 

( 2 - 88 ) 


where 


Lq ■ the energy storage Inductance 
pQ ■ output power 
n * efficiency 
Ej • Input voltage 
Eq “ output voltage 

Eq saturation voltage drop of the power transistor 
Eq ■ forward voltage drop of the diode 


The expressions for Tq|^ and TF1 are generated based on the flux and energy 
conservation principles given below (refer to Figure 2.13): 

Flux Conservation 

T0N(E1 - EQ) ■ TF1(E0 ♦ ED - El) (2-89) 

Energy Conservation 

nCl(£l • , EQ) ton? ♦ rflUO ♦ K • £1) TFl^ ■ ( 2 _ 90 ) 

‘■'p • P 

For the computation of the approximate steady-state solution, the following 
state transition equation Is solved: 

X* - e X* 4 V U 


- 58 - 


Therefore: 


X* - [I - e]-^ y y 

By imposing the property of the state transition matrix, e, noted above, the 
state variables X<|. X 2 and X 3 can be computed Independently of X^. Therefore, 
the solution of the steady- station expression may be partitioned such that only a 
3x3 matrix need be Inverted to solve for the X^ * ^3 states. The existing 

voltage restriction on X^ (tj^) [capacitor voltage In Figure 2.13] Is used to 
solve for X^. The following expressions are the closed form solutions for the 
approximate steady-state vector: 


XI 

X2 

X3 

X4 


tO-^22) ^ ♦ fi2 VU23/DEN 
[♦21 0-^u) VU 23 /OEN MODE • 1 

0 MODE • 2 


(2-91) 

(2-92) 


^<♦21 ^32 ♦ai 0"^22)) ♦ (^32 ♦ ^12 ♦si) ^2)3/DET 


♦ WO-^33) 


(2-93) 


ET - ♦l^j XI - ♦l^j X2 - ♦l^j X3 - d\j U1 - d\j 

- dl,,j U3 - dl,,,, U4 The voltage restriction (2-94) 

on X4(t|^) [e^^ (Figure 2.13)] 


Inhere 

yy ■ yy 

DEN ■ (1-4u)0-422) - 4i 2 ^21 
OET • 0-#8j) OEN 

MODE ■ I 1 continuous Inductor current operation 

I 2 discontinuous Inductor current operation 

The performance of PAS2 (Appendix B) , Volume II, with these computation 
procedures was checked out. 
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Continuous Inductor Currtnt Mode of Opention 

In ordtr to make It complete in its capability of analyzing 
converters, PAS2 was modified to Incorporate the continuous Inductor 
current mode of operation. The appropriate subroutines requiring modifica- 
tions were: 

1 ) PAS2 

2) STATE2 

3) PSIM2 

4 ) FFUNC2 

5) 6AMM2 

Linearized System Analysis 

The eigenvalues of the Jacovian matrix characterize the linearized 
system stability and indicate the system transient behavior, i.e., damping 
and response time. 

The system eigenvalues, as functions of the circuit parameters, are 
computed numerically. In order to characterize the stability boundaries 
(if relevant) of the linearized system, these eigenvalues are parametically 
plotted on the Z-plana and analyzed with respect to stability and response 
criterion. 

Eigenvalues near the positive real axis of the Z-plane indicate that 
the system has a long time constant; eigenvalues equal to zero in discon- 
tinuous current mode indicate that the inductor current state variable 
sampled at -jy cycle is insensitive to the state variable in the previous 
cycle. The indictor current translates to zero following any small signal 
perturbation. 
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2. 3. 2. 2 Boost Converter with Constant Off-Time Duty Cycle Control 

The primary structure; and analysis approacn,of the boost MS program for 
constant off*t1me duty cycle control are the same as for the constant frequency 
boost PAS program except for differences In the evaluation of the converter 
switching times. 

Off-time. Tp, as defined In the context of previous converter 

analyses, Is the total time that the switching transformer remains off during 
a switching cycle. Therefore: 

♦ ^F2 (2-95) 

In constant off-time control, Tp Is held constant while Tp^ and Tp 2 are allowed 
to vary. Arbitrarily, Tp has been assigned the value determined for the constant 
frequency ‘'?iost converter operating 1n MODE 2. 


The conservation of power for the boost converter may be represented by 
the equation: 


.2 _ 2 T, (to ‘ to - Ej) 
tH rrn'i-' 


(2-96) 


Substituting the expression Tp ■ T^^^ + Tp In the above equation results 
In the following quadratic equation for T-„: 

n Ej (E, - tpXto ♦ to - E,) T,* - 2 L, P, (E, ♦ Ej - E,) Tgo (2-97) 

- E 1. ^ «0 * 'O ■ 'l> Tf - 0 

Therefore, solving this quadratic equation for T^^^ gives the following solution: 

Too * ” * 
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where 


QA ■ n E, (E, . Eg)(Ejj ♦ Ep - Eg) 

^ (£p ♦ Ep . E,) 

QC • QB • Tp 

and where the positive sign has been chosen before the square root since 
QB < /(Qb2) ♦ 4 • QA • QC 

As in constant frequency boost control, the clock pulse Initiates the 
period. Following the computation of the other switching times can be 
readily computed: 



»i ■ 'o> 

• Eo * Ep -"j, 

(2-99) 


“ ^P * ‘ “^Fl 

(2-100) 


“ ^F 

(2-101) 


The testing sequence for the constant off-time duty cycle scheme and the 
different Indu^ wor MMF modes of operation are Illustrated In the computational 
flow chart presented In Figure 2-14. As In the revised constant frequency 
boost converter PAS program, the error criterion for MODE 2 operation Is that 
Tp2 must be greater than 1* of the switching period. 

The boost PAS2 program (Appendix B, Volume II) is written such that one 
computer program package may be used to analyze both duty cycle control schemes 
operating in either continuous or discontinuous Inductor MMF mode. 
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2.3.3 Buck Boost DC- DC Converter 


2 . 3 . 3 . 1 Buck Boost C onverter with Constant Frequency Duty Cycle Control 

Figure 2>15 preset . he circuit configuration of the buck-boost 
converter with a constant frequency duty cycle scheme. 

As opposed to the Inductor current state variable representation utilized 
In the previous BUCK and BOOST power converter analysis programs, the In- 
ductor flux, which Is continuous, replaces the discontinuous Inductor current 
as a system state variable. Therefore, the state and Input vectors, X and U 
are defined as: 



( 2 - 102 ) 



(2-103) 
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-30 S 




The constraints on this system are: 


1) • E, 

(2-104) 

i) XjU*) • riux(tj) • 0 

(2-105) 

3) . T, 

(2-106) 


In order to utilize the Inductor flux as • state variable, the magnetiza- 
tion characteristics of the Inductor have been quantified. The flux Is ex- 
pressed by the following equation: 

Flux . 

where w s permeability or the ratio of the flux density to the magnetizing 
force 

A = cross-sectional area 
i s length 
N s number of turns 
1 £ current 
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For convcnienct In deriving the stete equations, a constant k Is defined: 
k E ^ 

Therefore, the expression for Flux may be written: 

Flux • kNI - ijf- 

The following expression for k can then be utilized for the resulting equation 
set: 


where k Is In units of webers per amp*turns. 

In the derivation of the state equations. It Is convenient to define the 
current i In the following manner: 

k’’ 

1 • Kp Flux , where Kp -||j- 

Ourlng the period, when the switching transistor Is on and the diode Is 
off, the following set of equations describes the power converter operation: 


• 

"c 





(2-107) 

Flux 

*■0 ”l 




(2-108) 

• 


w 1 



(2-109) 

®R 





• 

. .ri 4,^ 

^2 1 

A. V ♦ 

*'‘3*^0 

Flux 

®c 






* CTC •* * 0 : Sr • '1 * El 

1 

[Rj^R|k2 ) 

•h 



. J "3 J 
♦t-SJ- RJ-'q 




(2-110) 
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(2-ni) 


'L^S 




During the period when the switching transistor Is off and the power 
diode Is on, the following equation set describes the system: 


'c • - *c * s 


(2-112) 


Flux 


- ♦ V q (2-113) 


1 1 1 
spi; “c ■* tpi; ' • ips^ «R (2-114) 


where 


/I ^2 1 ^3 

• <E^ * E^TK^Ifjtti-lRi+Rjn' • • *c 

r 1 ^2 1 ^3 1 ^3 T 

" fjTR^RpfRJDC^^R^TT ' rp^ Epi^ SJ 'olS-f’"* 


1 


1 


1 ^3 


R, RjR, 

* ''c 


(2-115) 


(2-116) 


k-1 

KJ" 
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During the Tp2 pe»*1od when both the switching transistor and the power 
diode ere off, the following state equations describe the converter system 
operation. 


Flux 




(2-117) 


(2-118) 




(2-119) 


«c * 








l^'^c 


(2-120) 


( 2 - 121 ) 


In matrix representation, the above state equations for each time period 
may be expressed In the general form: 


x»FlX+^U i* 1,2,3 
where for Tqj^ 


flu 

0 

0 

0“ 

. G1 • 

”0 

0 

0 

0“ 

0 

^22 

0 

0 


•’21 

0 

«>J3 

0 

fl31 

0 

^33 

0 


0 

0 

0 

0 

_fl41 

«42 

^43 

0 


^> 4 , 

»U2 

*’«3 

0 
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«nd for T 



The state transition equations are of the form: 


*k+l ■ ♦ Oi u 1 • 1.2.3 
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and 

♦3 • 

♦3,, 

0 

0 

o“ 

03 • 

n 

0 

0 

o“ 



0 

0 

0 

0 


0 

0 

0 

0 



♦^31 

0 

♦^33 

0 


0 

0 

0 

0 



_43^1 

0 

♦^43 

1 


0 

«42 

0 
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Switching Times Analysis 

The approximate steady-state switching times may be derived for the constant 
frequency control case. The boundary condition inherent with this system Is that the 
time Tp Is a constant, that Is: 

Tp ■ Tqi^ ^ ^ ^F2 * constant (2-122) 

Due to the conservation of flux In the Inductor, the following restriction 
applies: 



Finally, an expression Involving the conservation of power Is needed to complete 
the necessary Information for calculating the approximate times: 



E, (Ej-Eq) 

tut; — 



(2-124) 


Using these expressions. It Is possible to derive the following relationships: 


/ “ OTp ii; — 

^OH • V n E, it{- 

N, t,- E. 

■ Ilf ^ON 


(2-125) 

(2-126) 


Tp2 ■ Tp - Tq^ - Tp, 


(2-127) 


As In the previous material. Mode 1 Is the condition that Tp 2 " 0 
(Inductor current Is continuous), and Mode 2, the condition that Tp 2 ^0 
(Inductor current Is discontinuous). Assume that the circuit now operates 
with a load resistance such that Tp2 - 0. Allow the load to vary uniformly 
In a direction where Inductor current at the end of approaches zero. 
That value of load resistance at which the Inductor current at the end of 
Tp^ first becomes zero Is a critical resistance. (Equivalently, the value 
of resistance that first gives a non-zero Tp 2 Is a critical resistance). 
This critical resistance may be found In the following way: 


From the power conservation expression, we have: 

* n £j (E, - £q) 

From this expression, and using the expression for Tp.j, we have: 


* ^F2 


) 


T2 


" El(E,-£o> *L ^ON - sf <V 


2LV 


^»1 T 

5 — Q — I * ^ON 


(2-128) 


(2-129) 


The critical load resistance may now be found: 

N 




critical R. 


"'i (Ei-V^on 


(2-130) 


For MODE 2 , the switching times are given above. For MODE 1, 
howeve*^, the following times are more easily calculated using the equations; 






Tpi • Tp - 


(2-131) 

(2-132) 


Tp 2 • 0 (by definition) 


(2-133) 


For the actual PAS program, the criterion for MODE 2 operation Is given by 
the Inequality constraint: 

Tp2 > 0.01 Tp 

that Is, Tp 2 must be greater than 1% of the switching period. 
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Ste«dy»Stite Solution 

Propagating the state variables through a complete switching cycle 
results In the following expression: 

*(^♦1^ * li li IL til 51 il 51 fill H (2-134) 

i e x(t,^) 4 V u 


Since i* 

may be calculated using the following procedure: 


X* - [I - ♦]"’ V u 
where x* • [x* x^ xp^ end 


Ij - 4141 - fl 


'42 *2 


♦^43 *3 * “1 


dl42 “2 ' **^43 “3 


The exact steady-state solution for the state variables can be computed using 
a Newton- Raphson Iteration algorithm to find the proper Tq|^ end Tpi* and hence. 

x^, which satisfy the state matching conditions defined below. The matching 

k 

condition on the Integrator output voltage, and on the Inductor flux. 


“match- 


^match 


*4 


- xj • #34^ Jy 4 #343 *3 ♦ *4 ♦ ^^42 U2 • X4 


'41 ‘1 


'43 ‘3 


Ctch • *2 • *221 >1 ♦ ♦^2^ >2 ^ ®22$ U4 


Iterating on Tq|^ and Tp^ In order to drive these matching conditions to 
xtro results In the proper exact steady-state solution for x. 

Therefore, by defining the following variables 


jf(tj^) X 4 Tqj^) ■ #1 1^ ^ 51|j M. 

i<‘k> ‘ > • 42k 4<‘k’ ♦ 22* a 

£(tk„) .i(t,*T4).i2^KV4 0J^u 
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the state variable boundary conditions nay be expressed as: 



• 4141 *^(t,^) * 4142 *2^H^ 


♦ ♦ dl^2 U 2 + dl^3 Uj ■ E| 

(2-135 


• *hi J'l<‘k> * *hz l' 2 <‘k> * "^24 “4 ■ ® 

(2-136 

T ^ • 
^F2 

T T ^ T ^ 

‘p ' 'on 'fi 

(2-137 
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2.3. 3. 2 Buck-Boost Converter with Constant Off-TIma Duty Cycle Control 

Tht basic structure end analysis approach of the buck-boost PAS program 
for constant off-time duty cycle control Is the same as for the constant frequency 
buck-boost PAS program. The differences between the two converter schemes 
exist In the procedures for the computations of the converter switching times. 

Constant off-time, as defined previously, Is the total time that 
the swltchln^^ transistor remains off during a switching cycle. Therefore 

Tp 5Tp^ ♦ Tp 2 s constant (2-138) 

Arbitrarily, Tp has been assigned the value of the sum of T^^ and Tp 2 
determined for the constant frequency buck-boost converter operating In 
MODE 1. 


Conservation of power for the buck-boost converter Is expressed 
In the equation: 



^ ^0 V >’e 


(2-139) 


Substituting for Tp In the above equation results In the following 
quadratic equation for T^^^: 


^ V * ^ *'0 % V ■ ^ *'0 ^0 ° 


(2-140) 


Solving this quadatric equation for T 0 nj gives: 


. . S ^ ♦ ■> '•0 

Where the "+" sign has been chosen before the square root because 

>•0 ‘V “-O ♦ " Ei(Ei-EqK 2 Tp) 


(2-141) 


(2-142) 
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As In the constant frequency buck-boost control, the clock pulse Initiates 
the period. Following the computation of the other switching times 
can be readily computed: 


Tp, ■ Hiilal . T, 


0N 


T„ • ^F-Tpl 




(2-143) 

(2-144) 

(2-145) 


The sequence of testing for the duty cycle scheme and Inductor MMF mode of 
operation are Illustrated In the computational flow chart presented In 
Figure 2-16. As In the constant frequency buck-boost PAS program, the thres 
hold criterion for MODE 2 operation Is that Tp 2 Is greater than U of the 
switching period. 


The buck-boost PAS program (Appendix C, Volume II) Is written 
such that one computer program package may be used to anlayze both duty 
cycle control schemes operating In either continuous or discontinuous 
Inductor MMF mode. 
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2.4 Computer Programs . 

Computer programs have been generated using the three mathematical 
models described. Figure 2-17 shows the general Information flow chart 
for a composite subprogram Including both continuous and discontinuous 
Inductor current operation, and constant frequency or constant off-time 
control . 

The performance analysis computer flow diagrams are Included as 
part of Appendix A through C, Volume II, for the buck, boost, and 
buck-boost regulators, respectively. Although the programs are 
similar, that for the buck-boost Is most complex. In review of the 
presently available data. It seems reasonable that the buck-boost 
computer program can be modified so that It Is capable of performing 
the buck, boost and buck-boost DC-DC converter performance analysis. 



F1gurt?.-17 Information Flow Chart on a Composite Subprogram 



2 . 5 Computer Performance Analysis Results . 


The performance analysis computer programs provide the following 
results : 

• Linearized stability analysis. 

• Root locus analysis - up to 10 optional system parameters. 

§ Audio-susceptibility analysis. 

• Transient response analysis - linearized system. 

• Discontinuous or continuous inductor MMF mode of operation. 

• Load change analysis - linearized system. 

• Analysis of a non-linear transition response to a step input voltage. 

• Constant frequency or constant off-time duty cycle control . 

The computer programs contained in Appendix A through C (Volume II) 
have been verified experimentally ror the buck, boost and buck-boost 
regulators, respectively. 

A brief example of the computer results for a buck-boost converter is 
presented to illustrate the type of data available. 

Figures 2-18 to 2-21 show the root locus results of a buck-boost 
DC-DC converter where circuit parameters are varied in order to establish 
the relative ..ability of the converter. 

• Figure 2-18 C2 varied 

• Figure 2-19 R 5 varied 

• Figure 2-20 R4 varied 

§ Figure 2-21 CO varied 

Figure 2-22 and 2-23 show the small signal open-loop gain and phase 
relation-i.lip. Figure 2-24 shows the output transient response when the 
input voltage is changed from 24 VDC to 40 VDC. 

Figure 2-25 shows the output transient response when the output load 
resistance is changed from 49 ohms to 600 ohms. 
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FIGURE 2-13 EIGENVALUES AS A RIHCTIOH OF THE COMPENSATION CAPACITOR Cj 








FIGURE 2-19 EIGENVALUES AS A FUNCTION OF THE COMPENSATION LOOP RESISTOR R5 
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FIGURE 2-20 EIGENVALUES AS A FUNCTION OF THE AC-LOOP RESISTOR R* 
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FIGURE 2-21 EIGENVALUES AS A FUNCTION OF THE OUTPUT FILTER CAPACITOR Cg 























3.0 DESIGN OPTIMIZATION OF POWER CONVERTER 


3.1 Introduction. 

In converter design practice, it is often attempted to find the small- 
est possible magnetic core to accommodate the necessary windings in order 
to satisfy a given set of design requirements [g]. It is hoped, in return, 
that the weight and size of the overall converter would be reduced. One 
effective way to reduce the weight and size of the magnetic components 
(a way which many designers are pursuing nowadays) is to Increase the switch- 
ing frequency of the converter. To a certain extent, this approach is viable. 
Equations which govern the selections of design variables (such as magnetic 
components and capacitors), however, are nonlinear and interdependent due 
to the complex nature of various functions of power converter design. Employ- 
ing the conventional design approach, only a piecemeal, ^uboptlmum converter 
design at best can be achieved. For example, when the size of the energy 
storage inductor is reduced, the ac switching current component is invari- 
ably Increased. Consequently, a penalty is imposed on the weight and size 
of the l""->ut filter design in order to attenuate the ac switching current 
component which reflects back into the source. Similarly, as the ac 
switching current component is increased, larger output filter capaci- 
tors should be used to limit the output ripple voltage component. To give 
further illustration of the complex interrelations: When the switching 

frequency is increased beyond a certain range, the gain of weight- and 
size-saving of magnetic components is diminished because the magnetic 
core losses and the semiconductor switching losses are increased 
as a function of switching frequency. Thus, higher losses and heavier 
overall system design can result due to the Increase of weight ai»d 
size of package and heatsink. The goal of a minimum weight converter 
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design Is seldom achieved, despite the extensive end tlae-consuming trial 
and error design process. 

In power-converter design, the key to implementing a successful design 
optimization rests largely on the availability of suitable mathematical 
and computer programming technique Handicapped by a lack of suitable design and 
optimization tools, the tendency has been for a designer to rely on 
time-consuming intuitive and empirical methods resulting in a sub- 
optimum design. Such inadequacies invariably lead to penalties involving » 

equipment weight, operating efficiency or other performances. 

The philosophy of design optimization (s) of a power converter is 
briefly described in the following: 

3.1.1 Design Optimization 

A non-optimum design (Illustrated in Fig. 3.1.1) generally involves 

four design Ingredients First, a set of performance requirements 

such as output ripple factor and frequency-dependent source conducted 

EMI level, r ■■ (r, , r„, ..., r ), is given to guide the design. 

1 z m 

Second, a set of design constants, such as transistor switching times 
and rnaxlmum magnetic operating flux density, k * (k^, k^, ..., t 

is employed. These constants are known to a designer either through ^ 

manufacturers’ data sheet, or designer’s common sense and experience. 

Third, the objective of the design is to pinpoint numerically all 

« 

the unknown design variables such as the detailed magnetic core 
size, mean magnetic path length and other component sizes, 

X » (x^, x^, x^) • These three design Ingredients are then 

integrated together with the fourth ingredient, the nonlinear design 
constraints such as outvut ripple voltage constraint, EMI constraint. 
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DESIGN APPROACH WITHOUT OPTIMIZATION 



n > J --> INFINITE SET OF SOLUTIONS TO SATISFY 

ALL CONSTRAINTS qjlx.k.f) =0 
AND PERFORMANCE REQUIREMENTS 


Fig. 3.1*1 Conventional design approach 
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Due to the design complexities of the switching power converter^ the 
number of unknown design variables generally exceeds that of the constraints. 
Therefore there exists a virtually Infinite set of solutions which satisfy 
both the constraints and performance requirements. Different designers 
may come out with completely different designs based on Individual design 
experience and approacl^es. For example: Having selected a power circuit 

configuration, the designer picks a power-converter switching frequency 
through intuition or experience; then, using the performance requirements, 
he proceeds to obtain input filter design, output filter design, and energy 
storage inductor design. The same procedure may be repeated several times 
for different switching frequencies before a comparatively lower weight 
power-circuit design is achieved. Such a design is characterized by the 
designer’s subjective judgement. Despite the time consuming iterations, 
optimization of the overall power converter is seldom achieved. 

The optimum design approach takes advantage of powerful, computer- 

aided, nonlinear programming that integrates all the design Ingredients 

in the objective iunctlon - that is, the function to be optimized - defined 

by the designer. The objective function can be the converter weight, 

efficiency, or any other realizable physical quantity. The philosophy 

of design optimization (shown in Fig. 3.1.2) is obtained through the 

injection of an objective function f(x,k) into the ' onal design 

approach. The essence of the design optimization is when to realize a set 

of design variables using nonlinear programmln dque. This set will 

satisfy all constraints gj •0 and requirements , and concurrently 

optimize a certain converter characteristics, f(x,k), defined by the aesigner. By 

introducing the objective function Into the design 'irocess, the infinite set 

of design solutions presented in Fig. 3.1. *s reduced to a single solution 
set that Is an optimum of the objective function. 
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BY INTROD'jr.ING THE OBJE CTIVE FUNCT ION INTO THE DESIGN 
PROCESS > SINGLE SET OPTIMUM SOLUTION (Xj^/ X2^.../Xp) 
CAN BE PINPOINTED TO SATISFY ALL CONSTRAINTS Qj = 0 A 
CONCURRENTLY OPTIMIZE THE OBJECTIVE FUNCTION f(xJO 


Fig. 3.1.2 Op*:iroum design approach 
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3.1.2 Power Converter OptimlBetlon 

The piirpoee of this report la to demonstrate the uaefulneas of an 
Augmented Lagranglan Multiplier (ALAG) based nonlinear programming 
technique. This la used for a minimum weight design of the boost and 

a 

buck-booat power converters. Previous experiences of the authors with 
the use of AI.AG based programming techniques for power converter design 
optlmlxatlons [15] have encouraged such an approach. According to * 

expectations, reliable results have also been obtained for boost and 
buck-boost convertets. 

At the beginning, mathematical models are presented for the boost 
and buck-boost converters. Various design requirements and physical operating 
characteristics of these converters are summarized in the form of equality 
and inequality constraints. The minimum weight design requirement Is 
formulated as the objective function. 

A. Prob lem Form ulati on 

The circuit schematic, the objective function, and set of constraints 
are briefly discussed here for the boost and the buck-boost convertors 
shown in Fig. 3.1,3. » 

A .l. Objective F unction ; The objective function is formulated as a sum 

of various component weights which include: » 

(I) Core weight. 

(II) Winding weight. 

(ill) Capacitor weight. 
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(Iv) Source weight. 

(v) Package aii'' heat sink weights. 

A. 2 Constraint are ; A numbiar of equality and Inequality constraints 
form the constraint set: 

(I) The loss constraint which is composed of input filter copper 
loss, conduction and switching losses of transistor and diode, two-wlndlng- 
inductor copper and core loss, and output filter capacitor BSR loss. 

(II) Operating flux density constraint. 

(ill) Window area constraint. 

(iv) Parasitic resistance constraint. 

(v) Input filter peaking constraint. 

(vl) Frequency dependent source EMI constraint. 

B. Solution Methodology 

The stringent requirements for modeling the power converter design 
give rise to a set of very complicated nonlinear equations and an objective 
function. Obvl Misly, such a model does not lend Itself to a closed form 
solution; but one may use numerical techniques to arrive at an optimum 
solution. There are several nonlinear programming algorithms which provide 
convergence from a reasonable set of Initial guesses. The selection of such 
an algorithm depends on the characteristics of the problem at hand, and the 
availability or non-availt-blllty of a feasible starting solution. 

In the course of this research project two nonlinear programming 
algorithms were found to be appropriate for use in the minimum weight 
design optimization study. Both of these algorithms are based on 
transforming a constrained optimization problem into a sequence of 
unconstrained problems. The successive solutions of unconstrained 
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problama convarg* to a aolution of the constrained problem. The 
Sequential Unconstrained Mlnlmlsatlor. Technique T) [10,U] end 
the Augmented Legrenglen Multiplier Technique (ALAG) [l/!, 13, 14] are 
t«K> popular peremetrlc traneformatlon technlquea that have bean examined 

* during this study. The ALAG algorithm has been found to be the faster and 
easier technique for the problem at hand [l5]. The dealgn optimisation 

* results that are presented In this report are obtained using the ALAG 
algorithm. 

The following optimization approach Is proposed for finding the 
converter minimum weight design: 

(I) Fix the switching frequency. 

(II) Find all the circuit parameters which give the minimum 
system weight for the given frequency. 

(ill) Change the system frequency over a certain desirable range 
and repeat the process. 

C. A dv a nt ages of Computer Aided Design (CAD) 

Using tne Computer-Alded-Deslgn (CAD) approach, designers no longer 
have to relay on subjective and brute-force trial and error methods. 
Computer-aided design not only provides the optimum solution but also 
offers the following advantages: 

(1) The CAD approach is cos. effective, since the switching frequency, 
circuit components and optimum masaetlc designs, (down to the 
details of core size, mean magnetic path length, eve.) can be 
obtained In one computer run. This capability has a unique 
distinction over the conventional piecemeal suboptimum design. The CAD 
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approach Intagrataa tha intardapandant natura of tha varioua 
functlona of tha powar convartar. 

(2) By traating tha awltchlng fraquancy aa a paramatrlc constant 
in tha almulatlon procaaa, tha ovarwhalming computation tima 
and convergance difficultiaa which otherwiaa raault can ba 
reduced to a minimum. 

(3) Aaaeasment of tradeoffs between converter weight and loss as 
function of switching frequency is immediately possible through 
the proposed approach. 
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3.2 


BOOST CONVERTER OPTIMIZATION - OBJECTIVE FUNCTION AND CONSTRAINTS 


In thla chapter, the formulation of objective functions and design 
constraints for the boost power converter are presented: 

(I) To use this practical example in order to Illustrate the 
power converter design optimization using the nonlinear programming 
techniques; and (2) To demrnstrate the minimum weight design of the 
switching power converter and its weight/loss tradeoffs. The circuit 
schematic with a two-stagic input filter is shown in Fig. 3.2.1. The 
energy storage inductor stores the energy when the transistor Q is 
on, then releases the energy to the load and recharges the output filter 
capacitor when transistor is off. The key operating waveforms of this 
circuit are shown In Fig, 3.2.2. In this figure: 

I^ ■ input average DC current, 

■ output average load current, 

2d ■■ peak to peak ripple current in L^. 

The waveforms are employed to facilitate drivations of the objective 
function and the constraints. 

3.2.1 Unknown Design Variables. 

There are 22 design variables for this boost power converter 
including RLC component values, and details of magnetic design such as 
core cross-sectional area, mean magnetic path length, number of tums 
and winding area. The transitor switching frequency and converter 
overall operating efficiency are two other Important variables. 


- 99 - 




^ 1 * ^ 2 * ^5 * 
R 


I 


t 


"3 

h. tj 

i. 


C3, C^» Cg : 
^1 * ^0 * * 


Zt I Zai Ze : 


Nj. Nj : 
Acj.Acj.Acj: 

F : 

cff : 


DC winding r*«lstane«n of Inductor* Lj^, L 2 » r*ip«ctlv*ly. 
Input fllt«r damping raalstor. 

Input fllt«r Inductor*. 

Energy atorag* Inductor. 

Filter capacitor*. 

Core croaa-aecclon area of Inductor* L^i L 2 and L . 1 
reapectively. 

Mean magnetic path length of Inductora L^i L 2 and L^i 
reapectively. 

Number of turns on Inductors L^, L 2 and reapectively. 
Winding areas per turn for L^, L 2 and L^i respectively. 
Transistor switching frequency. 

Overall operating efficiency. 


3.2.2 Design constant*. 

Design constants are obtained either through manufacturer's 
specification* or designers' own experiences. Numerical value* In MKS units 
are given in the parenthesis. 


P 

D, 


ST 


BE 


'SR 


WindlnS Pitch fetor - 

Core window fill factor: (0.4). 

conductor resistivity: (0.172 x 10“^). 

Core densitjr: (7800). 

Conductor density: (8900). 

Maximum operating flux density: (0.4), 

Weight per farad: ( 0 ^ 3 , Dj^^, Dj^^/210, 1100, 72). 
Transistor saturation voltage drop: (0.2SV). 
Transistor emltter-to-base voltage drop: (0.8V). 
Transistor turn-on rise time: (O.lSvs). 
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‘sr 


'ND 


TD 


*RE 

•Si 


I TrMtlttor turn-off fall tlMj(0.2 no), 
i Dlodo conduction voltago drop: (0.9V). 

Dlodo turn-on rloo tiiMt(0.03 uoX 
t Dlodo turn-off foil tint: (0.05 mo). 
t Dlodo turn-off rocovory time: (0.03 mo). 
: Hoot oln.. volght donslty: (15.4 w/kg). 

: Sourco wolght donolty: (30.8 w/kg). 


3.2.3 Powor Convorter Porformonco Roquiremento. 


Tho porformonco roquiremento opocifiod bolov will bo omployod In tho 
next ooctlon to fomuloto doolgn constrolnto. 

Input voltogc (28V). 

Eq : Output voltago: (37. SV). 

Pq : Output powor: (70W). 

S : Froquency dopendent source conducted Interference: (O.IA) . 

This specification limits the maximum percentage of the 
switching current being reflected back to the source en- 
suring that the source is not significantly disturbed by 
the switching action downstream. Referenced here is 
mil - sta 461, whose characteristic curve is shown as 
follows : 
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t output ripplo factor: (IX). 

Tho output ripplo Cocto** dofinod oo 
output ripplo focton (?X) ■>. 
pook-to-DOok output, ripple voltoao , 

noainol dc output voltogo 

PBl : Znput-filtor rotonont poaking liMitt (2), 

Tho input filter pooking at ito rooonont froquoncy 
ohould bo linicod in order not to degrade the etebility 
and the audioaueceptibility of the converter. 


3.2.4 Objective Function. 

The objeccive function ia defined oe the total weight of the con- 
verter, which ia tho aum of varioue component weighta including: 

(a) Core Weight: W1 - Dj(A^Z^ + A 2 Z 2 + A 5 Z 5 ) 

(where AZ ■ core volume) 

(b) Winding Weight: WTW - ^C2**2*^ 

(where 4F^/^ ■ mean length per turn of the winding) 


(c) Capacitor Weight: WC 


PO 


°K3^3 °K4^4 


(d) Source Weight: WS ■ effK 

Po “ 

(where • Input power) 

Pod - efO 


(e) Heat Sink Weight: WH 

Po 


etf Kl 


(where ■ Pq " total looe) 

Oblective Function : F - WI + WTW + WC + WS + WH 


(3.2.1) 


3.2.5 Design Conetrainte, 

In this study, the design effort is carried out using appropriate 
models which portray the physical characteristics of the boost power 
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converter. Som of tho Boro olgnlf leant eharactarlstlea arat powar loaa, 

corn wltidov area, eora flux danaity, and nagnatlc vlndlut raalacaneat 

Incluaion of thaat ,;haractoria£lca raaulta in a vary eonplleatad 

oat of nonllnaar conotralnta. Mathanatlcal raatriet Iona for thaao contrainta 

nan ba found In Appandlx D. 


(a) Loaa Conatralnt t C(l) ■ 0 

C(l) - P()(^ - 1) - ?IP - PQ - PD - POP - PCAP 

whara PIP ■ input filtar copper loaa 


( 3 . 2 . 2 ) 


off E. 


(R| Rj) . 


PQ ■ Tranaiftor aaturatlon loaa 4- Base 

drive lose 

4- tranaiator turn on switching loss 
■f transistor turn off switching loss 




.ff EjEo 


.ff EjE(, 




eff 


2L5E0F 


+ 



eff E^ 


^i <gQ-V 

2Lj^P 


4 




4- 


9 


(3.2.3) 


4 






C(6) - N 2 A 2 - 


eff£^ 2 


i<Eo- 

2SV 


(c) Window area conitraint; C(7 


All the Inductor windings must be accommodated within the physical 
confinement of the available core window area. All cores employ a 
toroidal configuration with square cross section area. 


w 


Z, A-0.5 

w 


IVcsl*’'* h 

w 


(3.2.10) 


(3.2.11) 


(3.2.12) 


(d) Parasitic resistance for L, , L,, Lr*.C(2) ■ C(3) - C(12) ■ 0 


C(2) - R^Aj,^ - ApF^/aJ^ Nj^ 


(3.2.13) 


C(3) - R 2 Ajj 2 - ApFpi^ N 2 


(3.2.14) 


C(12) - RjAjjj - ApFjj/T:^ N 5 


(3.2.15) 



(•) Input fllfr p»«kln& con«traint;C(4) ■ 0 

Thip conatralut It laportant in dattralnlng tht audio- 
autcaptlblllty performanct and the control loop ttablllty. (16] 


C(6) 


(PEI)^ 


1 + 


-U-^C 

3 


3 





4^^ 

44. 


(3.2.16) 


(f) Output ripple conatralntt C(ll) ■ 0 

Output ripple factor (in percentage) Is expressed as 


'R 


Ej(Eo-Ej)- 


effEjEfl 




"6^ 




c(ii) - V. - 


E,(Eo-E,)- 


ett Ej E„ 


2L5Eo‘F 


(J.2.17) 


(3,2.18) 


(g) Frequency dependent source PlI constraint; C(13) > 0 
This constraint limits the maximum percentage of the 
switching current being reflected back to the source. 

The input filter must be designed to satisfy the following requirement: 
Required attenuation at switching frequency 

EMI requirement 

Fundamental component of the switching current 
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C(13) 




''A . 


“1 


where 




A • 


it^jF 


Hin 


»(Eq“E^) 


tO.3 


D - R, 


(3.2.19) 


(h) Additional inequality constraints; C(14) , C(15), C(16), 

C(17) > 0 

These constraints are needed to confine some of the variables 
In reasonable ranges In order to facilitate program 
convergence. 

C(14) - 0.97 - eff >0 (2.2.20) 

C(15) - - R 2 > 0 

C(16) - - 1.0 « 10 “^ > 0 

C(17) - C, - 1.0 X 10 ”^ > 0 


(3.2.21) 

(3.2.22) 

(3.2.23) 



3.3 INTRODUCTION TO NONLINEAR PROGRAMMING AND AUGMENTED LAGRANGIAN 
(ALAG) PENALTY FUNCTION METHOD 

Moat optimization problcma ariaing from practical power converter 
applicationa are aufficiently complicated to defy cloaed-form aolutiona. 

To numerically realize an optimum dealgn, one haa to reaort to nonlinear 
programming algorithms which can provide faat convergence to an optimum 
solution from a reasonable guess of the starting point. The non- 
linear programming problem (NLP) of extremizing (maximizing or minimizing) 
a function of n variables, while requiring other functions of the same variables 
to satisfy either equality or inequality constraint relationships, is called 
constrained NLP. The problem is to maximizing or minimizing a function of 
n variables without regard to side conditions or constraints is called an 
unconstrained NLP. Wltile there exist numerous methods of nonlinear programming, 
the effectiveness of each method depends greatly on the particular multldlmen- 
sioital problem to which the method Is applied. The availabilty of numerous 
efficient numerical methods for solving the unbonstralned optimization problem 
has motivated the design of algorithms that transform a constrained problem 
into a sequence cf unconstrained problems such that the successive solutions 
of the unconstrained problems converge to a solution of the constrained problem. 
These transformation methods Implicitly Incorporate all the constraints into 
the objective function that is to be optimized. The algorithms based on the 
transformation approach are conceptually simpler and easier to implement 
than the algorithms that handle the constraints directly because 
of the relative ease of extremizing an unconstrained problem com- 
pared to a constrained one. 
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3.3.1 Nonllnaar Programing Penalty Function Mathod. 

The Panalty Function Technique employa tha aforemntioned traneforaation nathod 
[ll], Let us define problem PI as the original constrained NLP problem and P2 
as the transformed unconstrained NLP problem. 


PI: Minimize objective function f(x) subject to: (3.3.1) 

1. ) Inequaltly constraints g^(x) > 0. 1 - 1,2...., P. 

2. ) Equality constraints hj(x) ■■ 0, J " 1,2,..., q. 

P2: Minimize A(x, w", g, h) , m - 1, 2. . . . (3.3.2) 

Where x - vector of n unknown variables. 


A(x,w"') ■ new objective function formed by augmenting the 

original objective function f(x) with weighted terms 
(penalty terms) that depend on the constraints g and h. 
w” • controlling weighting factor, an penalty term, a vector 

of Lagrange Multipliers, 
m ■ number of iterations. 

The essence of transforming the constrained NPL into an unconstrained NLP 
is that by gradually removing the effect of the conctralnts in the new obJec> 
tive function (by controlling the weighting factor w"*) it is possible to 
generate a sequence of unconstrained problems that have solutions converging 
to the solution of the original constrained problem. 

That is; lim A(x , w® g, h) - f(x ) ■ 0, 

nr*-« (3.3.3) 

After inlteratlons, the variable x approaches the optimum x*^ 

In effect • the influence of the constraints on the augmented objective 

function la relaxed and, in the limit, removed, and the augmented objective 

function A(x, w*", g, h) converges to the same optimum value f(x*) of the original 

objective function. 
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3.3.2 S«qu«ntlaily Unconatralncd Mlnlmliatlon Tachnlqua (SUMT), 

SUMT w«« davaiopad* valldatad, axtandad and raflnad by Placeo and 
NcCornlck {10,111. Thla amthod raplacaa cha conatralnad problaa PI daflnad balow 


by a aaquanca of unconatralnad alnlnixatlon problaaa aa daflnad in P2. 

PI: Minlmlsa f(x, y) aubjact to: (3.3.4) 

1. ) y, *) > 0, 1 - 1, 2, .... P 

2. ) hj(x, y, *) - 0, J - 1, 2, .... q 

P2: Minimize P(x,r. ) « f - r. f in g. + J (3.3.5) 

* 1-1 ^ J-1 J 


where P(x. f|^) la the penalized objective function, la a monotonlcally 
decreasing sequence tending to zero, and x Is an n-’dlmenslonal vector re- 
presenting the design variables to be optimally selected. In the power 
converter design optimization, the components of x are values of R. L. C. 
and the design details for magnetics such as core cross-section area, mean 
magnetic path length, wire size, number of turns on magnetic winding, etc. 
y represants the vector of constants related to component characterlstlca 
such -is winding and core densities, transistor and diode switching times, 
the Intended maximum operating flux density of given magnetlca, etc. 

z represents the vector of performance requirements to be met by 
optimum design such as the maximum output ripple, E3fl requirement, out- 
put power. Input filter peaking limit, etc. f(x,y) represents the 
objective function (such as the total converter weight to be minimized. 

The basic Idea of SUMT Is to solve a sequence of unconstrained 
problems like P2 whose solution approaches the solution of PI. Considerable 
computational difficulties have been experienced with the SUMT algorithm. 
The most serious handicaps are summarized below. The contours of P(x,r|^) 
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corraspond to Incraaalngly stoap sldad vallaya aa tha eontrolllnt paraaatar 
dacreaaas»and thaHaaalan of tha function baconaa prograaalvaly aora ill- 
condltlonad an 0 and tha optima aolutlon x* la approachad. Aa a 
raoult, tha aaarch dlractlona nay bacona nlalaadlng. Tha rata of convar- 
ganca dapenda on tha Initial valua of and tha mathod of raduclng 
Finally, moat of the information about the topology of f(x,y) and P(x,ry^) la 
dlacarded from ona atage to tha next even If some type of extrapolation la 
Incorporated In the algorithm. The attempts to overcome theae computational 
difficulties have resulted In several modifications of SUMT. The ALAG 
penalty function technique resulted from such efforts to Improve the compu- 
tational method. It has gained recognition as one of the most effective 
methods for solving constrained optimization problems. The algorithm based 
on this method converges at a super linear rate; the computational 
effort per Iteration falls off rapidly; the Initial starting point need 
not be feasible; and the transformation function Is defined for all values 
of the parameters. 


3.3.3 ALAG Penalty Function Technique. 

The Augmented Lagranglan Penalty Function for PI Is obtained by 

combining the Powell-Hestenes [12,13] penalty function and the Rockefellar 
penalty function [14] as In P3 In the equation below. 

P3: Minimize ip(x, X,o) ' 

*1 , 2 
where \|)(x,X,o) - f(x,y) - I [X^hj - ^jh^] 


I 


1-1 


■ 

X, 1 

CM 


1 

- ^ 

°i 

»i 




- 1 


(3.3.6) 


where 




»1- o. 


^1 

mlnKgj^ * 




Vi 


X^c E^, 0^£ a^t 


Vi 
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Th« 6^s and o^t are controlling paranatara, wharoaa tha SUMT algoritha 
haa only ona controlling paranatar Tha paraaater o^ia changed only 
whan tha rata of convarganca la not aatlafactory and Q^la changed In 
ovary Iteration to enforce conotralnt oatlafaetlon. Tha Inportant 
feature of thia approach la that l/o^» which correaponda to the controlling 
paraaeter r^ln SUMT, la not required to tend to aero for convergence 
of the algorithm. In the penalty function t|i defined above, the term ^ 
repreaenta a penallalng threahold for the 1th Inequality conotralnt. 

Increaalng o ^ to enforce footer convergence reducea the penalty threahold 
level and leada to cloaor conotralnt aatlafaction. Vfhen the Inequality 
conotralnt g^ > 0, then (or 9^) la relaxed to aero; othervlae, it la 
changed to make the correapondlng conotralnt active at the current oolutlon x. 

3.3.4 Comparison Between SUMT and ALAG 

The SUMT package was utilized In the Initial phase of this design 
optimization. Considerable effort was spent In computer coding and Im- 
plementation of the SUMT package. Although optimum solutions 
were reached, the results were less satisfactory and were sensitive to 
Initial guesses. Therefore another software package, ALAG, requiring only 

alight modification of the original SUMT Code, was adopted to explore an 
alternative means of design optimization. The trials on the ALAG algorithm 

ware quite successful. The two algorithms are compared briefly below [IS]. 

(1) Since SUMT requires first and second erdar derivatives for the 
constraints and the objective function, many hours of data preparation 
are needed. The ALAG algorithm needs only the first derivatives. 


(2) Th« Initial starting point la SUMT ohould be strictly faaslbla 
with respect to the inequality constraints in ordar to gat proper conver- 
gence. The ALAG algorithm does not require the starting point to be faaslbla. 

(3) SUMT requires about tan times more computer storage end CPU 
time than the ALAG algorithm. 

(4) ALAG converges much more readily than SlRff. 

(5) Computational comparisons between SUMT and ALAG were obtained 

f 

using Buck converter as an example [15). The comparisons are shown In Table 3. 3-1. 


Table 3.3-1 COMPARISONS BETWEEN SUMT AND ALAG 
USING BUCK CONVERTER DESIGN OPTIMIZATION PROGRAM 




SUMT 


ALAG 

COMPILING TIME 
(SECOND) 

4.61 

(USER PROGRAM 
ONLY) 

6.54 ENTIRE PROGRAM 

EXECUTION TIME 

197 



6.34 

TOTAL TIME 

201 + 

? 


13 

INPUT CARDS 

2719 



1299 

KILOBYTE SECONDS 

165216 



10652 

TOTAL RUN CHARGES 

$16.65 



$4.29 


Since the ALAG was fotmd generally superior to SUMT In several 
aspects, it was adopted to implement the computer-aided design 
optimization for the remainder of the study. 
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3.4 IHPLEMENTATION OP COMPUTER CODE AND USER'S GUIDE 


Th« nonlinear programing packago AUG la daalgnad with tha uaar 
convanlanca In mind. Thara la ona main program and nlna aubroutlnaa In 
tha whole package. A Hat of computer programa la given In tha Appendix P. 
Tha nlna aubroutlnaa are: 

AUGA 

AUGZ 

QNTA 

MUDA 

MUDB 

MUDE 

BQMA 

BQMB 

ALAGB •* User's supplied 
subroutine 

All subroutines except ALAGB are supplied in the AUG package. 

The user also needs to supply the main program Including the objective 

function, constraints and their derivatives. The flow-chart of general 
optimization sequence using a multiplier algorithm Is shown In Figure 
3.4.1. 


3.4.1 User Supplied Main Program and Parameters 

The main program basically supplies the controlling parameters, 
the Input data ^uch as design constants, performance requirement. Initial 
starting point), variable and constraint scaling factors, and the out^ 
put Information with prlnt-out format. Detailed descriptions of the 
main program and the subroutine AUGB are given In this chapter and 
Appendix E. Both of these programs may serve as a limited user's 
application guideline. 
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( Stort ) 


Initialize algorithm 
k ■ 0 

I 


Colculote search 
direction r*^ 

I 

kjj^l 

N 


Convergence or 
Stopping conditions 
satisfied 


Find ! 


Lg(x*''*'*) • min + »r*^> 


± 


Calculate search point 
status at 


Output search point 
stotus (subject to 
satisfying conditions) 


Output 

termination 

status 


( stop ) 


Conditions f(^ 
jpdatlng multipliers^ 
^gnd penalty weight] 
satisfied 


N 


update multipliers 
land penalty weights 


Figure 3,4.1 General Minimization Sequence Using a 
Multiplier Algorithm 
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Th« uMr auppUad controllint paruwtara ar« «xpl«ln«d In thn 


followlnti 

Controllint Parmaf ra 

N An Intagar aot to tha numbar of varlablaa, N > 2* 

H : An Intargar aet to tho total numbar of conatralnta* M > 1. 

K : An Intagar oat to the total number of aquallty conatralnta. 

EPS A raal array of N alamenta uaad In variable convarganca criterion. 

EPS(I) ahould ba aet auch that EPS(I)/X(I) • AFMIN, where X(I) 
la a real array of N alemanta In which tha Initial aatlmata 
of tha aolutlon la aat. 

AKMIN A raal numbar In which tha relative error toleranca required In 

tha conatralnt raalduala muat be aat. 

MAXFN : An Integer In which the maximum numbar of calla of ALA6B (uaera* 

auppllad aubroutlne) on any unconatralned mlnlmlaatlon muat ba aat. 

IPRl : An integer controlling the frequency of printing for ALAGA 

subroutine, IPRl la usually set to 1 for printing. If IPRl ■ 0, 
then no printing. 


IPR2 


IW 


MODE 


DFN 


; IPR2 ■ MAXFN + IPRl, an Integer controlling the frequency of 
printing from QNTA, the minimization routine. 

: An Integer giving the amount of storage available in CCMMON 

statement. 

; An integer controlling the mode of operation of ALAGA. A normal 
setting is that MODE - 1. 

: DFN is aet to zero. 

The important controlling parameters the user needs to change for 


different programs are N, M, K, EPS, AKMIN, MAXFN. 


3,4.2 Ut«r Supplied Subroutine ALAGB. 

ALAGB ie the only ueer eupplied eubroutine. Thie eubroutine g^vee 
Che infomecion ebouc conetreint equatione and their firet order deri- 
vativea. Tlta following atepa are uaed to prepare the eonputer code 
for thia aubroutine: 

Step 1 

Define array variablee x(l) to x(N), were N ia to the nunber of the deaign 

unknownat 

Step 2 

Manipulate the conatraint cquatlona auch that they are aiaplified aa much 
aa poaalble. 

Step 3 

Examine the conatrainta obtained in Step 2 , and aaaign a name to 
{.^^on conatant terma in order to further almplify the conatraint 
equations. Step 2 and 3 will, save computer dato preparation time. 

Step 4 

Now transform all the constraint equations in terms of the array variables 
X(l) to X(N) and the constant name created in Step 3. The computer pro- 
gram layout is such thst the Itiequality constraints come after the equality 
constraints. The constraints are designated from C(l) to C(M)i where M 
is the total number of equality and inequality constraints. 

Step 5 

Take fl;'at derivatives of the objective function and the constraints with 
respect to their corresponding variables. For example, assume constraint 
C(i) contains one variable X(j), then the derivative will be designated as 
GC(J,i). For the objective function, the derivative will be designated 
as G(j), where j is one of the variables contained in the objective function. 
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Sf p 6 

Th« eoapuCPr Cvd* for tho oubioutliio ALACB will h«vo tho following Inyouti 
OOMION 8TATIHEIIT 

EQUATE XI ■ X(l) (Thla onabloa tho uaoro to uaa XI Inataad 
of X(D) 

OBJECTIVE FUNCTION F 
CONSTRAINT EQUATION C(l) 

DERIVATIVES OF OBJECTIVE FUNCTION G(J) 

DERIVATIVES OF ALL THE CONSTRAINTS GC(J»1) 

RETURN 

END 


The user la referred to Appendix B for more progrpnmlng detalla. 


Example 

Conalder tho EMI conatrei't: 


C(13) - 






Eo .(E„-Ej) 

where A ■ — = Sin 


tr^LjF 


E. 


'c'0-5 


D - R, 


Step 1 ; Aaslgn variable name: 


^io“ 

*15- 

h, 

PEl" 


'20 


L 

C 

C 

R 


1 

3 

4 
3 
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SfP 2 ; Manipulate the aquation: 


C(13) 


_fo_ siB 4, Vc I -1 

/^^l2000j ^ PEIRj ^ 


Step 3 ; Aaalgn constant nasia: 


let XM26 • 


A J-L]i 

'/*■ 2000 

^ < 


!o_ ’<v*i> 
.2, ““ -i7~ 


XM27 - 4n^F^ 


](M28 


2sF 

PEI 


Step 4 : Constraint equation: 

Now C(13) becomes 

^7 

- 1 + (XM26)(XM27)X^X^qX2q (XM28 “ D > 0. 


Step 5 t Take first derivatives 
GC(7,13) - ^.0(XM26)(XH27)(XM28) 


-(XM26)(XM27)Xj^qX2q 


^ 7 ^ 10^20 

*15 


GC(10,13) - (XM26HXM27) X^X2q (XM28 - 1) 


GC(15,13) - -(XM26)(XM27)(XM28) 


*7*10*20 


"15 


*7 

GC(20,13) - (XM26)(XM27)X,X.n(XM28 T“ - 1) 

/ lu Xj5 


STEP 6: ■ Set up computer code for subroutine >ALAGB. 
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3. A. 3 Initial Starting Point and Scaling T«chnlqu« 
3,A»3vl S<l«ction of Initial nfrting point 


For a complicated nonlinear optimisation problem with 20 or 
more varlnblea, proper aelectlon of the Initial starting point plays 
an Important role In the speed of convergence and accuracy of the 
aolutlon. As a general guideline « the Initial starting point should 
be selected such that the value of each equality constraint la as 
small as possible. This point must also satisfy the Inequality 
constraints In order to stay In the feasible region. A properly 
selected Initial starting point will speed up the rate of convergence. 

The time and effort spent to choose a good Initial point prior to 
running the program is well worth the result. For a practical 
problem, choice of a good initial starting point can usually be 
based on the designer's past experience, or on some simplified 
design guidelines and equations. 

3.4 .3 . 2 Variable scaling technique and convergence 

In a switching power converter design, the values of the design 

variables are scattered over a wide range. The capacitance may be 

in the order of 10 , for example, and the switching frequency In the 

order of 10^. This wide scattering of values Is one of the primary 

causes of convergence difficulty. Therefore, a variable scaling technique 

is provided in the computer program to scale all the variables between 

■•6 3 

values of 1 and 10. For example: if Xj^ “ 0.5x10 and X 2 ■ 0.8x10 , 

-7 2 

then one can use VSCAL (1) ■ 10 , VSCAL (2) - 10 , so tnat: 

*1 ^2 
VSCAL (1) " vscSTTi) “ 

where VSCAL (1), VSCAL (2) are scale factors for the respective variables. 
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For reasonably acceptable accuracy, the tolerance for variable 
convergence (EPS) Is set around EPS 

to exist from the Iterative computation through the variable 


-6 “7 

10 to 10 . For the program 


convergence criterion. It must satisfy the following requirement: 


(k) 


(k-1) 


< EPS 


This requirement states that the largest difference between two values 
of any variable from consecutive Iterations must be less than the tolerance 
required. It should also be mentioned that the program can exist via a 
constraint convergence criterion. 


3.4. 3. 3 Constraint scaling technique and converRence 

It Is very unlikely that the Initial starting point can satisfy 
all the constraints to the extent that each equality constraint 
residual Is smaller than the constraint tolerance and each Inequality 
constraint Is also satisfied. If the starting polng did satlslfy all 
the constraints, then, of course, the problem would already be solved. 

In reality, the constraint values based as Che Initial guess can 
vary over a wide range. Since conditions where certain constraint values 
may be so large that Che effects of other constraints are obscured 
should be avoided. It Is desirable to scale each constraint by such 
a factor that Che effect of violating any given constraint Is of the 
same order of magnitude as the effect of violating any other 
constraint. Unfortunately there are no universal guidelines for 
selecting the constraint scaling paremerters. It has been 
observed that faster convergence can be achieved by the proper 
selection of these parameters; however. Improper use of 
constraint scaling can cause divergence problems. Experience shows 


4 
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2 -2 

that by scaling tha constraint values In a range between 10 and 10 • 

■•3 

end setting the constraint tolerance around. 10 , the progran can achieve 

a faster rate of convergence. 

For 0n acceptflble dccuracyp the constraint tolerance AKMIN is set 

around lO”^ to lO”^ using the scaling technique. Vfhenever the maximum 

scaled constraint violation AKK' ‘ is less than AKMIN, program 

convergence is reached. This stopping criteria can be put in a more 

concise way in the following. 

; Constraint value for ith constraint in iteration k 

SC^ : Scale factor for the 1th constraint 

(k) 

WW^'^ ^ : Scaled constraint violation for 1th constraint In Iteration k 


AKK 


That is 


WW 


(k) 


(k) 


SC. 


(k) 


: Largest scaled constraint violation in iteration k, that is, 

AKK^*^^ - max 

^ 1 


Whenever AKK 
is terminated. 


(k) 


< AKMIN, the convergence is reached and computation 


The program sometimes can also be run without using the constraint 
scaling technique. Experience shows, however, that by using the con- 
straint scaling technique the program can be brought under better control. 


3.4.4 Stopping Criteria for Computation 

For normal exit, there are constraint convergence criterion and variable con- 
vergence criterion as mentioned In the previous subsections. In roost 
cases, constraint convergence is deemed more desirable. The accuracy of 
the result however depends on how to choose the constraint tolerance and 
variable tolerance. In some circumstances, the solution via variable 


convergence criterion la sufficiently accurate to be acceptable. Infinite 
looping and obnormal exit are also possible as shovm In the flowchart 
In Fig. 3.4.2. 

Discussion of flowchart 
EXIT 1 

This exit means the objective function has been evaluated a number 
of times equal to the user's supplied paramteter MAXFN. The solutions 
from this exit are In moat cases not accurate. The user may Increase 
MAXFN to get proper convergence and more accurate solutions. 

EXIT 2 

This exit means the largest scaled constraint violation Is less than 
the constraint tolerance. The solution from this exit Is deemed most 
desirable. 

EXIT 3 

The exit means the largest difference of variables between consecutive 
iteration is less than the variable tolerance. Depending on the exit 
condition, the solution from this exit is often acceptable. 

LOOP 4 

In this ted tape loop, the program Is never converged. The user 
must set execution time limit or printing page limit In case endless 
loop occurs. 

3.4.5 Checklist for Computer Printout 

The following checklist is provided for the user to assure the final 
solution Is accurate and acceptabl' . 

(1) Check If the solution C Is In the feasible region, that Is, If the 
Inequality constraint residuals are all greater than zero. 
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it (2) Check the accuracy of Che aolutlon obtained. That ia: 

check if largeat acaled conetralnt violation ia leaa Chan AKMXN. 

(3) Check if MAXFN 18 REACHED, if the exit ia normal or abnormal. 

(4) Check if the aolution X ia reaeonable uaing common aenae and 
previoue experience. 

If the aolution obtiined is not accurate enough then one can follow 
the flowchart as shown in Fig. 3.4.3, by changing the starting point (using 
^ the final result of the previous run) or readjusting the scaling factors 

and rerun Che program. 


4 ■ 


€ 


% 
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CHECK MAIN 
AND ALA6B 
STARTING POINT 
PARAHCTER 

VARIABLE 

SCALING 

. „ I : 

CONSTRAINT 

SCALING 

'I 

SET TIME AND 
PAGE LIMIT 
RUN PROGRAM 

INTERPRET 

RESULT 


RESULT 

ACCEPTABLE 


TERTURB 
STARTING 
PDINT ^ 


STOP 


INCREASE 

MAXFN 

MAXFN 

INCREASEI 


MAXFN 

.REACHED. 


Figure 3,^1. 3 Flowchart of Effective Programming Approach 



3.5 OPTIMIZATION RESULTS OF BOOST CONVERTER 


The translator switching frequency Is a critical paraoMter In the 
minimum weight design of switching power converter. In the course of 
optimization the frequency is held constant for a computer run and a 
sufficient number of runs are collected carrying the frequency over a certain 
range of Interest. Several distinct advantages can be obtained In this 
approach : 

(1) By treating the frequency as a constant In each computer run, the 
nonlinear optimization problem Is simpler and reaches convergence easily. 

(2) Important design Insights can be obtained when the weight and 
loss are plotted against frequency. Instead of identifying a single 
optimum switching frequency and a minimum weight design as in our earlier 
optimization attempt [8], the curve presented here provides a range 

of frequencies in which the system weight Is minimized In all practical 
sense. The curve will also provide Information regarding sensitivity of 
the converter weight as a function of the switching frequency. 

(3) The trade-offs between weight and loss as a function of switching 
frequency can be evaluated readily. This Information can be used as a design 
guideline for the weight /efficiency optimization. 

By treating the switching frequency as a constant In the optimization 
process, a set of design data as a function of switching frequency Is obtained 
by varying the frequency between 20KHz to 120KHz In a lOKHz step. The design 
parameters specified In Section 3.2 are employed to make these computer runs. 
Detailed design results Including the detailed loss and weight breakdowns oi 
various components are collected and tabulated In Table 3.5-1. The mlnlaum 
weight converter design data Including the details of magnetic design 
for each chosen frequency are shown In each column. This helps in the 
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Table 3.5-1 Boost Converter Optimization Results 


•••••■■•nt of Ch« opciwia cooponont doolgn •• a function of awltehlng 
froquoncy. Tho daalgnatlons uCllicad In tho loan and waight braakdowna 


are shown 

below. 

PQ 

m 

Total power diaaipation in the transiator 

PD 

m 

Power dissipation in the diode 

PCAP 

m 

Power dissipation in the output filter capacitor 

PMAG 

■ 

Total magnetic losa (core loss -f winding loss) 

PI 

m 

Total loss 

US 

m 

Source weight 

UH 

m 

Packaging Weight 

UI 

m 

Magnetic core weight 

UC 

SB 

Capacitor weight 

UU 

m 

Winding (copper conductor) weight 

WMA6 

m 

WW + WI 

WT 

m 

Total weight 


In order to gain more design insights, the total-weight/ total-loss 

and the component weights/loss breakdowns are plotted against the switching 

frequency in Fig. 3.5.1 and Pig. 3.5.2, respectively. Sunsuirlzed In the 
following are several Important observations from the table and curves. 

(1) The curve of total weight versus frequency exhibits U-shapa 
characteristics. Tha. converter weight Is heavier at both low 
frequency end and high frequency end. 

(2) The total weight of the converter reaches its minimum value in 
the frequency range from 30KHz to 50KHz. 

(3) The U-shape curve is also observed by plotting the total loss 
characteristics against frequency. The rapid decrease of the 
total loss at lower switching frequencies is caused by Che 
reduction of winding losses of the magnetics, meanwhile the Increasing 
total loss at higher frequencies Is caused by Che higher switching 
















lo«««s and much rapidly Incraaalng matnaclc loaaaa. Aa a raault of the 
Incraaaa of total loaa at high fraquanci«ij(, the aourca weight and packaging 
weight alao lacreaae rapidly. The weight reduction due to the decreaelng 
magnetic component weight aa frequency Increaaeala leaa pronounced than 
the Increaae of packaging and aourca weight. Therefore a U-shape curve 
of total weight vs. frequency la formed. The U-ahaped curve Implies that 
there exists an optimal switching frequency. 
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3.6 DESIGN OPTIMIZATION OF BUCK-BOOST CONVERTER 




Th« Buck-Boost switching povsr converter Is chosen ss a second sxaaple. 
The same procedure used for the design of Boost converter Is now applied 
to Buck-Boost. The circuit schematic and problem formulation are stated 
in the following. The Input-output relationship and derivations of the 
constraints are given In Appendix F. The results of design optimization 
are demonstrated In Section 3.7. 

3.6.1 Circuit Schematic, Design Variables, and Waveforms 

The circuit parameters and design unknowns are shown In Fig. 3.6.1. 
There are 24 unknown variables including the details of magnetic design. 

This circuit contains a two-stage input-filter, a two-winding energy 
storage Inductor, a power transistor, a diode and an output filter. When 
the transistor is turned on, the energy from the source Is stored In the 
two-winding energy storage Inductor ; the output filter capacitor 
supplies power to the load. When the transistor Is switched off, the energy 
previously stc i^d in the inductor is dumped out to the load where it also 
replenishes the output-filter-capacitor energy. 

The operating waveforms are shown in Fig. 3.6.2. These waveforms are 
used in derivations of design constraints such as the output ripple factor, 
and EMI constraint, etc. The notations marked on the waveform are defined 
as follows: 

2d ■ peak-peak ripple current through transistor Q 
■ average input current from source E^ 

Iq ■ output DC current 

n - prlmary-to-secondary turns ratio of energy storage inductor, 

(- 1 In the design example presented In Section 3.7). 
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transfomer 

R3 ; Input filter resistor 

ly I2 : Input filter Inductor 

Lp ; Primary Inductance of two winding 

Inductor 

Cy C5 : Filter capacitor 

A2> Ap : Cross section area of Inductor or 

transformer 

Zi^ Z2/ Zp j Mean magnetic path length 

tiy ^2* Np : Number of turns of the winding 

Ac2> Aq 2> A(-p; Winding area per turn 

F : Switching frequency 

eff : :^yerall efficiency 

Figure 3 . 6.1 Circuit Schematic and Design Variables of 
Buck'Boost Converter 


- 135 - 



3.6.2 Objective Function end Constraints 


The notatlone used previously for the Roost converter ere adopted here 
for Buck-Boost. 

3.6.2. 1 Objective function: total 

Core Weight W1 - + 

where AZ ■ Core volume 
Winding weight WTW ■ 

where 4F^v^ ■ mean length per turn of the winding 
Capacitor weight ■ 

**0 ^0 

Source weight « WS ■ eff k" * w**®*^® power 

Pn (1-eff) 

Heat sink weight WH > 

**0 

where “ ^0 " power loss 

Objective function • W1 + WTW + WC + WS + WH (3.6.1) 



weight 

*2h ♦ *p^> 

t Ac2N2*^ 


Input filter copper loss 


0 


eff E, 


(Rj + R2> 


- Transistor saturation loss -f Base drive loss 
+ transistor turn on loss + transistor turn off loss 


VsT /-‘VbE 


off E off Ej 


T F 
SR 
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PCAP - Output filter capacitor ESR loss 
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C(14) - 0.97 - eff >0 
C(15) - >0 

C(16) - - 1.0 X 10"^ > 0 

C(17) - - 1.0 X 10’^ > 0 


( 3 . 6 . 19 ) 
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3.7: OPTIMIZATION RESULTS OF BUCK-BOOST CONVERTER 


By Creating cha awiechlng fraquancy aa a conacant In aach opcialca- 
clon run, a aat of convarcar paraaatar data la obcalnad. 

This set of parameter data represents the optimum converter design 
for Che specified switching frequency. A sufficient number of runs 
are executed by varlng the Irequency between 20 KHz to 120 KHz In a 
10 kHz step. Detailed optimization results, following the afore- 
described design process are collected and tabulated In Table 3.7-1. 

To facilitate comparison of optimal converter designs between the 
boost converter and the buck/boost converter, the same Input-output 
requirements, design constants, and converter performance specifi- 
cations are used. (Reference to Section 3.2 for detailed Information.) 

The turn-ratio n ■ N /N * 1. The designations employed in Table 3.7.1 

1 O 

are the same as those In Table 3.5-1. To provide more design Insights, 
the weight and loss breakdowns are plotted against the switching 
frequency in Fig. 3.7.1 and Fig. 3.7.2, respectively. 

The difference between converter optimization results for the boost 
converter and buck/boost converters are summarized as follows: 

(1) The Buck-Boost converter Is heavier than the Boost converter. 

In order to have the mlnlnum weight design the Buck-Boost conver- 
ter has to operate at a higher frequency than Boost power converter. 

(2) Switching losses of semiconductor devices are higher for the 
Buck-Boost converter. This Is logical since the switching current 
amplitude Is considerable higher than that of the boost converter 
for the same Input and output voltage and .the same power level. 

(3) The magnetic component for the buck/boost converter are generally 
larger In size and heavier in weight. 
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(4) Th« aagnatic losa«* (PMAG» cor« lost •¥ winding loss) for ths 
buck/boost convsrtsr srs doainstsd by ths winding loss in low 
frsqusnciss. Ths PMAG loss chsrsctsristic falls rapidly «s ths 
switching frequency increases. The high magnetic losses in low 
frequencies cause severe weight penalty. It is clearly demonstrated 
in Pig. 3.7.1 and 3.7.2 that in order to minimise the converter weight/ 
loss, it is desirable to operate the converter frequency about 
80KHs ~ lOOKHz. For the minimum weight/loss boost converter design, 
however, the optimal frequency rest about 40KHs * 60KHz. 
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Results 


145 



Figure 3.7.1 Weight Breakdowns for Buck-Boost Converter 
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3 . 8 CONCLUSIONS AND SUGGESTED FUTURE WORKS 

3.8.1 Conclusions 

Nonlinesr programming techniques have been successfully employed to im- 
plement the minimum-weight design of switching powe* converters. Two different 
computational algorlthmv - ALAG and SUMT-based on the penalty function method 
were compared and their figure of mertis were assessed. For power converter 
optimlzation» the ALAG package was deemed more effective than its counterpart 
the SUMT package, when the computation time, ease of coding, and rate of 
convergence are concerned. 

Adopting the ALAG rountlne, a cost-effective computer-aided design approach 
is presented which provides a minimum-weight converter design down to the details 
of component level and concurrently meets all power-circuit performance require- 
ments. This computer-aided design approach provides important design insights 
which helps to assess the following important design concerns: 

(1) The trade-offs between weight and loss as the switching frequency is 
Increased . 

(2) The optimum converter design down to the details of component 
levels. 

(3) The optimum component designs as a function of the switching frequency 
and their relationships to the overall system optimization. 

(4) The significance of the U-shape curves representing total -ws ight/ 
total-loss _rsus frequency as observed in the collected sub- 
optimization runs. This allows the designer to easily identify the 
optimum switching frequency or a range of frequencies over which 
the total weight/loss Is minimum in tlie partlcal sense. 

(5) Impact of various critical component characteristics, such as mag- 
netic losses, switching losses of semiconductor devices, to the over- 
all system. 
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(6) The optimal converter topology for a given application. 

Employing the nonlinear -program baaed optimization technique, the power 
converter designer can conceive the overall optimum system design taking 
into consideration the powt r-clrcult relate^ performance requirements with 
Che design objective of either minimizing weight, loss, or any other physical 
realizable quantity. It thus sets the stage for a more scientific design 
approach Instead of subjective brute-force, trlal-and-error, piecemeal design. 

3.8.2 Suggested Future Works 

The investigations of complex converter upClmlzaCion problems using 
nonlinear programming techniques have shown marked success. Demonstration 
of the buck converter optimization in the previous modeling and analysis 
phases sponsored by NASA, the half-bridge converter optimization sponsored 
by NAVY, together with the boost and buck/boost converter optimization 
presented in this report have collectively provided clear evidence that 
a large scale converter optimization is feasible using NLP techniques; 
yet, the d.-velopment of such a tool has not reached the stage of maturity 
where it can be widely used. Presently, it takes a person with considerable 
insight to the nonlinear programming algorithms, and with sufficient 
converter design experience, to make the program converge. It is our 
belief, however, that the afore-described NLP techniques could be made 
easier and more systematic than they are now. The following tasks are 
suggested as means fur Improving the NLP techniques to make them a more 
universal converter design tool with wide user applicability. 
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(1) Systematic way of Improving initial starting point. 

(2) Meana of optimising variable scalingiand constraint scalings. 

(3) Maans of optimising convergence stopping criteria, 

(A) Improved method of formulating nonlinear constraints to enhance 
convergence . 

(5) Establish conditions for convergence. 

(6) Program cransportability . 


4.0 INVESTIGATION OF CURRENT- INJECTED MULTILOOP CONTROLLED SWITCHING REGULATORS. 

4.1 INTRODUCTION 

In recent years, vast amounts of Interest and research In universities 
and Industries have been directed toward development of a multi-loop, multi- 
state control scheme which could be applied to switching regulators. This 
collaborated effort has resulted In astonishing Improvements of stability 
and dynamic performance of switching regulators. 

A host of control schemes has emerged many of which employ the principle 
of current-injected control [17,18,19,20,21]. These control schemes share, 
the common property of transforming a switching converter from a voltage source 
Into a current source. This control concept has exhibited many desirable 
properties such as Inherent over-load protection, stable and equal load sharing 
when several power converter modules are In parallel, and t. i system response. 

Illustrated in Fig. 4.1.1 Is a buck/boost converter employing current- 
injected control. The control Is Implemented by sensing the output voltage Vq 
of the converter and the Instantaneous current Ip through the power switch. 

The duty cycle signal is terminated when switching current ascends and Inter- 
sects the threshold voltage Vj^ (dc error signal) determined by subtracting 
Vg from the reference voltage E^. Since the switch current waveform sensed 
by the current transformer contains both the dc bias current component, and 
the small amplitude ac modulation signal (to be used for additional error 
compensation), the control thus provides Inherent transistor peak-current 
protection (from the dc-current component) and Improved dynamic performances 
(from the ac modulation signal). 
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Presented In this report is the modeling and analysis of the current- 
injected control system. The modeling approach employed in the present 
paper, a departure from previous efforts ( 17,20), provides additional 
insight to the current-injected control characteristics which failed to 
be manifested in the previous modeling and analysis efforts. To facili- 
tate comparison between the method presented in this report and the ap- 
proadi employed in the earlier attempts, a brief review of the earlier work 
is provided. The concept of major loop and minor loop was employed in the 
previous works (17,20). By considering the dc-feedback and compensation 
network being the major loop, and the ac-(swltch current) feedback as being 
the minor loop, the minor loop was lumped into the power stage in the 
process of modeling. The multi-loop converter was thus reduced to a single 
loop system as shown in Fig. 4.1.2. The ac feedback loop which contains 
the switching current information is embedded in the "new" power stage. 

A A 

The transfer function Ox' the new power stage has a surprisingly simple 

form (only a single pole and a single zero) . The authors feel that while 
this modeling approach offers a vray to examine certain small signal charac- 
teristics of the system, little information is provided regarding the rela- 
tive stability of the system (the concept of the gain margin and the phase 
margin). Even though the dc loop can be opened, the ac loop is inherently 
closed in the "new" power stage model. A true open loop characteristic, 
where both the dc and the ac loop are opened, is thus not accessible in 
this modeling approach. 

The modeling and analysis approach of the multi-loop current injected 
control presented in this report eliminates the aforementioned modeling 
dilemma. Following un approach similar to that described in the 
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author's prcvioua work [22,23,24|2S]» tha aaall signal aadal of tha con- 
vert ar In Flg.4«l*l la darivad as shown In Fig. 4.1.3. Tha saall signal 
■odal has tha following faaturaa: 

<1) Tha power atage has thraa Inputs and t%R> outputs. 

Tha three Inputs ara: 
line disturbance 

A 

load disturbance 

A 

duty cycle disturbance d 
Tha two outputs are: 

the output voltage Vg 

the switch current ip 

A 

(2) The error processor senses the two modulation signals 1^ and 

A 

v^. The transfer function represents the gain of the ac loop 
and represents the combined gain of the dc loop and the 
compensation network. 

(3) The duty cycle modulator Is represented by a describing 
function Fj^. 

Employing the above described small signal model one can readily examine 
the following performance characteristics: 

(1) The control-to-output characteristics v^/v 

2 2 

where A ■ S + 2CWgS + Wg. (This Is the characteristic examined 
In the previous papers [l7,20] which exhibits a single-pole and 
and single-zero). 

(2) The open dc loop characteristic 

WDcfDl'«'*VAc'D2>- 

(3) The open loop characteristics (open both dc loop and ac loop) 

S “ A ^AC^D2^* 

The open loop characteristic Gj Is used to examine the relative 
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•eablllty of th« lystm. 

(4) Th« audlosMcopelblUty choractorlatle • v^/v^ 

®A " 1 + Cj { ^ VaC^D2^ “ ^DiVaC^U21 }. 

(5) Th« output lapodanco charactorlatlc - ^o^^o 

^o " 1 + Cj { ^ ^U12 ^^'*’ VaC^D2^ ” ^U22^AcVd1 }. 

Modeling of the power-stage, error processor, and pulse modulation 
Is presented In chapters 4.2, 4.3, and 4.4, respectively. Various open 
and closed-loop performance characteristics are evaluated In chapters 
4.5 and 4.6. Effects of dc-loop and ac-loop gain, and of compensation 
networks, are also discussed. Finally, guidelines for selecting control 
circuit parameters are provided. 


j 
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4.2.0 CUIRENT- INJECTED BUCK/B008T POUEE STAGE MODEL 


4.2.1 Power Circuit Deecriptlon 

The function of the dc-dc converter ie to proceee end trenefer 
electric power from an unreguleted input to a regulated output Vq. 

The output voltage of the two*winding buck-booet converter ehown in 
Fig* 4. 2. 1.1 can be either greater than or leas than the input voltage, 
depending on the duty cycle of the awitch and the turn ratio of the 
atorage Inductor. The magnetically-coupled wlndinga, provided by 
the energy-atorage inductor, allow Input/output laelatlon and multiple 

outputs. Proper choice of the inductor's turn ratio alee can alleviate 

the difflcul*”' of implementing the extreme duty-cycle eonditlom due to vide 

rangea of input and output voltages. 

The energy exchange transpires in the power stage (Fig. 4.2.1.1(a)) 

in the following fashion. The state of the duty-cycle drive d(t), 

shown in Fig'. 4.2.1.1.(b), determines the instantaneous position of 

the switch. The high-level of d(t) indicates the conduction, or on-state, 

and the low-level determines the off-time, or off-state, of the power 

switch. During T^^>a voltage approximately equal to v^(t) (neglecting 

the losses due to the parasitic resistance in the primary winding) is 

established across N_ which causes a current i (t) to increase as 

P p 

illustrated in flguro 4.2.1.1.(c). Occurring slmultaneeusly, a voltage 
Vg(t) is induced across by transformer action, but no conduction 
is allowed because of the reverse biased diode. As the off-time, T^^^ 
is initiated, the energy associated with is transferred to by an 
ampere-turn redistribution. As a result, the current reduced 



to Z9V0 and igCt) la f oread Co a magnitude needed to maintain a con- 
clnuoua MMF flow through the inductor. The output voltage In Figure 
4.2.1.1.(e) la kept nearly conatant due to the large capacitive filter 
at the output which absorba the pulaatlng currant lg(c) and delivera a 
dc current with minimal ripple to the load. 

4.2.2 Analytical Impleirv* cation 

The awltched dc-to-dc converter, aaaumed to be nondlaalpaClve, 
is nonlinear In nature. The baalc dc-to-dc voltage conversion la 
achieved by repetitive switching between o number of linear networks 
switches and diodes. The number of linear networks In one switching 
cycle is determined by the mode of operation of the Inductor's magne- 
tomotive force, MMF. If the MMF Is continuous as shown lu Fig. 4. 2. 2.1. (a), 
the power stage model has two linear networks corresponding to a mode 1 
type operation. For a discontinuous Inductor MMF operation as shown in 
Figure 4.2.2.1.(b) the power stage model is composed of three linear 
networks corresponding to a mode 2 type operation. 

Each linear circuit model Is described by a set of linear state- 
space equations. For the current-injected buck-boost power stage the 
state variables (Independent variables) are customarily the magnetic 
flux, 0, and the capacitor voltage, v^. The total number of storage 
elements determines the order of the system. 

4.2.2. 1 State Space Averaging Technique 

To derive a linear model for the power stage, the averaging tech- 
nique Is used. (27) Employing the mode 1 operation as an example, the 
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power stage is modelled by two intervals of operation, and ^Fl, 
respectively. 

(1) interval T ^ (ii) Interval T_, 

on Kx 

X ■ A^x + Bj£ X ■ AjX BjU (4. 2. 2.1) 

X “ + EjH Z. ■ ^2— 

The principal of the average method is to rf'rlsce the state-space 
description of the two linear circuits by a single state-space descrip- 
tion which represents the approximate behavior of the system through 
one cycle of operation. Taking the average of both intervals and 
summing the results yields the following line.;r time-varying continuous 
system: 

X - d(A^x + Bj^u) + d' (A^x + B^) (4. 2. 2. 2) 

jr - d(C^x + Ej^) + d' (C^x + E 2 U) 



where Is the period of the switching cycle. 

The basic requirement for the average method is that the effective 
filter comer frequency of the switching converter be much lower than the 
switching frequency 127]. 

The linear time-varying equations (4. 2. 2. 2) can be rewritten in the 
following form: 

X ■ Ax + Bji (4. 2.2.3) 

■ C X + E^ 

where A ■ dAj + d'A^ (4. 2. 2. 4) 

B - dBj^ + d'Bj 

C - dCj^ + d'C^ 

E - dEj^ + d’E2 

- 162 - 
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4.2. 2. 2. PTturbation 

To study th« •■•11-slgnal behavior , tha linear tine-varying aquatlona 
(4.2. 2. 3) are perturbed. The introduction of input variations and duty-cycle 

variations in turn perturb the output and state vectors. Tha perturbed 

a 

input vectors are: 

A A 

u ■ lU + ^ and d ■ D d (4e2e2s5) 

where U and D are the steady-state valuea and u and d are small perturbations. 
These perturbations In turn lead to following: 

A ^ 

X - X + X and Z “ 1 Z 

where + X steady-state values and x and y are small perturbations. 

With the corresponding perturbations substituted Into equation (4. 2. 2. 3) the 
basic model becomes: 

■ X - + BU + + Bu + [(A^ - A 2 )X + (B^ - B 2 )Ujd (4.2.2. 7) 

(dc term) (line (duty ratio variation) 

variation) 

+ [(Aj^ - A2 )x + (Bj^ - B2)u]d 
(nonlinear second order) 

Y + i - CX + EU + Cx + Eu + [(Cj^ - C 2 )X + (E^ - E 2 )U]d (4.2. 2. 8) 

(dc term) (line (dute ratio variation) 

* variation) 

+ l(C^ - C2)x + (E^ - £2)0] d 
(nonlinear second order) 

The perturbed state-space description is nonlinear owing to the presence 

A A A ' 

of the product of time-dependent quantities x and u with d. 

4. 2. 2. 3. Linearization 

Since the ac variations are very small in magnitude compared to their 

f 

steady-state value, the following small-signal approximations can be made: 
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(4.2. 2. 9) 




u 


1 




<< 


A ^ 



1 


AAA AAA 

Where u - (u^,U 2 ), V - (U^.U^), x - (Xj^.Xj), e» 

Using the approxlmatlone (4.2.2.9)» the nenllneer k^teond erder tens In 
equations (4. 2. 2. 7^ and (4. 2. 2. 8) can ba neglected, resulting In a linear 
system. Separating the steady-state (dc) and dynamic (ac) parts of the 
linearized system, the final state-space model Is acquired. 

Steady-state (dc) model: 

X - - a'^BU (4.2.2.10) 

Y - CX + EU (4.2.2.11) 

- (E - CA”^B)U 


Linear Dynamic (ac) model: 

s 

X » Ax + Bu + [(A^ - A^X + (B^ - B 2 )U]d (4.2.2.12) 

^ - Cx + Eu + l(C^ - Cpx + (Ej - E 2 )U]d (4.2.2.13) 


4 .2 . 2 . 4 . Transfer Function Representations 

In small-signal analysis, particular Input /output relations 
(transfer functions) are needed to construct the basic building blocks 
necessary to fully describe the power stage model. To find the input- 

A 

to-state variable x and Input-to-oatput transfer functions, one assumes 
the ac duty-ratio variation is zero. The dynamic trtodel described in 
equations (4.2.2. 12)- and- (4.2.2.13) can be simplified as follows: 

e 

A A * 

X ■ Ax + bu 

AAA 

■ Cx + D£ 

Taking the Laplace transformation of the previous equations, the 
following relations are obtained: 
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i 


■x(t) - X(o) - iU(t) -t- Bu(s) 

1 

And 

i<«) - C[«I - A]~^Bu(a)-t-Du(«) 

- ICI*I - A1"^B + Dju(«) (2.2.1A) 

« . __j * 

Lcting x(o) ■ o, the flrot equation becomes jc(b) ■ ~ A] B^ls). 

To find the duty cycle-to-state variable x and duty cycle-to-output 
transfer functions, one aesumes the ac variation of u Is zero. Equations 
(4. 2. 2.1?) and (4.2.2.13) yield the following: 

e 

X - A^ + l(Aj - A2>X + (B^ - B2)Ujd 2.2.15) 

i - Cx + l(C^ - C 2 )X + (E^ - E 2 )Ujd (4.2.2.16) 

A A 

And the resulting duty ratio modulation d to state-variable x and duty 

A .. 

ratio modulatlcn d to output ^ transfer functions are: 

1 

-X - (SI - n [(A, - A,)X 4 (B, - B-)U] (4.2.2.17) 

d(s) ^ ^ ^ 2 - 

Z(®> _1 

- - - C(SI - A) ^ ((A- - A,)X + (B, - BJUJ + (A.2.2.18) 

d(8) ^ ^ 12- 

I (C^ - C2)X + (E^ - E2)U] 

• ' These transfer functions will be used as building blocks to construct 

the power stage transfer functions to be presented in the following 
sections. 

4.2.3 Power Stage Analytical Model 

The objective of the power stage model is to develop a group of 
transfer functions that describe the low-frequency behavior of the 
switching circuit. The model developed is comprised cf three Inputs 
^ and two outputs. From Fig. 4.2.3. 1 the three inputs are the supply 
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i 




voltage the output current 1^» and the duty ratio d. The two outputa 

A A 

aupplled are the output voltage v^, and the awitching current 1^. 

The two-winding buck/boost power atage ia ahown in Figure 4.2.3.2(a). 
The two linear equivalent circuit models for the continuous MHF case are 
illustrated in Fig. 4.2.3.2(b) mid (c). The power stage model contains 
an ideal switch and a diode. The storage Inductor is a linear core cir- 
cumscribed by a primary and a secondary winding with Inductances Lp and Lg 

2 

respectively, where Lp ■ (Np/N^) L^. Also described by the model is the 
winding resistances Rp and R^/ The output filter is represented by a capa- 
citance and an equivalent series resistance, ESR. On the output a current 

A 

source 1^ Is employed to represent a disturbance Injected to the converter 
from the load. 

The power stage model Is composed of two Independent variables. The 

state variables for both linear circuit equivalents are the magnetic flux 

^ of the core shared by the primary and secondary windings N_ and N , and 

* s 

the capacitor voltage v^. During the Interval T^^, the power stage Is des- 
cribed by the following: 


4 

X 


where , 


X 




“ 

t ^ ■ 


• I " 

’^o“ 


f 

^0 


1 

p 


(4.2. 3.1) 
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For the Interval the power stage takes the following form 

X - A 2 X + B 2 U Z • C 2 X + E 2 H (4.2. 3. 2) 
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Appedlx A gives the derivetlon of the two equlvelent llneer taodeis used 
to describe the power stage during and 

4.2. 3.1 Average Model . 

Continuing the process of characterising the small-signal model 
needed to describe the dynamic, ac behavior of the power stage, the 
average model la: 


X ■ Ax -f Bu 


jr ■ Cx + Eu 


(4. 2. 3. 3) 


A - 



- D* 


"s + 


w 


-D'(Rj,//Rl) 

-D D* 

(Rj,+ Rj^)C ~ (Rj,+ Rj^)C 



<V'V 

s 

D(Rj.//R^) D(Rj,//Rj^)l 

-R^ * V 


c - 





D(Rj,//R^^) D'(Rg//Rj^) 
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(D D’) Rg//Rj^ 


where D ie the eteedy-stete duty ratio end D' ■ 1 > D. 
To simplify our matrix expreeslons, assume: 

(1) Rj, « Rj^ so that 


R- + 




( 11 ) ^ ^ 

Also by definition the following conditions are used; 


<‘> >-e ^7 


(ill) u) 


2 A 1 

o “ L C 
e 


(iv) D = 


A T 

A on 


(11) R 


A Rr. 

« d -2 


r A % 
C --2- 




and D ■ 


. A "^Fl 


where T Is the switching period, 
P 


so that D D* 


T + T T 

:ss.lhi -!i. 1 

Tp Tp 


Applying these simplifications the average model becomes: 




- D* 


D 

«■ 

A - 

L (D')2 
e 


B • 

«P 




- 1 


A 

1/C 

•M 


D' 



U 


(4.2. 3. A) 
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N„R., 



• S C 



L D' 


c 

e 




E - 


ON 



-L^ “ 


0 0 

P 




4.2. 3. 2. Linearized Powtt Stage Model 

A ^ 

Line voltage variation v^, duty cycle varlatlf^is d, and output 
curnnt source disturbances are now Introduced Into the circuit such 
that 


V “V *f V 
I I 1 


1 6 0 + i 

0 o 


d ■ D + d 


. • 4 _ I * 
d - D - d 


(4.2. 3. 5) 


where is the dc Input voltage and D and D* Is the steady-state on-tlme 
and off-tlme duty ratios# Employing the small-signal approximation, as 
discussed In Section 4. 2. 2. 2, the second-order nonlinear term may be 
neglected, and a linear system obtained. 

4. 2. 3.2. (a) Steady-state (dc) model. 

Using method of section 4. 2.2. 3, the steady state output vector Y 
may be seen to be: 

- (E - CA“^B)U 


(4. 2.3.6) 
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4, 2, 3, 2(b) Th» Input-to-output tfn»fT function . 

Th« Input-to-output transfer function duo to sn input variation 

A A 

u aaauning tha duty ratio d (ac) variation is saro is as follows t 


z<«) 


Vo(a) 


S<*> 


* + Eju(a) 


^ F “o (vH 
♦ *?) 


«c ( 


S *f 2,S + W 2^ 

1 O / 


-DD'Nj 

W 


Rj,Ca -f 1 

k 4 


(4.2. 3.7) 


-f 2;o) s + <D ^ 

o o 


where z. ■ 




\ D’ 


4.2.3.2(c) Duty cycle-to-output transfer function 

A 

The duty cycle-to-output transfer function assuming u ■ 0 Is: 


d(8) 


_v^) 

d(s) 

y») 

L d(») 


-1 


- C(SI.,') ((A^ - A2>X + (Bj - B2)U1 + 

I (Cl " C2>X + (E^ - Ej)])] 


-17.3- . 






y(n) 


d(«) 


DD (S^ 2Cw S 4 b> 

' 0 0' 


1- 


s 


u C 

‘ O 


♦ » ♦-r 

• D 


VqNs 

Vp 


(r^cs + D + l) 
(s^ + 2C« S + uj) 

' o o ^ 


1 

~r 

D 


(4.2. 3. 8) 


For detailed derivations of the transfer functions, refer to 
Appendix I. 

2 2 

Further aimplifcatione are made by letting A ■ S 4- 2Cb>oS -f and 
then forming gain expressions. THe gain blocks in matrix form due to a 

A A A m 

variation u “ [v^ i^]^ (4. 2. 3. 7) results In Che following: 


\ , A 1 

Z(*) ■ T 


ull 


u21 


ul2 


u22 


(4.2. 3. 9) 


-L 

A 


N 


s D 


Np D' “o 


(Rj,CS + l) Rj, + zjS + Wo^Z2 J 


-DN_u 

b O 


(RI.CS-H) («jCS+l) 


(4.2.3.10) 


From the duty cycle-to-output transfer function (4. 2. 3. 8), the matrix 
simplification yields the following: 
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(4.2.3.11) 


z(«) 


4 


'D1 

^2 


d(s) 


2 % 




4 


Vi 


I •* 




S X , - *^C 
o 


(Rj^CS + D •♦• 1) ^ 

V ^ -F 


d(«] 


(4.2.3.12) 


The block diagram for cha amall-aignal approximation of tha power 
stage la ahown In Fig. 4.2. 3. 3. Tha analytical axpraaalona for tha 
and Fjj_ funcClona are glvon In equations (4.2.3.10) and (4.2.3.12). 










4.3.0 error processor MODEL 

4.3.1 Introduction 

The error processor, EP, Is a feedback compensation network, It 
processes multiple-input control signals derived from the power stage, 

and delivers the required analog Information to the pulse modulator. 

Fig. 4.3. 1.1. Illustrates the analog-signal error processor employed 
V In the present analysis. Prom a small-signal viewpoint, the EP is 

a linear network consisting of two control loops. The signals sensed 
In these two loops are the converter output-voltage v^, and the primary 
switching-current 1^. 

The loop sensing v^ Is the same as any conventional dc loop, where 
the sensed v^ Is processed by an amplifier with reference voltage E^^, 
to generate a dc error signal v^. The dc error voltage after processing 
through an Integral plus lead-lag compensation network establishes as 
the threshold level for the switching-current Information derived from 
the ac loop. 

The integral plus lead-lag compensation network Is employed to 
^ shape the frequency response to Improve the converter stability and 

dynamic response. 

The ac loop which senses the collector current of the power 

t 

transistor serves two functions. The first function is to transform 
the primary switching-current Into a proportional voltage signal. 

This voltage signal Vg^ Is then compared with v^ resulting in a control 
that turns off the power switch when v^. The second function of the ac 
loop is to derive the low-frequency modulation signal or error signal for 
additional loop compensation. 
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.'3.2 Error Proc«*«or Ac «nd Dc Loop Galna 

Fig. 4.3. 2.1 Is a functional equivalent of Che elreult 
shown In Fig. 4.3. 1.1. In Fig. 4. 3. 2.1. the two error signals, namely 
and are subtracted to form the positive input to the threshold 

detector. The negative input of the error processor is now replaced 
by a zero reference. 

The authors feel that the modified error procesaor better serves 

to perceive the modeling effect for the following reasons: 

(1) As stated earlier, the ac loop contains two types of 

information: The large amplitude switching-current 

waveform is used to implement the analog-to-dlgltal 

conversion^. Such a function is considered part of the 
* 

pulse modulator instead of Che error processor. The 
small-amplitude low-frequency modulation signal (similar to 
the modulation signal sensed by the dc loop) is brought together 
with the error signal from the dc loop and compensation loop to 
provide the total state-feedback compensation for improved 
stability and dynamic responses. For the purpose of modeling 
of the EP, the large-amplitude switching-current information 
should be extracted from the ac loop and Incorporated into the 
pulse modulation model (presented in Chapter 4.4). 

(2) Modeling the EP is difficult because of the unconventional circuit 
implementation of the threshold detector. In the conventional 
design, the threshold detector is implemented with one fixed 
threshold voltage Input, the other input containing the error 
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Information. With the current-injected mode of control, both 
Inputs to the threshold detector contain error signals or ■ 
delation signals. In order to model the whole converter system, 
the pulse modulator should take the combined error signal from 
the EP output and convert It Into a pulse-width modulator duty- 
cycle signal. In order to Implement such a "single-input single- 
output" pulse modulator model, various error signals derived from 
the multiple feedback paths should be combined to form a compo- 
site error signal containing error Information from the dc, ac 
and compensation loops. 

(3) For reasons mentioned above, the original EP circuit is modi- 
fied to that of Fig. 4. 3.2.1. where the positive Input, v^, 
to the threshold detector contains error signals from all three 
feedback loops. The error voltage v^, is the output of the 
multi-loop EP and also serves as the input to the pulse modulator. 
The functional equivalence cf Fig. 4. 3. 1.1. and Fig. 4. 3. 2.1. can 
be justified in the following way. Mathematical)' < >t.’ any v^^ and Vg^ 
waveform, the block diagram of the circuit in Fig. ^(.3.2.2. (a) is 
equivalent to the analytical model In Fig. 4.3.2.2(b). The equivalence 
Is proven through the following arguement. For the pulse modulator control, 
the following equations hold true: 


If Vjj >. Vg„, then v^„ - 1 


(4. 3. 2.1) 
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To fora a composita arror algnal aa tha Input to tha thraahold 
datactor, tha dc arror voltage and voltage are combined.^ The 
following relations hold for the modified circuit of Fig'. 4. 3.1.2: 


''X " "'SW ^ ° ^TH " ^ 

“ “ ^SW ^ ° '^TH • 0 


(4. 3. 2. 2) 


Equations (4. 3. 2.1) and (4. 3. 2. 2) are exactly the same; therefore, 
Fig. 4. 3.1.1 and 4. 3. 1.2 are functionally equivalent. 

In performing the small-signal analysis of the network shown in 
Fig, 4. 3. 1.2, the reference is replaced by a short circuit. Apply- 
ing Klrchhoff's current law at node A yields the following: 


i_ + 1_ ^ A— 

Rj^ R 2 R 3 


v V. 

Rl R 3 


(4.3. 2.3) 


where v^ is the voltage across R 2 , is the differential input voltage 

A 

of the operational amplifier, and v^ is the output voltage. 


As for ncde B the Klrchhoff equation is 


' 1 

+ 4T- + BCt 

■ V 

— + flC, 

*x 

Rj 3J 

X 

i^x ’J 


R, 


(4. 3.2.4) 


'X ■ *4 * ic7 • (4.3.2. 5) 

The oen-loop response of the operational amplifier is A(s), where 
V » - A(s)vjj ^ (4. 3. 2. 6) 
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A(«) - -K 

1+8 /Wn 


(A. 3. 2. 7) 


Substituting squstlon (4. 3. 2. 6) Into (4. 3. 2. 3), and (4. 3. 2. 4) 


results In the following set of equstlons: 


" 1 . _L . 

''e R, Rj 


R^ A(s)R3 


(4.3.'2.fl) 


-V 







1 

A(s) 


+ 

— + »c. 



R, 


(4.3. 2. 9) 


Eliminate the node voltages and by simultaneously sclvlng 

A A 

(4. 3.2.8) and (4. 3. 2.9) for and v^. A single expression In 
terms of and results: 


1 


A(s)R. 


1 + ^ + ^ 
'^l ^2 


1 

A(s) 




— +sC- 
ZX 3 


R R3 

— + + 1 

Rl R2 


(4.3.2.10) 


Since the ac and dc expressions ace low-frequency models, the operation- 
al amplifier gain Is A(s) = -K. With a very large K, equation (4.3.2.10) 
can be simplified to the following form: 

^ ' - (^.3.2.11) 
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Th« ac loop ia daacrlbad by tha following: 

A 

whara n la tha turn ratio of tha currant tranoforawir. 

* } _ j 

'^aw " aw *^SW " p n 

From Flgura 3.1.2 the following equation la derived: 


(A. 3. 


(4.3. 


Substituting equations (4.3.2.11) and (4.3.2.13) Into equation (4.! 
yields the following: 

1 




Zx 


‘'SW : 

~ n ^p 


(4.3.2 


The dc loop gain for the EP Is defined by the following: 

A 


DC V 


S-o 


1 

’ 1 + sR^C 

1 



a 

GP 

•^1 

(€^+€ 3)8 

l + sR, 
4 

[V 3 ' 




, " 1 ^. . 

d 




(4.3.2 


where GP - 


Ra Rm 
1 + — + 

Rl R2 


The ac loop gain Is given by the following: 


r- T 

■ r 

P 

. ^sw 

n 


Vq • 0 


(4.3.2 


Fig. 4.3.2. 3 Illustrates the error processor block-diagram. 


.12) 

.13) 

.14 

1.2.14) 

.15) 


.16) 


.17) 
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4.4. 


PULSE MODULATOR CHARACTERIZATION 


4.4.1 Pulse Modulator Descriptions 

The pulse modulator converts the analog error signal from the 
output of the error processor into a modulated pulse-train in order to 
provide proper duty-ratio control of the power switch. Ths pulse modu- 
lator analyzed in the present chapter utilizes a constant off-time 
control. 

Fig. 4. 4. 1.1 illustrates the pulse modulator (PM) in simplified 
block diagram form. The two inputs of the threshold comparator are the 
dc error voltage v^, and the ac voltage signal v^^ (proportional to the 
current through the power switch). The dc error voltage v^ is a floating 
threshold level proportional to the difference between the desired and 
the measured output voltage. As the error between the desired and actual 
voltage Increases, the level of the voltage v^^ rises. Since the ac voltage 
Vg^ Is proportional to the current waveform through the power switch, peak 
current protection Is an Inherent feature of the current-injected control. 

UUen the condition Vg^ > v^^ is satisfied, the threshold detector 
output V.J, is momentarily driven low. This signal instructs the digi- 
tal signal processor (DSP ) to turn off the power switch. Therefore by 
limiting the level of v^^, the peak current protection of the power 
switch can be achieved. 

For constant off-time implementation the DSP output, D, is latched 
low for a predeteradned amount of time. At the end of the off-time 
period, the power switch is commanded on and the cycle repeats. 










4.4.2 SMll>8ignal Kodtl 

To study ehs sasll-signsl bshavior of tho pulss Modulator, an 
analytical nodal is davalopad using tha describing function tachnlqua. 

Pig. A. 4. 2.1 shows tha block dlsgran of <:.ho analytical nodal at 
staady-stata without a low fraquancy dlaturbatica. Aa ahown In 
Fig. 4.4.2. l.(b), whan Is graatar than connancaa. S|| la 

tha rising slops of tha switching-currant signal during tha 
psrlod Sp Is tha falling slops datarnlnad by tha sacondary 

currant of tho two winding Inductor roflactad back to tho prlnary aids 
during the period T^p^. The trajectory formed by and 8p Is pro- 
portions} to the inductor MMF and is used to >}atarnlne continuous or 
discontinuous current operation. 

Fig.’ 4.4. 2. 2 shows the equivalent circuit waveforms for the per- 
turbed case. The small-signal disturbance is represented by an ac 
generator placed in the dc loop as Illustrated in Fig. 4.4.2.2(a). 

For Che perturbed case, the magnitude of the low-frequency ac 
signal is assumed to be sufficiently small and the rising slope 
S^, and falling slope Sp, are considered constant and unaffected 
by the disturbance. With this assumption, the perturbed switching- 
current Information Vg^, along with the error voltage can be 
modeled. 

The development of Fig. 4.4.2.2(b) is « follows: 

(1) Given a slope determined by the pov. ge's supply voltage 

Vj, tho prinary inductance Lp, and an initial dc ctartiag value 

large enough to guarantee continuous-current operation. From 

tha starting point a line is drawn with k s slope 
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At the intersect Ion of the perturbed end the evitchlng- 
current Information begins. After s set spsn of tine 

equal to the period expires, the cycle repeats. 

(ii) The initial point for the next cycle is detemined by subtracting 

the magnitude '^OFF particular point in time. 

(iii) From the new initial values (i) and (ii) are repeated. 

Fig. 4.A.2.2(< > shows the input to the analytical threshold 
detector. Using the relation ** development of 

Fig. 4.4.2.2(d) can be understood. Fig. 4.4.2.2(d) illustrates the 
perturbed duty-cycle. 

The small-signal behavior of the PM can now be modeled. A set of 

equations are developed for the darkened trajectory v.^ of Fig. 4.4.2.2(c). 

v^ is used because the trajectory describes the low-frequency modulation 

affecting the PM gain. The i.uditional information in v.^, that deals 

with the dc component of the switching waveform of the converter is not 

relevant in the modulation of the dutv cycle signal, and therefore has no 
contribution to the PM small-signa^ gain. 

4.4.3 Formulation of the PM Transfer Function 

The output of the PM can be expressed in a Fourier aeries in the 

form: 

d(t) ■ D + sin <i)t + bj^ cos ut + ••• 

The input of the PH can aloo be expressed in a Fourier series in 
the form: 

v^(t) ■ V + sin wt + dj^ cos wt + ••• 

The describing function F|^ of the pulse modulator is defined as: 


(4, 4. 3.1) 




A 3<*^) _ 


(•i 

(Ci*+ <*1 )** 




tan 


-I h 


-I **l' 
can 


f # 

The reaultant waveform la a function of and Sp (both aaaumad 
poaltlva) and the ainuaoidal dlaturbance K aln wt. During the period 
Tqpp at(N-l)^^*' cycle has the form: 

v^(t) - Sp(t-t 2 ^_j^) + A aln ut - A aln (4. 4. 3. 2) 

During the period at the cycle has the form: 

v,^(t) - ^ i»)t - A aln (4.4. 3. 3) 

The on-time ATj^ at the n^h cycle U2^»t2u^2^ expreaaed through 

the following development: AT^ can be derived from equatlona (4.4. 3. 2) and 

(4.4. 3. 3) by evaluating the equations at t - t^^, where 

''T<‘2n> ■ *F“^2n-'2n-l> * * "‘2n ' * “4n-l (4.«.3-‘) 


''i‘'2n> ■ S»<‘2n+1-‘2D> * * **" “‘2n ‘ ^ “‘2n+r 

From Figure 4.2.2, Tp 4 Tgpp - t2n"'2n-l’ 


Let 



^2n+r^2n-l • 


(4. 4. 3. 5) 


Set (4.3.3) equal to (4.3.4) and solve for AT^. 

SpTp + A .in - A .in + 

A aln - A .In 

uATj. AT 

(SF + Sj,)Tp-Sj, AT„ - -2A sin -j- cos m(t 2 „_i + — ) 


(4.4. 3. 6) 


(4.4. 3. 7) 
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r 


AT, 


(Sp + S^)Tp 


N • - A(i) CO* <i>t2n-l 






A(it 


N ‘ 1 - ~ co« wt 

‘N 


( 4 . 4 . 3 . 8 ) 


2n-l 


The duty-cycle output d, nay be expressed In the following fom: 

d(t) A D + Sj^ sin (dt bj^ cos wt + ••• (4. 4. 3. 9) 

The coefficient of sin wt In equation (4.4. 3. 9) can be expressed as: 


(1) 

‘l • V 


r p^l /^2n+l 

I sin ot dt -f I sin wt dt ••• + sin «t dt + ••• 


«- o 


"2n 


] 


W 1 

ir A 


- tj I 

I A sin o)t dt + I 


A sin (dt dt + • • • 


o 


1 1 


N 


n A ^ ^ 2 n+l 

n-n 

AT^-Tp 

where A2„^., ^ “ 4 n ^ “^ 2 n+l^ 


(4.4.3.10) 


AT„-Tp 


- (ATjj-Tp)A sin “(^20-1 2 — ~ ^ 

- A(ATjj- Tp)|sln <*>t 2 n-l^** cosidATj^ + % cos(dTp] 

ajkiwApj f (4.4.3.11) 

A first-order approximation for is: 


+ cos sl»dAT|j+*l slnidTp]| 


A 2 n+i = A(ATj, - Tp)sin <*>t 2 „_i + j (AT,J - Tp^ )cos «t 2 „_^ (4.4.3.12) 
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Eliminating second order effects: 


^2n+l ” ‘^2n-l 


From equation (4.3.7) » 


*1 " w A ^ " AT„ “^2n-l 

n»o N 


(4.4.3.13) 


(4.4.3.14) 


where K is the ratio between the switching frequency and the modulation 
frequency. In order to simplify the analysis, it is assumed that K is an 
integer. 

If AT|j is small compared to the disturbance period, i.e. if K is very 
large, it follows that 


'i ■ s 1 




ST * 


t'U l-Slcoa'iit 




• sin ut dt 


^ Vu, 

- — ■ — ^ I cos (i»t sin u)t dt ■ 0 
" . ®F I 


(4.4.3.15) 


Equation 4.4.3.15 is orthogonal. Therefore the coefficient a^^ is zero. 

The coefficient of cos wt in equation (4.3.6) is: 

rr*"! . (^3 f^2a+l “I 

bj^ ■ ^ I I cos <Dt dt + I cos wt dt + • • • + j cos wt dt + • • • 

I- o 


cos wt dt 


■ ^ (A sin wtj^ + A sin wt 3 - A sin ut 2 + •** + A sin 


- A sin wt» + • • • 1 


where, 

B 


2n+l ■ * “*2n+l ’ * "hn 


■ “ “* 2 *' 211+1 ■'■‘2n**‘" 2 *'2n+l"'2n* 

- 2A coe to (^2n-l'*’^^N”T’ ^ ^ 2 — 

A firs', order approxinatlon for equation (4.3.15) la: 

A” - T 

“20+1 ■ “‘ 2 n-l 


(4.4. 


^ (4T, - tp)co. 


From equation (4.4.3.13) 
N 


1 A^t 


n«o N 


(4.4, 


With small compared to the modulation period equation (4.4.3. 
be expressed as follows: 

2 ”/ to r 




Ao) 


0) 1 


TT S 


F 


1 - , - /t— (1 cos ti»t) 
Sp Sj, 


cos ut dt 


2 7to 


Ab) 


2S 


(1 'f'cos 2(i>t)dt 


1 +^ o 
‘’N 


N 


.\(0 




(4.4 
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3 . 17 ) 

3. IB) 

,3.19) 
19) can 


.3.20) 




Again* whan Bodallng tha analog voltage v^, tha Infonution that naada 
to ba axtractad la contalnad In contalna tha lotr^fraquancy 

Information neceaaary without tha hlgh-fraquancy avltchlng algnal of 
. From aquation (4.3.1) tha amall-algnal Input voltaga to tha pulaa 
modulator la: 

v^Ujn) - SyTy-A .<n + A .in 


U) 


0) 


■ Vr + “ “• 1 <'2n*‘2n-l>*‘" 2 


(4.4.3.21) 


wTp uTp 

In equation(4.4.3.21) sin — can be approximated by -y reaultlng in: 


'i ’ ^F^F ^F “<'2n-l * T > 


(4.4.3.22) 


From equatlon'(4.4.3.22) the fundamental of the small-algnal average, .of 


v^ can be taken. 

T 

Ao) Tp cos u(t ~Y ^ 


(4.4.3.23) 


The describing function for the pulse-width modulator i.s: 


4 iifX . >“ 


M 


(4.4.3.24) 


t t 

Substituting equations (4.4.3.15), (4.4.3.20) and (4.4.3.23) into (4.4.3.24) 
the PM gain Is: 


A(i) 


IFmI 




>I AiJTp ' 


(4.4.3.25) 


(iiTp 

Phase 1. Fj^ - — , 


■ (S^ + S,)Tp 




(4.4.3.26) 
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The slope* and Sp can b* defined by the system parameters. 

During T 

• on, 

dl„ 


IpRp 


+ L, 


P dt 


(4.4.3.27) 


A s' 

Let -rr- ■ K where S' is the slope of during T . Then: 
dt " p ^ on 


Vi - IpRp + LpS' 


and 


'N 


Vi - IpRp 


(4.4.3.28) 




The magnitude of IpRp is much less than Vi, so (4.4.3.28) is approximated 
by: 


'N 


(4.4.3.29) 


The current transformer used to sense the switching-current reduces the 
sensing current by a i:n turns ratio resulting in the following 
expressions: 


'N 


nL„ 


(4.4.3.30) 


The signal sensed by the comparator is a voltage proportional to the 
current passing through R^^. Therefore 


“ -nr s, 


R • S " 
N SW 


During Tqppi Sp msy be found through the following: 




'^^^OFF'" ®F ^OFF 


(4.4.3.31) 

(4.4.3.32) 

(4.4.3.33) 

(4.4.3.34) 
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lh« tubacltutlon of oquatioM (4.4.3.33) and (4.4.43.34) into aquation 
(4.3.32) rasulta in the expreeelonat 


^ON " ^OPF 

Tqn 

^ ^OPP 


Applying aquation (4.4.3.31) to Sp • we obtain 

R, 


(4.4.3.35) 

(4.4.3.36) 


Np Vq 


s' ■ V R 

®P Tqpp ''i nLp Ng nLp ‘'sw 


(4.4.3.37) 

Substituting equatlona (4.4.3.31) and (4.4.43.37) Into the PM deacrlblng 


function the following can be obtained: 

.^£p 

2 ^ 2 


nL„ 


'ON 

*‘0FF 


+ 1 


(4.4.3.38) 


‘OFF 


J— 


"Tp 

e 


(4.4.3.39) 


where << 1. Tp the off-tlme period of the awltch and la the 

period of the low-frequency modulation. Therefore the F^^ gain nuiy be 
approximated by: 


2nLp 1 

* V,R„.. T«„ + T, 


(4.4.3.40) 


I*'SW ‘ON OFF 

Expressing the F({ describing function In terms of the duty-ratio Instead 


of Tqjj .-nd Tqpp* equation (4.4.3.40) becomes: 




2nLp 

1 

’'I'sw . 

N V 1 

1 +_P _o 

*OFF 

^ ^ N V 
S I J 


(4.4.3.41) 
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Fig. 4. A. 2. 3 la tht block dlagran for tho •Mll-aigna' gain froa tha 
coapoaita arror procaaaor algnal v,^ to tha pulaa aodulator output d. 


d(t) 



(■) 


Fig. 4,4. 2. 3. Pulaa aodulator gain block 
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4.9.0 PERFOIMANCI ANALYSES AND TEST VERIFICATIONS 

4.5.1 Block Dlasroa 

Tho objocclvo of Chi* **€11011 1* Co lncorpor*C* fh* power *C*e*» 

*n*log error proceeeor* end pule* nodulecor enell-elgnel aodele Inco e 

cloeed-loop block dlegren. The block dlegrea aodel will be ueed Co 

exemin* Che open-loop end cloeed-loop perfoniince* of ch* ewlcchlng 

reguleCor. Fig. 4.5. 1.1 ehow* ch* block dlegrea for Che buck/booeC 

regulecor derived from Fig. 4. 2.3.3, Fig. 4. 3.2.3, end Fig. 4. 4. 2. 3. 

The enell-elgnel enelyel* will Include conCrol-Co-ouCpuC 

reeponee, dc open- loop behevlor, oyeCem open- loop reeponee eteblllcy, 

V (*) 

cloeed-loop reeponee audlosuecepClblliCy, and ouCpuC Impedance 

In Chi* chapCer, a Bode ploC showing Che gain and phase of 
each particular transfer function Is presented. The analytical 
curves are compared with the experimental measurements to verify the 
small-signal model. The analytical curves are presented by solid lines 
and the experimental results by Che following: 

(I) Gain curve - x 

(II) Phase curve - • 

The gain end phase plocs for che smell-signal analysis employ Che 
following converCer parameters end operating conditions: 
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Convrf r Paramef r« 


OpTatina Condition 


L 


S 


R 


S 


R 


«L 


464.94 mlcrohanriat 
101.64 mlcrohenrlaa 
.293 ohms 

549.70 microhenries 

I. 58 ohms 

441 microfarads 
.07 ohms 
15 ohms 
86 turns 
40 turns 

II. 47 kilo-ohms 
.10.0 kilo-ohms 
23.3 kilo-ohms 
105.0 kilo-ohms 


11.8 nanofarads 


■ 22.6 picofarads 

CTN ■* n - 200 turns 
- 250. ohms 

u 

o 


Vj - 18.5 

- 11.4 
o 

D - .57 
d' - .43 

^OFF - 11 


vjlts 

volts 


0 microsscs 


2021 rad /sec 


4.5.2 Control-Co-Output Reapont* 


V 


Tha control-to-output reaponaa la tha opan-loop gain froai 

A 

to tha output v^(a) In Fig. 4. 5.1.1. The oontrol-to-output Cranafar 
function denotes Che power stage characteristics plus ac feecback. 
Figure 4.5.2. 1 illustrates the gain block diagram for the control- 
to-output response, derived from the system block diagram. 

From Fig. 4. 5.2.1 the control-to-output transfer function la 
derived, yielding the following results: 


<^^(8) - 


^ ^M*^AC^D2 


(4.'/. 2.1) 


Where A - s +2 s+u) 2 

Substituting equations (4.2.3.12), (4.3 2.17) and (4.4.3.41) Into the control- 


to-outpjt transfer function the following equations can be generated: 
2 


Vn 

0 o 


V^(8) 

^) 


“W 8+1) 



DL 

R 1 


" 


8 + 

R +t:t 

~-i 



A 

e D* 




(s^2CV+“o'>*VaC^ 


(Rj^C s + 


(s^ + 2 fl+at ^) 


(4. 5. 2. 2) 


Rc D 

Employing the Inequality (R^ D*") ^ equation (4. 5. 2. 2) is simplified 


V (s) -F V w ^ (Rj-Cs+l) 
o^®^ . ^ro o 
n — ; — 


DL« 


s- 1 


v^(8) 


DD 


F 2 

(l+p)s + 




(■4, . 


■-S*' “■ 


8+« +F 
o 


2 4 


D+1 ***o 

V O’ 


(4.5. 2. 3) 


where F ■ F„F 


'^0^8 


IfAC Rj^Np 
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Fig. 





Equation (4.5.2. 3) may be further simplified to read 


v^(.) 


■Wo 

(D* +F)D 


(RcCa + 1) 


DL, 


s-1 




2^0 + 


FD’ 


■•S 


. 2 . FD' 

•■^“o 


(IHl) 


( 4 . 5 . 2 . 4 ) 


Equations (4. 5. 2. 4) is the control-to-output transfer function with 
the angle expressed in phase delay. 

4. 5. 2.1 Test Verification for Control-to~Output Response 

Fig. 4.5. 2. 2 demonstrates the method used to measure the control- 
to-output response. An ac signal, A sin wt, is injected into the dc loop. 
The magnitude A is adjusted to provide the optimal signal-to-noise ratio 
( a clear and stable read out). Excessive amplitude of the injected 
signal could result in distortion of the modulation waveform 
and should be avoided. Channels A and B are connected as shown in Fig, 

4. 5. 2. 2. Channel B minus Channel A results in the desired control-to- 
output response., , 


Fig. 4. 3. 2. 3 shows the gain and phase of equation (4. 5. 2. 4) for the 
dc-dc converter. The solid-line projection represents the analytical 
solutions for the control-to-output response. 


4.5.3 DC Open-Loop Behavior 

The dc open-loop response is the open-loop gain derived from 
opening the dc feedback loop. Due to the nature of the current- 
injected mode of control , the true system open-loop response cannot be 
measured. The dc open-loop response can be used to indicate the 
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FREQUENCY (HZ) 


Fig 4 5.2.3. Theory end measurement of control-to~output characteristics 
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r«lftClv« ttability of th« tyfeon. Fig. 4. 5. 3.1 illuotratM eh« block 
dlagrui tot tho dc open- loop bthavior. 

Rafarring to Flg> 4. 5.3.1 tha dc open- loop cranafar function can 
bo axpraaaad In tha following fora: 


"Vdc^DI 


(4. 5. 3.1) 


-DL A + F^^,F^2 

The general expreaalon can be formulated by aubatltutlng equatlona 
(7.2.3.12), (4.3.2.16), (4.3.2.17) and (4.4.3.41) Into equation (4. 5. 3.1) 


The following formulation yields: 

+ F V(ju)^ f 1 fl 

■ 


l+aR^Cj^ 



f^l^3 ] 

1+8R^ 

c.-fcJ 


1 y 




(8-+2Ca,^8+u.^- ) +F„F^c ^ ^Rj^Cs-HHlJ + 


8 "rZCu) 0+U) 
o o 


(4.5. 3. 2) 


f Rcl D Vs 

rT 1 and F » F^F^g , equation (4. 5. 3. 2) becomes: 




•4F)D(C,+C,) <RdCB+l)(j^^ 8-1) 


l+sR^C^ 




. j. 2 . FD' (D41) 
8 + 2C‘*U'’’f+D' L ' ■ ■*' "*“• • - 


(4. 5.3.3) 











Equation (4. 5. 3. 3) la tha de opan-loop tranafar function for tha 
cur rant-lnj acted converter. 

4. 5. 3.1 Teat Verification tor DC Opan"»Loop 

Tlia dc-opan-loop mcaaurenenta are taken by Injecting an ac algnal 
Into the dc loop. After traverelng the loop, gain and phaae of the 
ac elgnal la measured. Fig. 4. 5. 3. 2 llluetratee the measuring technique 
employed for the dc-open-loop response. 

The analytical and experimental waveforms for the dc open-loop gain 
and phase are shown In Fig. 4.5. 3. 3. The analytical waveforms are 
represented by the solid line and the experimental gain measurements 
are denoted by the "k" sign. The phase measurements are represented by 
the sign. 

4.5.4 System Open-Loop Response 

Th<> system open-loop response characterises the converter's 
stability. Because of the multiple feedback paths of the converter, the 
loop is opened at a node coimson to all the feedback paths. Analytically, 
the loop can be opened in this fashion. However, experimentally this 
characteristic can not be measured. In current-injected con- 
trol the ac signal and dc information are fed Into a comparator. 
Physically, the only location common to the dc and ac loop Is the output 
of the threshold detector. Unfortunately the comparator output Is a 





OPENOC 








logic l«v«l signal. Evan for tha aodifiad aquivalsnC circuit as shown in 
Fig. 4. 3. 2. 1| the signal contains both digital and analog inforaation. By 
opening tha loop at of Pig. 4. 3.1. lor v^of Fig. 4.3. 2.1 tha waasuramant of 
small-signal gain and phasa would losa its physical meaning. For this 
rasson, only tha analytical solution for tha system open-loop response 
is given. Fig. 4. 5.4.1 shows tha block diagram for tha system open- 
loop response. 

From the Fig. 4. 5.4.1 tha system open-loop response can be derived. 

The resulting equation is given below: 

®OL “ T ^*^DC^D1 ■*'*^AC^D2^ (4. 5. 4.1) 


Substituting the converter gain expressions (4.2.3.12), (4.3.2.16), 
(4.3.2.17) and (4.4.3.41) into (4. 5. 4.1), the following formulation results: 


8^-*-2C0) 8+U»^ 

o o 


(€^+€3)8 


1 f sR, 


Cl 


\ 



GP 

f^l^3 1 


R, 

jcj+CjJ 


1 


Vo-o 





r 


+ F ® ® - 

ACRj^Np 1 
w 


o^+2C(i» s+w^^"' 
o o 
+ ?n 


(4.5. 4. 2) 


with « 1, equation (4. 5. 4. 2) is simplified yielding: 

*4. 
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Mg. 






1 


°ol’ 


dd‘ 


o o 


-C.P 

l + sR^Cj ] 

(Cj-Kj)R^S 

> 

Id-iR^ 

FlS 1 



[^. il 


.2 . 

R-D* 

. + 12:. 

l^ J 

Wo 

m ▼ 

[4 '”1 



(4.5. A. 3) 

Equation (4. 5. 4. 3) la tha ayatam op( n-loop tranafar function In 
ganeral form. The axpraaalon will be further simplified In chapter six. 

The analytical waveforms for the gain and phase response are 
shown In Fig. 4. 5. 4. 2. It Is interesting to point out that the system 
open-loop gain approaches a constant value and the phase delay 
approaches zero as frequency increases to a high value. This character- 
istic has made the current Injected control unique in comparison with 
any other type of control. Hlgh-galn, wlde-bandwldth and stable operation 
could be accomplished simultaneously. In the analytical model, the pulse 
modulator gain Is assumed constant with no phase delay. However, it Is 
pointed out In reference ( '] that the reduction of pulse modulator 
gain together with an Increase of :he phase delay was observed in 
laboratory measurement. Having Incorporated the non-constant gain and 
phase delay of pulse modulator model In the system open-loop character- 
istic, the gain characteristic would carry a ceftaln slor<s and the 
phase delay would not be approaching zero. 
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4.5.5 Audlosusceptiblllty 

Audloausceptiblllcy la the cloaed-loop, input-to-output reaponae 
for the dc-dc converter. The audloauaceptiblllty characterlatlc la 
uaed to evaluate the rejection rate of the propagation of a alnuaoldal 
disturbance from converter Input to output. Figure 4. 5. 5.1 lllustratea 
the block diagram for the closed-loop analysis of the current- Injected 
regulator. The expression obtained by examining Fig. 4. 5. 5.1 la the 
following: 


^ull r, VaC^DzI ^DiVac''u 21 

‘ L ‘ J' 

v~7^ * F F F F P F (4.5. 5.1) 

1^ ^ ^ ) !rAc'^D2 I ‘^diVdc 

A A 

In equation (4. 3. 5.1), the denominator can be expressed In terms of 
the system open- loop response After substituting equation (4. 5. 4.1) 

Into (4.5, 5.1) the following expression is given; 


Vo( 8 ) 

7^) 


^ulL^^^ ^M^AC^D2^ ~ ^Dl^M^AC^uZl 

U + <=0L> 


(4. 5. 5. 2) 


The general form for the closed-loop transfer function can be 
expressed by substituting equations (4.2.3.12) , (4.3.2.16), (4.3.2.17) 
(4.4,3.41) and (4.5, 4.3) into equation (4. 5. 5. 2). 


v^(s) [ N.D 

Vj(s) ■ ■ “o ( V®*!) 




8^+2 CU) 8+ (a) ^ 1 

o 0 I 

D’ J 


VacV 


V!s 


^5j^(«cC+l) (-iq^s-l) (R^C4l) 
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D’ 


► • 


*!S 
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i. 
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Wo 

rr- 


DD' 


(s^ + 2Cw 8+w 

O 0 


-GP 


(C^4C3)RjS 

1+sR^ 

f^iS 1 

[C1+C3JJ 




+ 


Wo 


8 ^+ 


RlD* 


+ 2CW 



(4. 5. 5. 3) 


Equation (4.5.5, 3) will be simplified In chapter six. 

4. 5. 5.1 Test Verification for Audlosusceptlblllty 

The audlosusceptlblllty measurements are taken by Injecting an ac 
signal In series with the Input voltage of the converter, accomplished 
by Injecting the perturbation through a transformer whose secondary 

winding Is In series with the supply voltage. Figure 4. 5. 5. 2 demonstates 
how the measurement technique Is Implemented. 

Fig. 4.5. 5. 3 Illustrates the gain curves for equation (4.5.5, 3 ) 
for the current-injected regulator. The analytical curves are 
portrayed by the solid lines and the experimental result 16 
represented by the x curve. 

4.5.6 Output Impedance 

The output Impedance Is employed to measure the dynamic perform- 
ance of a switching regulator subjected to sinusoidal load disturbance. 

A switching regulator with zero output Impedance represents an Ideal 
voltage source. In a linear system, the output Impedance Is often used 
to analyze transient response. When a switching regulator Is subjected to 
a small step-change in load, the output voltage normally varies only 
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Fig* 4 * 5 . 5 . 3 . Theory and measurement of audlosusceptiblllty characteristics. 




•lightly. If the duty ratio la regarded •• a conetent during the load 
tranelant, the linear average model that characterlsee the converter 
at a given quleecent point remains valid. (2A,25t26]. 

The output Impedance of the converter Is defined ae the ratio 

%(■) 

, where l^(s) Is the sinusoidal disturbance at the converter 
o' 

output. The output Impedance, characteristic of an open- loop regulator 
usually has Its maximum value at the output filter resonance frequency. 
This undesirable result can be reduced by effectively designing the 
feedback parameters of the control loop as will be shown In chapter 
six. The block diagram for the output Impedance Is Illustrated In 
Fig. 4. 5.6.1. 

Using Mason's gain formula or a block reduction method the ratio 

V (8) 

« ' ; ■ - V reduces to the following form: 
io(8) 


^ul2 VaC^D2^ ~ VaC^u22^D1 
i^(s) A (A + + ^aC^D2^^ 


(4. 5. 6.1) 


The general form for the output impedance transfer functions can be 
evaluated by substituting equations (4.2.3.12), 4. 3. ...16), (4.3.2.17), 
(4.3.41), and (4. 5. 4. 3) into equation (4. 5. 6.1). 


0 U) 
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L ^ Vp 

L C ■*’d' 
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° A^+ 


M 0 o 

-GP 

1+SR, C, 1 
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(C^+C 3 )R^S 

1+SR^(C^// C 3 ) 
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V N . DL 

-Vac ^ “o 


N_u) ' S 

L P o 


s + 


2D' 

LgC, 


( 4 . 5 . 6 . 2 ) 
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4. 5. 6.1 T«»t VTlflcatlon for Output l»p«danc« 


Pit* 4. 5.6.2 llluttraeas th« Mthod iapl«Mne«d to Moouro tho output 
inpodonco. A voltago sourco vith • one-hundrod oha roolotor In oorioo 
■titutos tho currant aourca aquivalant. Channol A raada tha Input 
dioturbanca currant i^, and channol B aaaouraa tha parturbad raoponoa v^. 

In ganaratlng the theoretical reaponoe , tho synble is uaed for the load 

lapedance Inatead of R|^. defined by the following: 

"load • ('..5.6.3) 

The theoretical and aeaaured reeponae are shown together In Fig. 4.5.6. 3. 
The solid curve is the theoretical output lapedance and the curve marked "x" 
is the measured response. 
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Figure 4.5.6.3Theory and mcusuremenc of Che output impedance 
characteristic . 



4.6. PERFORMANCE EVALUATIONS AND CONTROL DESIGN 


4.6.1 Introduction 

The stability and dynamic performance of the current-injected switching 
regulator are examined. A critical look at the amall-aignal, open-loop, and 
closed-loop reaponaes allows for a qualitatjLve and quantitative understanding 

V 

of the control- to-output tranafer function , dc open-loop response 

V 

system open- loop characteristics snd the close- loop behavior 

''1 

Also presented In this chapter are parametric atudles on the behavior of syatem 
responses by varying certain "key" control variables In the feedback loops. 

These studies in turn egnerate information on how the ac and dc feedback- 
loops effect the stability, audiosusceptibility, andother Important design 
features. After collating all Che above mentioned information, some design 
guidelines are set forth for the dc and ac feedback-loop parameters. 


4.6.2 Discussion of Che Control-To-OutpuC Characteristic 


Referring to equation (4. 5. 2. 4), the following expression can be written: 
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(4.6. 2.1) 


Equation (4. 6. 2.1) is reduced in the following manner, assume the second order 
polynomial of equation (4. 6. 2.1) is represented by: 

(s + a) (s +6) ■ s^ + (a + $) • s + oB (4.6. 2. 2) 


If a > > 6 then equation 4. 6. 2. 2 yields the following approximation: 
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(■ a) (• 4- 0) ; <(> OS a0 (4.6. 2. 3) 
With th« ratultlng approxlutlon» th« ••cond order polynomial almply factors 
into the followings 
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Substituting aquations (4. 6. 2. 4) and (4. 6. 2. 5) for the second order polynoainal 
genersteB the simplified version given ly, 
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The justification for the approximation is proven by substituting the 
parameter values of the tested model into equation (4. 6. 1.6) 

A 

091 (s 4- 32.394K)(s - 47.873K) 

■^TT “ (s 4- 49.376K)(s 4- 304.86) (4.6. 2. 7) 


with o"49.376K and 0*304. 86| a > > 0, aquation (4. 6. 2. 6) is a good approxlmatiosji. 
The pole/zero plot for equation (4. 6. 2. 7) is given in Fig. 4. 6. 2.1. For small- 
signal analysis the power stage parameters are given in Chapter 4.5. 

The control-to-output rharacteristic does not include the dc feedback-loop 
parameters and com" .ion network parameters: howiver, it does include the 
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Pig. 4.6. 2.1 'v POLE/ZERO PLOT 




ac feedback loop parameters. The control- to~outpuC tranafer function conslata 
of two poles and two seroa. It la Interaatlng to compare thla tranafer func- 
tion with that developed at Caltech (following a different modeling approach) 
which consists only of a single pole and single zero [20]. Although these two 

transfer functions are similar In the qualitative nature, yet there exlats 
an Important difference which can be Illustrated by examining the asymptotic 
curves of the bode plots. In Fig. A.6.2.2, notice the dip In the gain 
curve around 6KHz before leveling off as the frequency Increases. Thla 
second-order effect could become more pronounced with a greater pole/zero 
separation. A similar effect can be seen In the phase plot of the 
control-to-output expression. This second order effect Is demonstrated 

as the load parameter increases forcing the current to approach the 
discontinuous mode; a noticeable change In the phase response of the 
control-to-output response occurs. Fig. A. 5. 2. 3 illustrates the 

- 15n. With - Aon the theoretical 
and phase as shown in Fig. A. 6. 2. A definitely has a second zero which I 3 
pulling the phase up. This second order effect around six to eight 
kilohertz cannot be explained by a single zero/slngle pole system. 

From the above discussion, one can conclude that the first-order 

f 

approximation may be valid for a given range of parameter values. 

It however, cculd produce gross error under other parameter values 
or operating conditions. 

A. 6. 3 DC Open Loop Characteristic and Compensation Network 

The dc open-loop response Is simply the control-to-output response 
multiplied by the dc loop gain (Fjjq) . The model for Fjj^ given by 
equation (A. 3. 2. 16) is a general Integral-plus-lead/lag network. 


control-to-output response for R^^ 
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Illustrated In Figures 4.6. 3.1 (a) ond (b) are two comtonly 
used integral'>plus>lead/lag networks derived Iron Pig. 4. 3. 1.1. 

Fig. 4.6. 3.1 (a) and (b) are qulvalent dc loop conpensatlon networks. 
The transfer function for compensation network A, as derived from the 
generalized expression* Is 


DC 


1 sR^C^ 


1 — f . . — ■ 
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(4.6. 3.1) 
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Where (Ip - ( 1 + H3 + p.) 

and compensation network B can be Coddled by the following expression: 
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(4.6. 3. 2) 
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GP ^ ^ 

With network A the lead/lag corner frequencies are Independently 
adjusted by changing the capacitance values Ci and C^ ; in network B the 
same is accomplished by varying the resistive values ^ and R^. Neither 
network has any particular advantage over the other. For the discussion of 
the dc open loop characteristics network A Is chosen. 

Compensation network A has a slngle-zero/two-pole transfer characteristic 
which provides an integral plus lead/lag compensation. Fig. 4. 6 . 3. 2 shows 
the gain response of network A with respect to parameter variations in 

resistor R^ and capacitors Ci and C^. Capacitor Ci alters the lead- 
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corner frequency. Larger values of generate a lower croeaover frequency as 
Indicated by the arrow In Fig. A. 6. 3. 2. Increasing forces the mid- frequency 
gain to increase while decreaalng values of reduces the lag croieover 
frequency. Similarly, increasing Gp (or reducing R^) causes the overall dc 
loop gain to increase. The dc loop response can be reshaped to any desirable 
characteristic by varying these parameters. The relationship of the open dc loop 
characteristic to the closed-loop system response will be discussed in 
section A. 6. 6. 

Using equation (A. 6. 3.1) as the expression for the dc open-loop 

response described by equation (A. 5. 3. 2) is simplified as follows: 

DL 
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(Rj,C 8 + IHR^C^s + 1)(^ 8 - 1) 


8 (R^C^s + 1) 
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2/1 p I>+1n 
“o <1 ♦ 


(A, 6. 3. 3) 


F 0) ^GPR^N 

where CONST - - 

The pole/zero plot for the DC open loop transfer function with the converter 
parameter values given in section A. 5.1 is shown below in Fig. A. 6. 3. 3. 

A. 6. A System Stability and Its Interpretation 

As discussed In section A, 5. A, the true system open-loop response 'cannot 
be experimentally measured, but the dc open loop response can be measured. 

In Caltech's modeling approach [20], the ac loop is lumped into the power stage 
and a small-signal system model is derived which consists of a "new" power stage 
and a dc feedback loop as shown in Fig. A. 6. A. 1. While the dc feedback loop re- 
mains in its original fora.', the "new" power stage in effect incorporates both 
the original power stage and the current-injected loop. 
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Caltech's modeling approach was legitlioata and waa expariiBentally 
verified. Nevertheless* the Interpretation of the model has to be 
exercised with great caution. It is obvious that the open loop response 
of Fig. 4.6.4. 1 is merely the dc open loop characteristics of the original 
system and not the true open-loop response since the ac loop is embedded 
in the "new" power stage model. The effect of ooenlng the ac loop cannot be 
examined using their model. To illustrate, Fig. 4. 6. 4. 2 demsnitrates a two-loop 
system illustrated in a general block diagram fashion. assume that the 

loop containing is the ac feedback loop and the loop containing H 2 Is the 
dc feedback loop. The following transfer functions can be easily derived. 


Close loop gain: 

''0 , 

2 

(4.6.4. 1) 


i + G (H^ + H2> 

Open loog gain 
at A: 

«A- 

G (H^ + H 2 ) 

(4.6.4. 2) 

Open loop gain 
at B: 


G H 2 / (1 + G Hj) 

(4.6. 4. 3) 


The system is unstable when the following Nyguist stability crlterlo'it 
is satisfied: 

G(Hj^ + H 2 ) - -1 (4.6. 4. 4) 

The left-hand side of (4. 6. 4. 4) is identical to the open-loop gain at 
point A. Examining equation (4. 6. 4. 2), provides both the condlton 

for stability, and information concerning the relative stability of the 
system (the phase margin and the gain margin). Examining the open- 
dc loop characteristic Gg, only when Gg - -1 is the condition for in- 
stability as given in (4. 6. 4. 4) revealed. Therefore, only at the point 
of stability is the following expression valid: 
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If, howsver, one employed the open-wc-lcop gain to Investigate the 
relative stability of the system (gain margin and phase margin) the conclu- 
Slone could be very misleading. Tlie point becomes obvious If one examines 
the dissimilarity of the open-dc-loop response vs. the system open loop 
response as Illustrated In Fig. 4. 5. 3. 3 and Fig. 4. 5. 4. 2, respectively. 

In conclusion, although the system open loop response could not be 
measured experimentally. It provides the only proper way of measuring the 
relative stability of the system via Bode analysis. It should be noted that, 
the characteristic equation for system open loop and open*~dc loop are, 
however, the same l.e. 

1 + G - 0 (4.6. 4. 5) 

A 

1 + G- - 0 (4.6. 4.6) 

D 

Equation (4. 6. 4. 5) and (4. 6. 4. 6) are Identical. Therefore the stability could 
be examined via characteristic roots or elgnevalues using either (4. 6. 4. 5) 
or (4. 6. 4. 6). 

4.6.5 System Open Loop Response 

The system open loop characteristics can be examined by the small-slgnal 
model, but experimental verification of that response is not permissible with 

the control hardware currently used. The problem is due to the Inability 
to retrieve the ac-modulatlon-slgnal information from a "logic level" 

signal produced by the output of the comparator in the error processor 
with currently available instrumentation. The comparator output, 
determined by the control Inputs (switching current and the dc error 
voltage), is the only location where both dc and ac loops can be opened 
to examine the system response. To check the system open-loop character- 
istics, the various gain blocks formulating the analytical expression 
are examined via the control-to-output response and the dc open loop 
behavior. The e..pression for the system open loop response G_. as given 
by equation (4. 6. 5.1) utilizes the compensation network "A: represented 
by equation (4. 6. 3.1). 
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(4.6. 5.1) 


A unique characteristic of equation (4. 6. 5.1) la that the number of zeroea 
la always equal to the number of poles regardless of the form of compen- 
sation network used. Equation (4.6. 5.1) has four poles and four zeroes; 
but If the order of the system changes due to the addition of a pole or a 
zero In the compensation network, the open loop transfer function 
will carry five poles and five zeroes since the second term of (4. 6. 5.1), 


F F 
AC D2 


, has an equal number of zeroes and poles (two zeroes and two poles). 


This second order expression Is determined solemnly by the power stage 
and is Independent of the form of the dc feedback loop. Equation (4. 6. 5.1) 
can be expressed In the following form: 
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In the general form, equation (4. 6. 5. 2) la too difficult to factorise. 
Possibly If bounds are given on parameter valuea, aome simplification can be 
done. Due to the complexity of the fourth order expreaitflon and not to 
take from the generality, a parametlc study of the system open-loop response 
along with other system behaviors will be discussed In the next section. 

Another salient feature of the system open-loop characteristic Is 
the gain which approaches a constant value at high frequency. A simplified 
expression for the constant gain at high frequency can be attained. From 
equation (4. 6. 5. 2) assuming s .*> «, equation (4. 6. 5. 3) yields: 


OL 




*^AC 


(4.6. 5. 3) 


In reality, however, the gain is expected to decrease at high frequency. 
This Is because the pulse modulator gain can no longer be regarded constant 
at such high frequencies. The analytical model presented here falls to 
take the high frequency effect of the pulsd modulator Into account. Due to 
the complicated describing function modeling employed here, the pulse 
modulation Is simply represented by a constant gain. 
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4.6.6 Parametric Evaluations 


A parametric study la performed on the small-signal behavior of the 

current-injected converter. The purpose cf the parametric examination is 

to key cn particular salient features of the control characteristics and 

design strategy for the curre.it- Injected control. Due to the complexity 

of each small-signal models, especially the audiosusceptiblllty and the 

output Impedance characteristics, the analytical expressions could nut be 

simplified without a loss In generality. For this reason, parametric 

study by ways of changing certain key design parameters including the 

T 

duty-cycle variation, , ac loop-gain, dc loop-gain via GP and 

the dc loop-gain via are applied to the small-signal expressions. These 
results can be used to establish design guidelines for the control loops 
In order to optimize the regulator performances. The small-signal 
characteristics examined are the following: control-to-output response, 
dc open-loop characteristic, system open-loop characteristic 

audiosusceptiblllty '^®/v^; and output Impedance, z^. 


4. 6. 6.1 Duty-cycle variation 

Figure 4. 6. 6.1 and 4. 6. 6. 2 demonstrate the control-to-output anc dc 
open-loop response with respect to changes in the steady-state duty-cycle 
ratio. As expected, an increase in the duty-cycle ratio forces the 
control- to-output gain and the dc open-loop gain downward. 

As shown in Fig. 4. 6. 6. 3 the duty-cycle has no drastic effect on the 
system open- loop phase delay. Therefore, a wide variation in duty cycle 
operation should not disturb the stability of the system. The duty-cycle 
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does have a more appreciable effect on the high frequency gain of the 
system open- loop reaponse. The larger duty-cycle forces the crossover- 
frequency to a higher value. Caution roust be exercised when examining the 
high frequency behavior uince the accuracy of the model degenerates when 
the modulation frequency approaches one-half the switching frequency. 

The quality of the system's audlosusceptiblllty demonstrates the 
converter's ability to reject Input noise. Figure 4. 6. 6. 4 shows the lover 
the duty-c'^cle the better the audlosusceptiblllty. Intuitively thlo Is 
understandable, because a low< r duty-cycle allows a smaller percentage 
of the total Input noise to propagate to the output. 

Similar effect is shown when one examines the output Impedance 
characteristic as a function of duty-cycle. The output Impedance reduces 
with the lowering of the duty-cycle as shovm In Fig. 4.6. 6.5. 

4 . 6 . 6 . 2 Variation in ac loop gain 

The ac-loop-gain block is comprised of a resistor, and a turn 

I 

ratio n, a current transformer in series with the power switch. F^^, ■ ^SW/ ' 
In Fig. 4. 6. 6. 6 the ac loop gain, forces the control- to-output gain to 

reduce for larger values in Rg„/n« With either lower or higher F^^ the 
effect of two-pole/two-zero characteristic becomes more pronounced. 

Examining the phase curve for F.„ ■ 3.0, one notices the phase curve rises 
after approaching -90®. This phenomenon will not be observed In the 
simplified single-zer ^/single-pole approximation with the present model, 
a second zero in the left-half-plane causes the phase to rise after 
reaching -90® cud then the second pole cancels the LHP zero's effect 
allowing the phase to continue downward to -180®. 

Increasing the ac loop gain decreases the dc open-loop gain as shown 
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in Fig. A. 6.6. 7. Th« higher dc-open-loop gain, in general* reaulta to 
better closed- loop responses such as audiosusceptibility and output 
impedance as will be discussed later. 

Illustrated in Fig. A. 6. 6. 6 is the effect of the ac-loop gain to the 
system open- loop characteristics. 

-:«2 ac loop gain F^^ has n very dominating effect on the 
system open-loop characteristic at high frequencies and virtually no 
effect at frequencies below the output filter resonant frequency. By 
decreasing the ac feedback two of the systeti open- loop zeros move from the 
LHP toward the imaginary axis. Further reduction in the ac gain results 
in a complex conjugate pair of zeros as exemplified by the second-order 
effect of the gain curve with F^^^ ■ 0.3. Continuing even further the 
conjugate pair changes into two positive zeros. The open-loop phase 
drastically changes when the zeros migrate to the right-half-plane. This 
detrimental phase characteristic which approaches -360" has a very in- 
desirable effect on the system stability. With positive zeros adding to 
the phase delay the phase margin is reduced to a negative value which 
results to an unstable system. 

The audiosusceptibility characteristic as a function of the ac loop 
gain is plotted in Fig. A. 6. 6. 9. The ac feedback loop provides an excellent 
parameter to adjust the low-frequency response of the audiosusceptibility. 
Since the passive filters in the regulator generally can provide adequate 
attenuation of disturbances at higher frequencies, the lower fr<$quency 
range within, say, zero to ten-times the output filter resonant frequency 
is the frequency range to control. For this reason the ac-loop gain will 
play a significant role in optimizing the regulators ability to attenuate a 
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small-signal sinusoidal disturbance propagating from the regulator input 
to its output. 

For smaller values of the audio response reduces particularly 

at low frequencies. The ac feedback gain can be reduced to the point 
where second order peaking effect appears around 5-6 KHz. The second-order 
peaking effect is caused by the emergence of two poles of the closed-loop 
response into a complex-conjugate pair. This condition causes severe 
degradation in the audiosusceptibility performance of the system. As 
discussed earlier, the system becomes unstable under this condition. The 
optimal condition for the audiosusceptibility response occurs when F^^ *0.6 
prior to the occurrence of the two zeros forming a complex conjugate pair. 
Similar effect is shown when the output impedance characteristic is plotted 
as a function of F^^ in Fig. A. 6. 6. 10. 

A. 6. 6. 3 DC loop gain variation via GP 

CP contains the dc gain resistor, R^, as shown in tlic following 
expression; 

CP R— 5 

As Rj decreases, GP Increases, forcing the dc-open-loop gain upward as 
demonstrated in Fig. A. 6. 6.11. The dc open-loop phase is not effected. 

Variations in CP have a significant effect on the low frequency end 
of the system open-loop gain as shown in Fig. A. 6. 6. 12. As GP approaches 
its theoretical limit of one, the system gain Increases. The second-order 
peaking effect occurred around 30KHz exceeds the theoretical bound of the 
model, since the model is only good up to half of the switching frequency 
(AOKHz/2). Examining the phase characteristics of the system open- loop 
in Fig. A. 6. 6. 12 no zero migration to the right-half plane is noticed as 
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mentioned previously with a variation In the ac-loop gain. 

Approaching the theoretical limit of GP Improves the audlosusceptlblllty. 
Again GP has a notable effect on the low frequency section of the audio. 

Fig. 4.6.6.13 shows u 25 db Improvement over the range of GP. Therefore by 
decreasing the dc gain resistor, the attenuation of noise from the Input 
to the output Is Improved. 

Fig. 4. 6. 6. 14 shows the output Impedance characteristics as a function 
of GP. 

4. 6. 6. 4 DC -loop gain variation via 

The dc-loop gain can also be varied by changing the gain feedback 
resistor, R^. But when is varied the corner frequencies of the pole and 
zero produced by the transfer function, F^^, changes. Therefore, the 
parametric study of R^ has been constrained to only take in the gain effect 
and not a shift the corner frequencies. This is done by holding the 
expressions and constant for any change in R^. Two sets of data 

are provided for variations in R^. The first set varies R^ from lOKG to 
150KG while the second set examines variations of R^ from lOOKD to 2MG. 

Figure 4.6.6.15(a) and (b) demonstrate for an Increase In R^ the de- 
epen- loop gain increases. The phase does not change because the constraints 
set previously inhibit the corner frequency of the compensation network 
from changing. 

From Fig. 4.6.6. 16(a) the variation in R^ has no degrading effect on 
the phase and affects the system open-loop gain only in the low frequency 
range. Apparently, the variation from lOK to 150K of R^ has a negligible 
effect on the system's stability. As R^ is increased further (R^ IMD), 
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two zeros migrate to the RHP forming a complex-conjugate pair and forcing 
the phase to approach -360* instead of 0* as shown in Fig. 4.6.6.16(b). 
This drastic change of the phase has detrimental effect on the system 
stability. In fact» for large R^, Che branch formed by R^Cj^ can be re- 
garded as open circuit. Figure 4.6.6.17(a) shows Che switch current wave- 
form of the stable converter employing the set of circuit parameter 
values given in Chapter 5, page 53. Figure 4.6.6.17(b) illustrates the 
switch current waveform of the unstable system when Che branch R^C^^ is 
opened . 

The effect of R^ to the audiosuscepCibllity curves is given in 
Fig. 4.6.6.18(a) and (b). Notice as R^ increases Che audiosusccptibility 
steadily improves until the two positive zeros are felt causing second- 
order peaking in the high frequency range. The complex zeros in die RHP 
eventually turn into real zeros. As the peaking effect is 

no longer shown and the audiosusceptibility is greatly improved. Never- 
theless. Che system becomes unstable due to the excessive phase delay. 
Similar characteristics are displayed for Che output Impedance curves 
as shown in Fig. 4.6.6.19(a) and (b) . 
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3PENLP PHASE (DEG) 

-360.03 -3C0.00 -2<i0.00 -180.00 -120.00 -6C.OO 
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JIO.OO -90.00 -70.00 -SO. 00 



Fig. 4. 6. 6. 4 Audiosusceptibility gain and phase characteristics for 
D - 0.25, 0.3, 0.4, 0.5, 0.6, 0.65, 0.7, 0.75. 
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6.5 Output Impedance characteristics for 

D - 0.25, 0.3, 0.4, 0.5, 0.6, 0.65, 0.7 
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Fig. 4. 6. 6. 8 System open-loop gain and phase characteristics for 
F._ - 0.05, 0.1, 0.3, 0.6, 1.0, 1.5, 2.5, 3.0. 
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Fig. 4. 6. 6. 9 Audiosusceptibility gain response for • 0.05, 0.1, 
0.3, 0.6, 1.0, 1.5, 2.5, 3.0 
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Fig. A. 6. 6. 12 System open-loop gain and phase characteristics for 
GP 0.05, 0.075, 0.1, 0.3, 0.5, 0.7, 0.9, 0.95. 
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Fig. 4.6.6.13 Audiosusceptibility characteristics for GP - 
0.075, 0.1, 0.3, 0.5, 0.7, 0.9, 0.95 
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Fig. 4.6.6.15(b) DC-open-loop gain and phase characteristics for 
- lOOK, 200K, 300K, 400K, 500K, 700K, 900K, 2M0. 
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Fig. 4.6.6.16(b) System open-loop gain and phase characteristics for 

- lOOK, 200K, 300K, 400K, 500K, 700K, 900K, 2MD. 
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6.6.18(a) Audloausccptiblllty characteristics for > lOK, 30K, 
SOK, 70K, 90K. IIOK, 130K, 150K. The product and 

R^C^ remain constant. 
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Fig* 4.6.6.19(a) Output Impedance characteristics for 

- lOK, 30K, 50K, 70K, 90K, IlOK, 130K, 150K. 
The product R.C. remain constant. 
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6.6.19(b) Output impedance characteristics for 

• lOOK, 200K, 300K, 400K, 500K, TOOK, 90CK, 2Mfl. 
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4.7 CONCLUSIONS 


Modeling and Analysis of a buck/booat regulator eaploylng the current- 
injected (current-programmed) control lo presented. The objective la to 
examine the small-signal dynamic performance of the switching regulator 
Including; 

a Open-loop characteristics and system stability, 
a Closed-loop characteristics and audlosusceptlblllty. 
a Response due to load disturbances and output Impedance 
characteristics. 

a Features of the current-injected control and design guidelines. 

In order to gain Insights of the performance characteristics of 
various parts of the switching regulator employing the current-injected 
control, the regulator Is modeled according to the three basic functional 
blocks: power stage, error processor, and duty-cycle pulse modulator. 

The power stage model presented in Chapter 4.2 consists of three inputs: 

A A 

disturbances from the line v^„ the loa.l i^> and the duty-cycle control 

A A ^ 

loop d; and two outputs: the output voltage v^, and the switch current 1^. 

The error processor (EP) model as presented in Chapter A. 3 contains the dc 
feedback, current-injected loop and the compensation network. The EP 
processes and compensates the signals from both the dc and ac feedback 
loops and provides the necessary error signal to the duty cycle pulse 
modulacor (PM). The PM takes the analog error signal from the EP and 
converts the error signal Into a series of pulse-width modulated duty- 
cycle signals. Presented In Chapter 4.4 Is the low-frequency model of the 
duty-cycle pulse modulator employing the describing function techniques. 
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In Chapter 4. S» Che anall-algnal nodel for Che power stage, the error 
proceoaor and pulse modulator, derived from Chapters 4. 2, 4. 3 and 4.4, raspec 
tlvely, are assembled. Employing this amall signal modal, the following 
open- and closed-loop characteristics are examined: 

• The control-to-output characteristic 

• The dc open-loop characteristic. (The dc feedback loop Is opened, 
while Che current-injected loop Is Intacted.) 

• The system open-loop characteristics. (Both the dc feedback loop 
and Che current-injected loop are opened.) 

• The closed-loop Input-to-ouCput response - audlosusceptlblllty. 

9 The closed-loop response due to a load disturbance - output 

Impedance characteristic. 

Experimental verifications are also provided for the above described 
analytical models with good correlations. 

In Chapter 4 .6, various regulator open-loop and closed-loop performance 
characteristics are examined In further detail. Effects of various key 
control parameters to the regulator performance characteristics are In- 
vestigated. Features of Che current-injected control are discussed. In 
particular, comparisons are made b<^tween the results derived In the 
current modeling and analysis efforts and Che earlier works presented by 
Caltech Cwhlch followed a different modeling approach) . 

Many interesting characteristics are unveiled upon careful exami- 
nation of the small-signal model derived In this report. Among those 
findings. Che following ones are part cularly noteworthy: 

(1) The ac feedback loop Is imbedded In the control- to-out) ’>t <.?»>:•• 
fer function which exhibits two-pole and two-zeto chay - 
tic Instead of a single-pole single-zero. However, the two-pole 
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two-zero transfer function can be approxlnated by a single-pole 
single-zero transfer function without encounter Ing gross errors 
under most circumstances. 

(2) Vlhen the system open-loop characteristic Is examined with 
both dc and ac loops opened, the two-pole two-zero (or single- 
pole single-zero approximation) characteristic ceases to 
exist. In fact the system open-loop transfer function exhibits 
astonishingly different characteristics. Very little resem- 
blance can be found when one compares the open-loop character- 
istic G.J, with the open-dc-loop 'characteristic 

(3) The stability characteristic for the current-injected control is 
quite unique comparing to other types of multi-loop control 
systems [2 2i 23] . The stability margin could extend beyond 90° 
with very high cross-over frequency. The stability of 

the system is particularly sensitive to certain control para- 
meter values. A positive zero could be observed In the open- 
loop characteristic if the control paraueters are not properly 


selected. 


5.0 MAPPS DEMONSTRATION PROBLEM FOR VSTOL EMERGENCY POWER SYSTEM 


5.1 INTRODUCTION 

The VSTOL Emergency Power System Is shown In Figure 5.1. A NIrCd 
battery source with an effective series resistance Rg Is processed by 
a boost converter to provide the bus power for emergency use when the 
generator failed. The battery voltage level depends on the number of 
Nl-Cd cells connected In series, with the number yet to be determined. 

The boost-converter power circuit consists of Inductance L with winding 
resistance RL, power switch Q, power diode D, and output capacitor C 
having an equivalent series resistance (ESR) RC. Since the losses In 
the converter are supplied from the battery source, and since the con- 
verter packaging weight (heat sink Included) Increases with converter 
losses, it follows that the combined battery and converter mechanical- 
structure weight becomes heavier If more converter loss Is allowed for a 
given output power. On the other hand, the converter component weights, 
I.e., those of magnetics and capacitors, tend to diminish with more allow- 
able losses. Consequently, for a given output power, there must exist an 
optimum converter efficiency at which the combined system weight Including 
battery, packaging, and converter components, is at Its minimum. The 
essence of the optimization Is to Identify the battery voltage level 
along with the detailed boost converter design so that the total system 
weight of the battery and the packaged converter can be minimized. 

The optimization effort starts with the Identification of key oper- 
ating waveforms of the boost converter as functions of line/load conditions, 
from which the voltage and current expressions are derived for all power- 
handling components. These expressions are then used to formulate component 
losses as well as other converter design constraints Including Input EMI, 
output ripple, and the proper design of the Inductor. The optimization 
objective of minimizing the system weight is then defined analytically. 


At this Juncture, two approaches become possible. One Is to rely on 
an established nonlinear programming routine such as the SUMT (Sequen- 
tial Unconstrained Minimization Technique) to numerically seek the 
optimum design and the attendant minimum system weight. Two basic 
variables In this approach are the Input voltage to, and the switch- 
ing frequency of, the converter. They are treated as unknowns, to 
be determined along with other variables In a single, large-scale 
numerical optimization. A second approach, less ambitious and chuv 
perhaps more practical, breaks the single large-scale optimization 
Into many smaller optimizations. The two aforementioned variables are 
given as known values, thereby greatly simplifying the problem and 
allowing a closed-form optimum VSTOL design for all other variables. 
Without resorting to nonlinear programming, optimum designs for other 
sets of known Input voltages and switching frequencies can then be 
similarly calculated, from which the true global optimum corresponding 
to a specific set can be Identified among calculated optimum designs 
by mapping the data of optimum designs for different sets of voltages 
and frequencies. A decision was made to adopt the second approach to 
the VSTOL system optimization, with the purpose of establishing an 
alternate means to nonlinear programnlng In achieving optimization. 

The chosen approach was successfully Implemented, and closed-form 

optimum system designs for any given set of Input voltage and switch- 
ing frequency were obtained. A computer program was generated to 
perform the straightforward numerical calculations of the analytically- 
derived closed-form solutions. A 3000-Watt, 0.5-hour emergency 
application was used as a demonstration example The closed-form 
solutions, the computer program, tpe calculated results and graphical 
representations, the discussion and the conclu'^^lon regarding the 
VSTOL emergency power system, are presented. 

5.2. CURRENT AND VOLTAGE WAVEFORMS 

The relevant boost converter waveforms are given In Figure 5. 2. The 
power switch operates with an on-time TN and an off-time TF. The 
switching between the two time Intervals are exaggerated In length for 
both switch Q and diode D to Illustrate for the switch fng-l os s model- 
ing to be formulated In a later section. 

The power-switch current IQ and diode current ID are shown In Figures 
5.2Aand 5.2B. The midpoint of each pulse current has a level of If. The 
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Figure 5.1 VSTOL Emergency Power System 
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average diode current Is Identical to the output load current lo. The com- 
posite waveform of Figure 5-2A and 5-2B Is given in Figure 5-2C which re- 
resents the Inductor current 1L as well as the converter Input current. 

The average converter Input current Is therefore II. The difference between 
ID and becomes the capacitor current 1C, as Is shown In Figure 5-2D. 

The switch and diode voltage waveforms are given In Figures 5-2E and 5-2F. 

5.3. BASIC VOLTAGE AND CURRENT EXPRESSIONS 
Let the following designations be made: 

D: duty cycle defined as (TM)/T 

d 1 : peak to peak ac component In the Input current 
F: switching frequency 
ICt output capacitor current 
ID: diode current 
1L: Input (or Inductor) current 

IQ: power-switch current 

II: the midpoint pulse current shown fn Figure 5-2A 

L: Inductance In boost converter 

PL: total loss In the converter 

PO: output power to the load 
TN; on- time of the power switch 
VI: Input voltage to the converter 

VD: output voltage to the load 

The following expressions can be formulated: 

(ID)ave « —§ 5 - (5-1) 


II - 


PO 

V0*O-D) 


(5-2) 
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(iL!n» . I [ ♦ C ^ 

dOrms . [(It)^*0-D) - (PO/VO)^]”’* • “TO *tlV^° ^ <*. 


Since the total Input power Cvr) (11) fs equal to the suf of output Po 
and loss PL, 


PO + PL 


P0*VI 

V0*(T-D) 


From which one obtains: 


(5. 


D » 1 - 


P0*VI 


(5. 


1 - 0 


P0*VI 


( 5 -' 


Substituting (5-9) and (5-10) Into (5-2) through (5-7) one has: 



(1L)«y« • 


(S-13) 


n * “cSt" * c ’ - -Ti^S^krr ^ 

(iL)m» . [(1l)J„ ♦ 

(ic)n« . -58- [ -yj2i2tl . ,!“•* 

5.4. FORMULATION OF CONVERTER LOSSES 

Let: A ■ Cross sectfonal area of fndactor core 

AC. ■ Winding area per torn of tndoctor 

BOC ■ Intended maxfmam dc operating flax density level 
ac flux-density excursion in the inductor core 

C ■ Capacitance at output 

F ■ Switching frequency 

FC ■ Pitch factor of inductor wfnding 

FW ■ Inductor core fill factor 

L • Inductance 

N > Number of tlioes* ured in inductor 
RC ■ ESR of C 

RHO “ Resistivity of iliductor wthdfng 
RI • Peak-peaR allowable ripple at tRe output 
TOF ■ Fall time of diode duriiig switcfting 

TOR • Rise tUne of diode during switclfing 

TSF ■Fall time of power switch during switching 

TSR ■ Rise time of power switch during swftcfting 


(5*14) 

(S-IS) 

(5-16) 


- 283 - 


VB " Forward base-afflitter voUaga of transistor 
VD > Forward drop of tha dioda 
VS ■ Saturation voltaga drop of tha powar switch 
Z • Maan langth of Inductor core 


With these descriptions defined, the following losses can be formulated. 


5.4.1 INDUCTOR LOSS 


The Inductor loss consists of copper loss and Iron loss. Based on 
the assumption of a toroid Inductor and a square cross section for the core, 
the copper loss becomes: 


PLC ■ (II) *1nductor winding resistance 
, 0.5 


[ - ] * 4 *RH0 * 


FC*N*A 


(5-17) 


The Iron loss per switching cycle can be expressed by the area of the 
minor BH loop. Such a characteristic Is shown In Figure 5-3. The flux 
excursion, 0, can be expressed as: 

0 • (V0-VI)*(1-D)/(N*F) 



- BAC * A 


(5-18) 


where BAC - 

The ampere-turns N1 corresponding to the Btt loop-width. H of the Inductor 
core can be expressed as Nf-Hz, where IT Is normally expressed In oersteds. 

Since H Is a nonlinear fui.ctlon of the switching frequency F, It can be 
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Since H Is a nonlinear function of the switching frequency F, it con be 
related to F as the following: 

H ■ 0.089 * F®*® (For Permalloy powder core only) (5-20) 


where H and F are In amp- turns/meter and hertz, resnectively. The loop 
width is 2 Ni , or 


2Ni ■ 2*H*Z 


(5-21) 


The iron loss becomes: 


PLI ■ 2*(Ni)*0*F ■ 2*BAC*H*A*Z*F 


(5-22) 


where (BAG) and H are given in 5-19 and 5-20, respectively. 
The total Inductor loss is therefore: 


PL -PLC PLI (5-23) 

where PLC and PLI are respectively given in (5-17) , >2). 

5.4.2 POWER SWITCH 

Losses in the transistor include the conv . less, the base loss, and 
the switching loss. The conduction loss PTC is: 




PTC ■ I1*VS*0 



1 


t* 


4 



I 


I 


I 

♦ 


t 


PTC . I1*VS*0 . S2^. VS • c 1 - 
Assuming a 10 to 1 base current drive* the base loss Is: 


PTB ■ 0.1 * II * D * VB - 0.1 * * VBE * C 1- 3 ( 


PO * VI 


The switching losses during turn-on and turn-off can be formulated with 
aid of Figure 5-2D. During turn-on* the voltage across the power switch 
decreases from (Vo * VD) to VS* and the current through the switch In- 
creases from zero to [ ( P04-PL)/VI - (d1/2) Assuming the rate of 
voltage or current change during the switching Interval TSR Is constant* 
it can be shown easily that the turn-on switching loss Is: 


PTR - (1/6) * (VO+VD-VS) * ( -SO+PL *^5^ *p 


During turn-off, the voltage across the power swttch Increases from VS 
to (Vo + VOI, and the current* decreases from [ (P0+PL)/VI ♦ d1/2 ] 
to zero. Therefore, during switching Interval TSF, one has: 


PTF • (1/6) * (VO+VD-VS) * (-2Jy^ + ^^) * TSF * F 

Realizing the opposite signs associated with the d1/2 term In (5-26) and 
(5-27),one makes the simplifying assooptlon of equal TSR and TSF* and the 
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sum of switching losses becomes: 


PTR ♦ PTF ■ (1/6)*(V(HV0-VS)*^§{£Ss— *(TSR+TSF)*F (5-28) 


By sunning up (5-24) (5-25) and (5-28) the total loss In the power switch 
becomes : 

PT - PTC + PTB ♦ PTR PTF 

■”^vr^ [l-75*pyglisq-] *(vs+0.1*VBE)+(V0+VD-VS)*(TSR+TSF)*F/6^g_ggj 

5.4.3 DIODE LOSS 

The diode conduction loss Is: 

PDC ■ II * VD * (1-0) ■ PO * VO/VO (5..30) 

Essentially following the same derivation for the pov»er-sw1tch switching 
loss, the diode switching loss becomes: 


PDR + PDF - (l/6)*(V0+VD-VS)*-2§y^ * (TDR+TDF)*F 

The total diode loss Is: 

PD - PDC + PDR + PDF 

- - + (V0+VD-VS)-£§y^ *(TDR+TDF)*F/6 


(5-31) 


(5-32) 
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5.4.4 CAPACITOR LOSS 

Tht capacitor loss PC Is: 

K • (iC)^ • RC . (-^)* * RC • [ - 1] 


(S' 33) 


4 


I 


I 

P 

•f 

r 


5.4.5 OTHER CONVERTER LOSSES 

The total convertor loss In the control circuit* housekeeping circuit* 
wires and connectors Is assumed to be 1.SX of the total output power* I.e. 

'OT • 0.015 * PO (5-34) 


5.4.6 TOTAL CONVERTER LOSS 

The total converter loss can be obtained by summing up (5-17), (5-22), 
(5-29), (5-32), (5-33) and (5-34): 

PL - (- ^ ■ — )^ * 4 * RHO * A — 

* 2 *Z 


[1- ggs^JJJpq-]*(VS+0.1*VBE)+(V0+VD-VS)*(TSR+TSF)*F/6 
+ + (VO+VD-VS) * * (TDR+TDF) *F/6 


. (.50., 2 . t « 


♦0.015* PO (5-35) 

With PO much greater than PL* very little error wdl be induced if the 
following simplifying approximations are made: 

P0/(P0*PL) *0.93 

2 2 (5-36) 

(PO+PL) - 1.156 *P0 


Substituting this* aDP'/oxIiMtlons 1nto(5-17)and(S-l9) •q.(5-32)btC0Mis: 

PL ■ PL1/PL2 
whtrt: 

PLl ■ 4.624*RHO* )*(VS-K).1*VB) 

♦ *(V(HVO-VS)*(TSR+TSF'MOR+TOF)*F ♦ 

+ 0.015 * PO 
*A * Z 


- * (TSR+TSF+TDR+TDF) * F 


+ (“5§-) * RC * (-Jf- - 1) 

. O.m ^ 

p,2 . 1 . JSi™ . ^ 


5.5. FORMULATION OF CONSTRAINTS 
The following constraints are to be observed: 

e The Input current ac component mast be limited to below 
a certain peak*to-peak amplitude. 

• The Inductor should not operate In the saturation region. 

• The output voltage ripple must be limited to below a cer* 
tain peak>to-peak amplitude. 

• The core must have sufficient window area to accomodate all 
the windings. 
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Thtst constraints art formuiatad as follows: 


5.5.1 PEAK-TO»PEAK INPUT AC COHPO^IENT 



Tht poak-to-paak ac componant has bntn txprtssed as d1 In tq/S-14). 
Lot the maximum amplitude be specified as TAC. then. 


IM > * [1 - (P0^)*V0 ] 


S.S.2 BELOW INDUCTOR SATURATION FLUX DENStTT 


This constraint specifies that at the peak Inductor current, the 
corresponding flux density In the core should be below a certain pre- 
determined level. (BOC). Consequently. 

N * A - L * IP/BDC > 0 


_ TA . JA A« _ PO+PL .VI * r, PO • VI 1 
IP - II + d1/2 ■ — yyi — + — * [1 - J ' 


PO * VI 


where N and A are the turns and the area of the Inductor core. 


5.5.3 OUTPUT VOLTAGE RIPPLE 

The peak-to-peak output-voltage rfpple Is caused by two components: 
the capacitive component due to the ampere-second prot.'^ssed by C and 
the resistive component due to the ESR RC of C. ft can oe shown that the 
constraint concerning the sum of these two components can be expressed as 



RI i (II ♦ 41/2) * RC + P0*0/(2*V0*C*F) 


RI > 

m 


pg»PL 

■7T" 


• RC ♦ * [1 




Where RI is the peak-to-peak ripple specification. In actual design, 
series-parallel combination of capacitors may be needed in order to meet 
the voltage and ripple requirements. However, the product (RC * C) will 
always remain a constant value regardless of the combination used. 

For example, three different combinations of capacitors are given in 
Figure 5-4. Capacitance C and resistance RC are associated with each in- 
dividual capacitor. It can be easily proved that the equivalent RC 
product for all three connections are identical, and are equal to RC * C. 
Let this constant value be designated G for a given capacitor type, one 
has: 


RT > PO*'PL * G_ * 1 * ri P0*VI T * / VI*G PO \ 

■ T" 2*T" ^ ^ WC ^ (S^do) 


m branchts 



Time Constant 
RC*C 


Time Constant 
RC*C 


Time Constant 
RC*C 


Figure 5-4 Different Capacitor r-r.nectlons with 


identical time constant RC*C. 



5.5.4 SUFFICIENT WINDOW AREA 

Assuming a toroid cor* and a square cross-sectional core area. then, 
for a mean length Z. the window radius will be (Z/(2*3.1416) - A®*®/2). 

The window area becomes 3.1416*[Z/(2*3.1416) - A®*®/2]^ . This 
area* multiplied by window fill factor FR. most be sufficient for all the 
Inductor windings, which has a total area of N*AC. Consequently. 


N *AC < t ♦ FW» ( j4y 


aO.5 ^2 


N * AC^ 


0.5 


Z 


“W 



(5-41) 


5.6. SIMPLIFICATION OF CONSTRAHtTS FOR CLOSED-FORH OPTIWrZATrON 

In equations (5-38) to (5-40), power loss PL Is always associated 
with an output power PO. Admittedly, PL Is a highly complicated function 
of many variables, as Indicated In eq. (5-37). If one seeks a closed- 
form optimization using the aforementioned constraints, the parameter 
PL must not be treated In accordance with (5-37). Rather, and Indeed 
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fortunattly so* llttlt trror will bt Introductd to thtst constraints If 
PL Is givtn a rtasonabla numarlcal valot* and ustd to dariva varlablas 
such as L» C* RC» N, A* ate. Thasa known valuas can than ba 
substitutad 1nto($-37)to raprasant a nora pracisa PL. In assanca* an 
Itaratlon of PL basad on a raasonabla approximation will ba utlllzad for 
tha saka of anhancing closad-form optimization. Nomarically* tha Itarativa 
approximation Is justiflad by considaring that from (S-38) to (5-40)* PL Is 
always assoclatad with PO In tha fom of (PO ♦ PL)/P0, which Is closa to 
unity dua to tha fact that PO Is ganarally much greater than PL. For 
example* a 10X error In assigning a numerical PL would only result In 
about n error for (PO PL)/hO. Consequently* the variables to be thus 
obtained are only subjected to a rather insignificant error* which In turn, 
will Introduce only a relatively small error for the precise expression for 
PL when the various variables are eventually deployed 1n(5-37). Consequently* 
a numerical PL • 0.07 * PO will Be used tn(5-38)to(5-41)1n quest for 
closed-form optimization* and equat1ons(5-38)to(5-40)become: 




(5-42) 


(5-43) 


C > 



1.07*P0*M 

71 — 




(5-44) 
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In eqs. (5-42) to (5-44). the toms lAC. BDC. and RI, Voltage VO. and 
power PO. are specified parameters. Consequently, given an Input voltage 
VI and frequency F. both L. N*A. C» and RC can be determined numerically. 
Substituting these values Into (5-37) reduces It Into the following: 


PL • 


K1*N*A°-® 




* K2 *A *Z ♦ K3 


(5-45) 


where constants K1 . K2. and K3 are: 


K1 


4.624 * RHO * (PO/VI)^ * FC 



(5-46) 


K2 


0.165 


* p0*6 * 


r :.. 


-VI)*VI/(V0*A*N) 


(5-47) 




(-P0 . 5^)*(vSi-0.1*VB)+^5§Y<TSR+TSF+TDR+TDF)*F*(V0+V0-VS)+0.015*P0 

K3 . ^ 


K4 - 1 - I S l^ V jE . ... VO g /^Vi * (tsR+TSF+TDR+TOF)*F - ■ 

5.7. BATTERY WEIGHT 

In order to formulate battery weight, let the following parameters 
be defined: 
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n ‘ Numbtr of Nl-Cd ctlls In strfts within tht bntttry 
RE Total effective resistance Internal to battery 
re : Effective retlstanca of each cell, I.e., RE-n(re) 

T : Anticipated emergency power utilization time 
VB Battery voltage before the Internal drop across RE 
Vk : Individual cell voU»ge, I.e., VB-n(Vk) 

With reference to Fig. 5-1, 


VB ■ VI + * n * re • n *Vk 

" ■ “vr * ^ ^ ^ ] 

The total power drawn from the battery becomes 
PB • PO ♦ PL + ( * n * re 

. PO . PL ♦ 

The voltage required from the Battery is: 

VB . p,/„ . VI 4 . Cl- ]■' 

. VI 4 -24^]-' 

(see eq. (5-56) later) 


(5-50) 


(5-51) 


(5-52) 


- 297 - 


Tht batttry tmptrt'hour ctpaclty rtquirtd Is: 


AH > 


PB * T 

ir*m - 


m • 1.0 for T-1 
m ■ 0.9 for T"0.5 


(5-53) 


whore m represents the reduction of cell cepecity due to more then 1C dischirge. 
Using the General Electric Nl-Cd Handbook as a reference, the Individual 
cell weight Is essentially a linear function of Its ampere-hour (AH) rating. 

Let Wk be the Individual cell weight, then. 


Wk - K * AH 


K*T*PB 


(5-54) 


where K 42 grams/ (arop-hr) for Nl-Cd ceT^s. Comb1n1ng(5-50)4nd (5-54) 
the total battery weight Is therefore: 


WB ■ n * Wk • 


K*T*PB 

TPBT 



(5-55) 
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Sinoj It Is ti«w known thtt rtslstanco (rt) Is Invtrstly proportional to 
tht coll size, end Is within o range of 10 - IS nllllvolts per C rate of 
discharge, one has 


rt« 


0.015*VI*» 

PI 


Substltut1ng(5'S6)to(5«5S)g1ves the total cell weight as: 




( 5 - 56 ) 


The total battery weight Is greater than the cell weight by an amount 
depending on battery packaging, with the difference likely to Increase 
with the number of cells used In series. Modeling the packaging weight by 
by the battery weight becomes: 

„ . . Cl - (5.57) 

5.8. OPTIMIZ/TION OBJECTIVE FUNCTION 

As previously stated, the objective of the optimization Is to minimize 
the weight of the battery-converter system shown In Figure 5-1. The total 
system consists of the following throe weight contributors: 

e Convert*: magnetic and capacitive components 

e Battery 

e Converter mechanical structure Incorporating thermal design 
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5.8.1 CONVERTER COMPONENT WEIGHT 

The converter magnetic weight Is: 

WM • 4 * FC * AC * DC* N * DI *A * 2 (5-51 


where these parameters were all defined In Section 4, with the two terms on 
the right hand side of (5-58) representing copper weight and Iron weight of 
the Inductor, respectively. The converter capacitor weight Is: 

WC - KCAP * C (5-5 


where KCAP is a constant expressed as grams per microfarad of capacitance 
for a specified capacitor type »nd voltage rating at a prescribed temp- 
erature. For this application, foil-tantalum capacitor Is tentatively 
chosen. Taking Into consideration the 400V rating required to pro- 
vide a proper derating for a 270V working dc voltage, the KCAP for this 
application Is taken as 3 grams /microfarad. From (5-58) and (5-59) the 
total converter component weight becomes: 


WCON « 4*FC*AC*DC*N*A°*^ + 01^*2 + KCAP*C 


(5-6i 


5.8.2 Bat tery weight 

The battery weight, WB, was defined In (5-57). 
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5.8.3 CONVERTER MECHANICAL STRUCTURE WEIGHT 


For the first order of approximation, the converter mechanical 
structure weight, taking thermal design Into consideration, Is assumed 
to be directly proportional to the converter output power PO and loss PL. 


WP ■ KMEC*P0 * KTHE*PL 


(5-61) 


KHEC and KTHE are proportionality constants. For this application, KMEC 
Is estimated at 5 grams/watt, and KTHE at 10 grams/watt. 

5.8.4 OPTIMIZATION OBJECTIVE FUNCTION 

The objective function "OBJF" Is the sum of eqs.(5>57) (5>60) and 
(5-61). Analytically, 


OBJF • WB + WCON + WP ■ VfT 

_ _K*T* PB * r, O.Ol5*m*(PO+PL) 

[1 - ‘ J " 


+ 4*FC*AC*DC*N*A°‘® + 0I*A*Z + KCAP*C 


+ KMEC*P0 + KTHE*PL 


(5-62) 


Strictly for convenience, the parameters fn(5-62)are clartfted again as 
follows : 


K Proportionality constant between Individual cell weight 
and Its ampere-hour ••atlng (42 grams/ampere-hour). 
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KTHE: Weight constant for thermal design (10 grams/watt of loss). 

T Anticipated emergency power utilization time (hours). 

PO Required emergency output power (watts). 

PL Converter losses (watts). 

: Cell capacity reduction when discharging more than 1C. where C 

Is equivalent to the ampere-hours rating of the cell. 

For a half-hour discharge application, m- 0.9. 

Vk ; Individual cell voltage (volts, l.OOV for a Nl-Cd cell). 

FC : Winding pitch factor (assumed to be 2.0). 

AC Inductor winding area Csqoare maters). 

DC Conductor density (8900 kg/ro for copper, or, 8,9 x 10 g/m ) 

N Number of turns of the Inductor winding (dimensionless). 

A Cross-sectional area of Inductor core (square meters). 

DI Core density (7800 kg/m^ for Iron, or, 7.8 x 10® g/m^) 

Z .• Mean core length (meters). 

KCAP: Weight constant for capacitors (3 grams /microfarad). 

C : Output capacitance (microfarads). 

KMEC: Weight con.stant for mechanical structure (5.0 grams/watt 

of converter output power) . 

Having clarified these parameters, they can be catagorfzed as follows; 

e Specified system requirements: T, PO, m. 

e Given cell characteristics; K, Vk. 
e Assumed constants: FC, DC, DI’, KCAP, KMEC. 

• Previously calculated parameter: C (see eq. (5-44). 

e Unknowns : PL, AC, A, Z. 


In ampere! 
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Separating (62) Into known and unknown terms • one has: 

02 1 

MT - * KCAP * C ♦ KMEC*PO] 

* o--21|^l^iE2tEki.)-\ KTHE] • PL 

♦ 4*FC*AC*DC*N*A°’® ♦ DI*A*Z 


(5^3) 


WT - [CONST 1] ♦ [CONST 2]*PL + 4*FC*AC*DC*N*A°‘® + DI*A*Z 

(5-64) 

where CONST 1 and CONST 2 are two constants representing, respectively, 
the quantities in the two square brackets on the right-hand side of (63). 

It is recalled that PL of(5-63)was expressed as a function of the other 
unknowns in (5-45). Corabining(5-45)and(5-64)and simplifying the result- 
ant equation, one has: 

WT ■ (CONST 1 + K3 * CONST 2) 

+ CONST 2 * K1 * N *A°*^/AC 

+ 4*FC*DC*AC*N*A°-^ + (DI + K2*A*Z*C0NST 2) 


(5-65) 


where Kl. K2. and K3 were gfven fn (5-46) to (5-49). Equation (5-65) is the 
objective function to be optfmfzed. The unknowns fn (65) are N, A, AC. and Z. 
For simplicity in programming, (PO+PL)/PB in(5.$3)1s treated as unity. 

5.9. REVIEW OF OPTIMIZATION CONSTRAENTS 

It is recalled that fn Sectf on 5.5 there were four constraints. 
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For a given VI and F, one constraint concerning the Input current ripple 
resulted In the numerical identification of Inductance L In (5-42). Another 
constraint regarding the output-voltage ripple made possible the design 
for the capacitance C In (5-44). The two remaining constraints on core 
window and output ripple are listed below for convenience: 


N*AC ,0.5 Z ^ A°*^ 

FW^ ■ FF + 




(repeat of (5-42)) 


N*A > 


L*IP 


“BDT“ 


(repeat of (.5-43)) 


Here, the known parameters are: 

FW Core window fill factor (assumed 0.35). 

2 

BOC Maximum operating flux density (0.35 Weber/meter ). 

L Inductance In Henries as calculated from the equal 

sign of (5-42). 


VI 

PO 

F 

IP 


Converter Input voltage In volts. 

Converter output power In watts. 

Converter switching frequency In Hertz. 

Peak inductor current In amperes defined In eq. (5-39). 


The unknowns are N, A, AC, and Z, which are the same ones In (5-65). Notice 
that the permeability of the Inductor core has yet to be defined. The 
reason for this seeming remission Is that It Is really not an Independent 
variable; 


Perraeabi 1 Ity 


L*Z 


N^*A 


(5-66) 
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Concisely stated, for a given set of VI and F, the purpose of this 
optimization is to Identify all variables N, A, AC, and Z such that eqs. 

(5-41 ) and (5-43)can be satisfied, and concurrently the 06JF of eq. (5*65) 
can be minimized. The Identification of these variables numerically 
defines the particular power loss PL per(5-4S)that will result In a 
minimum emergency power system weight for a given set of VI and F. The 
calculation can be repeated for other sets of VI and F, from which the 
global minimum for one particular set of VI and F can be numerically 
obtained along with the specific PL corresponding to this set. Substit- 
uting these values 1nto(S-52)and(5-53)w111 then yield the optimum bat- 
tery voltage level VB and the precise minimum battery ampere-hour rating 
AH. 


For simplification purposes, let 

KA » CONST 1 + K3 * CONST 2 (5 

KB - K1 * CONST 2 (5 

KC » 4 * FC *0C fS 

KD = 01 + K2 * CONST 2 

KE - L * IP / BOC (5 

KF * (ff* FW)’°‘® (5 


KG » ]/{2*K ) 


(5 




KH ■ 0.5 


(5-74) 


X ■ . y ■ N®*^, V • (AC)^’®, £ • I (5-75) 

Then, equat1on(5-65)becomes; 


WT ■ KA + KB * + KC * xy^v^ + KD * x^z (5-76) 

Equat1ons(6-41 )and(5~43)become: 

KF * yv - KG * 2 + KH *x < 0 (5-77) 

- KE > 0 (5-78) 

Using the method of Lagrange Multipliers, the optimization function h 
(x, y, z, v) becomes; 

h(x,y,z.v)-KA-«-KjJ*y^xv'^+KC*xyV+KD-x^z-o((x^y^-KE)- ^(KF*yv-KG*z+KH*x) 
where a and Bare the Lagrange Multipliers. 


ik 

»x 


KB*y^v'^ + KC*y^v^ + 2*KD*xz -2*®^ *xy^ 


4*KH - 0 


(5-79) 


■ 2*KB’*yxv'^ + 2*KC*xyv^ - 2*J<*x^y - /^*KF*v ■ 0 


(5-80) 


- 306 - 


»h 

dz 


(5-81) 


■ KD*x^ ♦ /« *KG ■ 0 


• -2*KB*yxv’^ + 2*KC*xyv - yj*0 


(5-82) 


The six equations, (5-77) to (5-82), can be used to solve for tjhe six variables, 
X. y, z, V, o'arid a: 


N - ( 


w *L*IP*FW 
"55TO — 


0.5 -1 
) *S 


(5-83) 


A ■ 


/ L*IP*AC 


0.5 

) * S 


(5-84) 


z . 2.I.( (4^* S-'’-5) 


(5-85) 


At* 0t75 pnr ^ »25 *0*25 ^0*5 ^ c f5«-86) 

u • 2 * • (jfTre ) * <TT> *‘- *s * 
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a - 12*FC*DC*(4*FC*0C - ) 


b • C0NST2*K1* - 24*FC*DC) 


c ■ 3 * (K1 * CONST 2)^ 


Values thus obtained for N, A, Z, and AC can be used In eq(5-45) to 
numerically determine the power loss PL. 





5.11. APPROACH FOR OVERALL OPTIMIZATION 


At this juncture It Is perhaps worthwhile to review briefly what 
has been accomplished. First, one has to realize that the optimization 
study so far has been predicated by a given set of VI and F. Based on 
such a given set. the progression Includes the following: 

• The numerical Identification of Inductance L 
through eq. (5-42). 

e The numerical Identification of capacitance C 
through eq. (5-44). 

e The detailed design for Inductance L through 
eqs. (5-83) to (5-91). Details Include N, A, Z, 
and AC. 

f The determination of loss PL through eqs. (5-45) 
to (5-49). Such a PL will give a minimum total 
system weight. 

f The battery voltage VB can be found from eq. 

(5-52). 

e The battery ampere-hour rating AH can be found 
from eq. (5-53) . 

f The minimum total system weight Is prescribed 
by eq.(5-64) . 

It Is Iterated that the m1n1mu..« weight of (5-64) applies only to the given 
set of VI and F. Different minimum weights corresponding to different sets 
of VI and F can be similarly generated. By plotting this weight as a 
function of F, with VI as the varying parameter, a family of curves can 
be displayed to Identify the global minimum for one specific set of VI and 


F. Tht global minimum roprtstnts tht optimum dtsign for all betttry 
and convtrttr paramctars, and the corresponding battery voltage level 
can be determined from the specific VI. 

.12. COMPUTER PROGRAM FOR DESIGN OPTIMIZATION CALCULATION 

The foregoing design optimization equations are transcribed Into a 
computer program for numerical processing. The program Is given as 
Appendix K. A total of 28 parameters are given to the program In 
statements 2200 to 2230. Their corresponding numerical values are 
provided In statements 2000 to 2030. For a cost-effective program,no tol 
erances are assumed for these parameters as was Indicated In statements 
2100 to 2130. The program Is designated as "V^TOL" In 2300. The 
output data of the design optimization program, as specified In 2400. 
Include the following: 

L : Inductance in henries 

C capacitance In farads 

RC : equivalent series resistance of C In ohms 

IP peak Inductor current In amperef 

N : nimber of turns of L 

2 

A : cross-section area of the core for L In meter 

Z mean core length In meter 

U : core permeability In gauss/oersted 

2 

AC winding area per turn In meter 
PL total converter loss In watts 

PLC : copper loss of Inductor In watts 

PL I : Iron loss of Inductor In watts 

PM total Inductor loss In watts 
PT total loss of transistor switch In watts 
PD total diode loss In watts 

PC total capacitor loss In watts 
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P : Othtr convtrUr lossts In watts 
NB nuRibtr of colls In strlos 

VB minimum bittory voUogt notr tnd of dischorgo 
In volts 

UM Inductor wolght In grams 

WP : convortor packaging weight In grams (excluding 
power components) 

WC : converter cooqionent and packaging wolght In grams 
WB : battery weight In grams 

WCON : converter weight In gramr* 

WT total power system weight In grams 
AH ampere-hour capacity required of the battery 

From statements 3000 to 3064, the program essentially reproduces the 
various equations presented In the text, with the following correspond- 
ences : 


Statement In Program 
3000 
3002 
3004 
3006 
3008 
3010 
3012 

3014, 3016, 3018 


Equation In Text 
5«>42 
5-44 
RC*C-G 
5-49 
5-46 
5-43 
5-47 
5-48 
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3020 

CONST 1 Of (5-64) 

3022 

CONST 2 of (5-64) 

3024 

5-89 

3026 

5-90 

3028 

5-91 

3030 

5-87 

3032 

5-88 

3034 

S-39 

3036 

5-83 

3038 

5-84 

3040 

5-85 

3042 

5-66 

3044 

5-45 

3045 

5-17 

30^6 

5-50, 5-56 

:047 

5-19 

3048 

5-52 

3049 

5-18,5-19 5-20, 5-21 

3050 

5-58 

3052 

5-59 

3053 

5-65 

3054 

5-23 

3056. 3058 

5-29 

mo 

5-32 

3062 

5-33 


5-22 
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3064 a-34 

3066*3200 Printout sptcl float Ions 


Tht nuRitrIca) Inputs to tha co^)uttr ar^ 'll van as follows: 


V2 

AO 

f 

At 

VO 

A2 

PO 

A3 

T 

44 

0 

AS 

iac 

A4 

PI 

A7 

vs 

A8 

VIE 

A9 

VI 


PI 

It 

TS8 

12 

UF 

13 


ISO 0 

8000 0 

270 0 

3000 0 

•3£«00 0 

. 1 081*04 0 

4 0 

4 0 

t.8 0 

1.2 0 

0 

3.1414 0 

.3S€*04 0 

.3SC-04 0 


T5R 

P4 

TOP 

33 

RHO 

(6 

rc 

17 

rV 

sa 

BDC 

B7 

K 

CO 

N 

Cl 

OK 

C2 

KCAP 

C3 

K.1EC 

C4 

PC 

C3 

31 

C4 

KTME 

C7 


.15E-0# !) 

.!S£*04 • 0 

.I728E-07 0 

2 0 

.3SEOO 3 

.3SE«00 0 

42 0 

.9E<00 0 

t 0 

3 0 

S 0 

8700000 0 

7300000 0 

10 0 


As previously outllnad. a sat of VI and F will be assigned for each 
computer run, and the corresponding printout represents the optimum-weight 
design for that particular set of VI and F. 
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5,13. COMPUTER PRINTOUTS 


A sample printout Is shown In Figure 5-5. The aforementioned twenty- six 
parameters are printed under the heading of “X-AXIS", and their numerical 
values under "OUTPUT". Since tolerances have not been given to the pro- 
gram, the number "0* will appear under TOL" and "- TOL", and the lower 
limit LIM" and upper limit "+ LIM" will exhibit Identical numerical 
values as those displayed under "OUTPUT". 

A summary of arguments and tolerances accompanies each printout. For 
example. In the sample printout the output voltage VO Is designated A2, 
and has a specified value of 270V with zero tolerance. The particular 
run Is based on an Input voltage VI of 150V, switching frequency F of 10 
kHz, output power PO of 3000W, and power utilization time T of 0.5 hour. 

(See underlined portion of Figure 5-5). The total optimum system weight, 

WT, Is calculated to be ]03.9 kilograms. 

Computer printouts for other sets of VI, F, PO, and F. 

VI • 25, 50, 100, 150, 200 V 
F - 5, 10, 20, 30, 40, SO, 100 kHz for each VI 
PO - 3000U 
T » 0.5 hr. 

The calculated results for total system weight WT are presented in Figure 5-6. 
A family of curves showing weight versus frequency, with the Input voltage 
as the varying Index, are plotted. These results will be discussed In the 
next section. 
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5-14. DISCUSSION OF CALCULATED RESULTS 


The following observations can be made based on Figure 5-6: 

1. The weight vs. frequency curves generally exhibit U- 
shape characteristics, with the shape more pronounced 
at lower Input voltages. Within quite a wide fre- 
quency range, from 20 to 60 kHz, the total system 
weight Is relatively constant for a given Input voltage. 

2. For each frequency, the weight reduces with an Increase 
of Input voltage. The weight reduction is most pro- 
nounced for voltage Increments at lower Vi's. The re- 
duction becomes diminished for voltage increment at 
higher Vi's. For example, at 30 kHz, a 16 kg reduction 
Is realized when VI is increased from 25 to 50V. The 
corresponding reduction Is only 4.5 kg for a 100-to-200V 
increase. 

These observations are discussed as follows: 

THE U-SHAPE WEIGHT VERSUS FREQUENCY 

The U-shape weight versus frequency can be understood by plotting 
the corresponding loss versus frequency and weight distribution 
characteristics for all power components Including the Inductor, 
capacitor, transistor, and diode, which are part of the computer 
printouts. These characteristics are shown In Figures 5-7 and 5-8. 
It is clear from Figure 5-7 the loss versus frequency curve also 
exhibits the U-shape. The diminishing total loss at lower fre- 
quencies fs caused by the copper loss of the inductor, which 
decreases as frequency Increases due to the smaller Inductor 
size for a given source-current ripple requirement. The Increasing 


- 317 - 















total loss at higher frequencies Is caused by the higher switching 
losses associated with the transistor and diode switches. The op- 
posing loss trends for the Inductor and the semiconductor switches 
as a function of frequency leave a minimum total los;. In the 10 - 
30 kHz frequency range. The weight profiles shown In Figure 5*8 Is 
directly Influenced by the corresponding losses. A monotonically 
reducing Inductor weight as a function of frequency Is argumented 
by Increasing weights from both the battery and the converter pack- 
aging, thus explaining the aforementioned U-shape. Since the In- 
ductor losses and the semiconductor switching losses are more pro- 
nounced when the Input voltage Is low and the corresponding rms 
and peak current in the Inductor as well as the semiconductors are 
high, the U-shape Is also more pronounced at lower Input voltages. 

The opposing trends of weight versus frequency for the Inductor and 
the battery plus converter packaging leave a relatively constant 
total system weight In the middle frequency range, from 20 to 60 
kHz. 

WEIGHT REDUCTION WITH INPUT VOLTAGE 

The converter loss and battery internal loss are much higher at 
lower Input voltages. Consequently, greater weight reduction can 
be realized in all system components Including magnetics, capacitors, 
battery, and converter packaging when a given voltage Increment is 
added to a lower input voltage. As the loss-related system weight 
diminishes, the battery and converter packaging weights, which 
are constant functions of the output power, become dominant in their 
weight contributions to the system. Consequently, little weight 
saving can be realized by increasing VI when the efficiency gain is 
no longer an important factor. This fact is amply substantiated by 
curves exhibited in Figure 5-6, where the weight reduction from VI“ 
lOOV to VI-200V is rather meager. 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


IMPACT OF DIFFERENT .CORE MATERIAL ON TOTAL SYSTEM WEIGHT 

In the analysis presented so far. the Inductor core Is assumed to be 

made of the molypermalloy powder material to achieve a minimum core 

loss, with an operating flux-density limit of 0.35 weber/meter^. 

However, In view of the high ration of copper loss PLC to core loss 

PLI, It becomes apparent that a different core material with a higher 

operating flux density level will result In saving of total system 

weight even though the attendant core loss will be higher. A Cut-C 

core of Deltamax (orthonol) laminations, with a conservative flux- 

2 

density limit of 1.2 wever/meter Is selected to replace the power 
core In the computer run to assess Its Impact on the total system 
weight. This change Is accomplished In the computer program through 
the following program edition: 


. ii,2G30 

2030 BATA 1.2.42,.?. 1 . J.5.3.9E6.;.8£b. 1 0 

■ i£,30t2 

30i2 ro=. 1A5*A1 .o0*\A2-A0)»A0/y3/A2/ Y5 
' 3012 Yi=4.84»A1".2(4*(A2-A0)»A0/r3/A:/lf5 
< E.J04? 

3049 XI=2*X04.089*Ar .6*Z9^K0»A1 
j0r49 XI=2*X0*2.40-*Ai .26*29>K0M1 


The flux density BDC In statement 2030 Is changed from 0.35 to 1.20. 

The core-loss description In statements 3012 and 3049 are also changed 

to reflect the core-loss profile of the new material, fvuns were executed. 

A weight versus frequency plot Is given In Figure 5-9. Comparing the 

curve for a given VI In Figure 5-9 and Figure 5-6, It Is clear that a 10 

kg weight saving Is accomplished for VI-25V when BDC Is Increased from 
2 2 

0.35 w/m to 1.2 w/m that the U-shape of weight versus frequency Is not 

as pronounced at the low frequency end of Figure 5-9 as that of Figure 

5-6. For a given VI, the Inductor weight saving as a result of a higher 

frequency Is less In Figure 5-9 In relation to that 1n Figure 5-6. 
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5-16. Justification For All Input Parameters 

Justifications for numerical values of all Input parameters 
to the computer are given In Appendix K. These parametrrs 
are defined In the previous sections. 


1 


5-17. CONCLUSION 

From analysis presented In this report, the following design-related 
conclusions can be made regarding the VSTOL emergency power system: 

1) As expected, the total system weight reduces with the 
converter Input voltage. However, a level of diminish- 
ing return Is soon reached for Vlfc lOOV. 

2) For a given watt-hour rating, connecting more cells In 
series to effect a higher battery voltage Is more costly 
than connecting fewer series cells. In view of the 
significant weight saving that can be realized by raising 
VI from 25V to 30V, and the rather meager weight saving 
that Is realizable by using a hIgher-than-lOOV Input 
voltage VI, It Is reconmended that the battery voltage 

VB be set within a 50-1 OOV range. 

3) Depending on the required output power PO, which deter- 
mines the power-switch current for a given VI, a 25V 
battery may be used In smaller power applications, e.g., 

PO i lOOO iv. 
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4) With the present availability of power switches and 
magnetic core materials, the total system weight WT 
for a given VI stays nearly constant within a wide 
frequency range. The range covers from 20 to 60 kHz 

In Figure 5-6 and 15 to 40 kHz In Figure 5-9. A recommend- 
able operating frequency Is In the vicinity of 30 kHz. 

2 

5) For a core with B0C“1.2 weber/meter , the optimum 
system weight In kilograms can be read from F1y;ire 
5-9 as follows for P0«3 kw and T*0.5 Hr: 

Switching Frequency (kHz) Weight (k g) 



VI -25 V 

WV 

100V 

150V 

5 

127kg 

113kg 

105kg 

102kg 

10 

121 

109 

103 

101 

20 

118 

108 

103 

101 

30 

119 

108 

103 

101 

40 

121 

109 

103 

101 

50 

123 

no 

103 

101 

100 

i:., 

115 

106 

103 


6) The use of higher flux-density core f.4ter1al reduces the In- 
ductor copper loss by an amount much greater than the atten- 
dant Increase In core loss, thus reducing the total system 
weight. 

7) While the Impact of core material Is not to be overlooked, 
the major weight contributor In the system Is the battery. 
Except at very low frequencies (5 kHz), the battery gen- 
erally Is responsible for more than 80* of the total sys- 
tem weight. 
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6. CONCLUSIONS AND RECOHMLNOATIONS 


The Modeling and Analysis of Power Processing Systems (MAPPS) 
Project has developed the mathematical models and computer software 
program to perform the calculation for Buck* Boost* and Buck*Boost 
1 DC-OC Converters In order to determine: 

• Performance Analysis 
e Design Optimization 

A computer-aided discrete time domain modeling and analysis 
technique for Performance Analysis has been presented which Is 
applicable to all types of switching regulators using any types of 
duty cycle controllers and operating with continuous, as well as 
discontinuous. Inductor current. State space techniques are employed 
to characterize coverters exactly by the nonlinear discrete time 
domain equations In vector forms. Newton's Iteration method Is 
employed to solve for the exact equilibrium state of the converter. 
The system Is then linearized about Its equilibrium state to arrive 
at a linear discrete time model. The stability nature and transient 
responses are studied by examining the eigenvalues of the linear 
system. Changes In eigenvalues due to system parameter changes can 
be plotted In the complex z-plane yielding an excellent design tool 
very similar to conventional root-locus plots. The analysis Is also 
' extended to determining the frequency related performance character- 

istics such as the closed loop Input-to-output transfer function used 
to determine the audlosusceptiblllty of the converter. The modeling 
^ and analysis approach makes extensive use of the digital computer as 

an analytical tool, replacing highly complex and tedious analyses by 
numerical method and making automation In power converter design and 
analysis possible. 
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6. CONCLUSIONS AND RECOMMLNDATIONS (Cont'd) 


In addition to Its particular utility at analyzing high-frequency 
control-loop related phenomena, the analysis also serves as a useful 
design tool which provides design guidelines for such Important control 
parameters as the dc loop gain, the ac loop gain, and the R-C compen- 
sation network of a two loop converter to optimize its transient 
response and to stabilize the system. 

To those working with switching regulators, converters, and systems 
comprised of these equipment, certain design and analysis Intricacies 
invariably make themselves felt throughout the equipment and system 
design and development stages. Empirical and Intuitive reliances often 
Intercede with the designer's desire to be "more scientific" and his 
commitment of being "on schedule". Handicapped by a general lack of 
established modling, analysis, design, and optimization tools. It has 
not been uncommon for a power processing designer to fulfill very 
little of the desire and/or the commitment. 

The cost and schedule plights that most equipment and system 
designers find themselves In have to do with at least one of the 
following entitles: power circuit weight/efficiency, control- 

related performance requirements, and trial -and-error power and 
control design Iterations. While power processing as a technology 
has reached the level of sophistication where the modeling, analysis, 
design, and optimization of these entitles should have been established, 
a survoy of literatures conducted at the Initiation of the MAPPS program 
had proved the contrary. In addition, the recent evolving trend of 
higher power and equipment standardization has further heightened the 
need for analytically-based design and optimization. 

In this regard, the program has accomplished the following objectives 

• The methodologies of power processing model 1ni|, analysis, 
design, and optimisation, are all established. 
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6. CONCLUSIONS AND RECOMMENDATIONS (Cont'd) 


• Application-oriented analysis, design, and optimization sub- 
programs for power-circuit design, control circuit design, 
and control performance analysis are becoming available. 

e Cost-effective system configuration study and system disturbance 
propagation are now feasible. 

Being government sponsored, all softwares ure available without pro- 
prietary complications. 

Continued MAPPS effort will aim at the following goals: 

• Analyze perTormance for commonly-used power processing equipment 
and selected systems. 

• Detailed power circuit design optimization to meet given power- 
related performance requirements for most commonly-used power 
circuit configurations. 

• Standardize control-circuit design to meet control -related 
performance requirements. 

• Provide cost-effective tools for the identification of optimum 
system configurations and system failure-mode analysis. 

The following is a list of basic tasks for future work: 

• Application Handbook for Performance Analysis. 

• Application Handbook for Design Optimization. 

• Continued Current Injection Multi loop Control Modeling. 

(a) Single Power Stage. 

(b) Parallel Modular Power Stages. 

• Continued Power Subsystem Optimization Techniques. 

• Development of Performance and Design Optimization for 
Series Resonant Inverter Power. 
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